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Abstract 

The state of the art concerning transverse-spin and transverse-momentum phenomena in hard 
hadronic reactions is reviewed. An account is given of single-spin and azimuthal asymmetries in 
semiinclusive deep inelastic scattering, e^e~ annihilation, Drell-Yan production, and hadroproduc- 
tion. The ongoing experiments and the main theoretical frameworks are described in the first part 
of the paper. The second part is devoted to the experimental findings and their phenomenological 
interpretations. A brief discussion of the perspectives of future measurements is finally presented. 
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1 Introduction 



For quite a long time the common lore in the hadron physics community has been that transverse 
polarisation effects are negligibly small in hard processes. In the last two decades a growing theoretical 
and experimental evidence has shown that this is not the case and that transverse-spin phenomena are, 
on the contrary, rather relevant in various high-energy hadronic reactions. 

The prehistory of the subject started in the mid-70s, when substantial single-spin asymmetries 
(SSA's) were found in inclusive pion hadroproduction at the center-of-mass energies of the Argonne 
synchrotron (few GeV) [H [21 [3]. At the same time, at Fermilab A hyperons produced in unpolarised 
pN collisions at y/s ~ 24 GeV and moderate transverse momenta Pt (below 1.5 GeV) were found to 
possess a large transverse polarisation [3], a result subsequently confirmed at slightly higher -y/s and Pt 
0. 

These findings stimulated both experimental and theoretical work. An experimental programme 
to investigate both longitudinal and transverse spin effects in high energy pp and pp scattering was 
proposed in 1978 at FNAL |6] and carried out more than 10 yeas later by the E704 Collaboration. 
Measuring inclusive pion production in collisions of transversely polarised proton and antiproton beams 
with an hydrogen target at the center-of-mass energy -y/i = 19.4 GeV and for P^ up to 2 GeV, the 
E704 Collaboration found single-spin asymmetries as large as 40 % in the forward region [3 [HI [HI [10] . 
More recently, the experimental collaborations STAR, PHENIX and BRAHMS, working at RHIC, 
have confirmed the early Fermilab findings on single-spin asymmetries in hadroproduction, pushing the 
frontier of the cm. energy to y/s = 200 GeV and covering wider kinematical ranges in P^ and in the 
Feynman variable a;^ [HI [El [131 [11 [15] . 

On the theoretical side, soon after the first experimental findings of transverse spin effects, Kane, 
Pumphn and Repko proved in a famous paper [16] that in coUinear perturbative QCD (applicable to 
high Pt) SSA's are of the order of as{mg/ ^/s) (where ruq is the quark mass) and therefore vanish in the 
massless limit. Also, in other important theoretical works [T71 [181 [IH] it was shown that non- vanishing 
transverse single-spin asymmetries may arise in QCD only if one consider higher-twist contributions 
(naively expected to behave as a power of (M/P^), where M is a hadronic scale. 

It took a while to realise that theory allows for transverse polarisation effects which are in some 
cases unsuppressed (for a review see f.i. Ref. [20]). In the early 90s various authors [211 [221 [23 [21] 
rediscovered the distribution of transversely polarised quarks in a transversely polarised nucleon first 
introduced by Ralston and Soper in 1979 [25] . This "transversity" distribution, usually denoted by 
hi{x) or by Axqix), is a leading-twist quantity that contributes dominantly to the double transverse 
asymmetry in Drell-Yan (DY) production. Due to its chiral-odd nature hi is not measurable in inclusive 
Deep Inelastic Scattering (DIS), where transverse SSA's are prohibited by time-reversal invariance at 
lowest order in [2S]- This argument, however, does not hold in semi-inclusive DIS (SIDIS), where 
at least one hadron in the final state is detected on top of the scattered lepton. In SIDIS processes 
no first principles forbid SSA's. Various theoretical proposals were soon put forward to measure hi 
[271 [281 [291 [301 [311 [321 [33] . In particular, Collins proposed a mechanism, based on a spin asymmetry 
in the fragmentation of transversely polarised quarks into an unpolarised hadron (the "Collins effect"), 
which involves a transverse-momentum dependent (TMD) fragmentation function, H^. This mechanism 
was originally proposed as a "quark polarimeter" , and could be conveniently exploited to measure the 
transversity function hi{x) in SIDIS. 

In a different approach, one year before the publication of the E704 results, Sivers had suggested 
that single-spin asymmetries could originate, at leading twist, from the intrinsic motion of quarks in 
the colliding hadrons [Ml [35]. The idea, in particular, was that there exists an azimuthal asymmetry 
of unpolarised quarks in a transversely polarised hadron (the so-called "Sivers effect"), and a new T- 
odd TMD distribution function, now commonly called Sivers function and usually denoted by f^, was 
proposed to describe the partons in a transversely polarised nucleon. 



Originally this mechanism seemed to violate time-reversal (T) invariance and it was demon- 
strated that had to be zero. Brodsky, Hwang and Schmidt [36| [37] proved however by an explicit 
calculation that final-state interactions in SIDIS, arising from gluon exchange between the struck quark 
and the nucleon remnant, or initial-state interactions in DY, produce a non-zero Sivers asymmetry. 
The situation was further clarified by Collins [38] who pointed out that, taking correctly into account 
the gauge links in the TMD distributions, time- reversal invariance does not imply a vanishing f^, but 
rather a sign difference between the Sivers distribution measured in SIDIS and the same distribution 
measured in DY. 

The phenomenological analysis of the E704 results |39l HQ] , however, showed that both the Collins 
and the Sivers effects are at work to generate the observed asymmetries, but a satisfactory theoretical 
description of the data is still missing today. 

New experimental opportunities came out in the 90's, when it was realised that high energy SIDIS 
experiments were needed to investigate the helicity structure of the nucleons. In fact a major event 
had focused on the nucleon spin the attention of the high energy physics community. In 1988 the EMC 
Collaboration at CERN, scattering a high energy polarised muon beam on a transversely polarised 
proton target obtained a totally unexpected result, namely that the fraction of the nucleon spin carried 
by the quarks was small, even compatible with zero, within the accuracy of the measurement. From the 
inclusive cross-section difference for parallel and antiparallel spins one could extract a linear combination 
of the quark helicity distributions Aq (or ^fi), defined as the difference of the quark densities for quark 
spin parallel or antiparallel to the longitudinal nucleon spin. Using complementary information on the 
quark helicities derived from the weak decays of the hyperons it was possible to add up the quark 
helicities, thus obtaining AE, the overall quark contribution to the nucleon spin. The result, which 
came to be known as the "spin crisis" was at variance with the current paradigm, i.e. the quark model 
and the beautifully simple explanation of the baryon magnetic moments. More than one thousand 
theoretical papers were written on the subject, and many experiment (SMC at CERN, E142, E143, 
E154 and E155 at SLAC) were proposed and executed to confirm the effect, to extend the result to the 
neutron, and to improve the accuracy of the measurement. The confirmation of this finding led to a 
growing attention to the other contributions to the proton spin, namely the gluon polarisation and the 
orbital angular momentum of both quarks and gluons, as well as to a deeper look at the QCD description 
of the nucleon and to the relation between hi and gi. Thus, new generation experiments, well suited 
to investigate SIDIS with both longitudinally and transversely polarised targets, like COMPASS and 
HERMES, were proposed and started their operations about 10 years ago. 

In 2004 HERMES [31] and COMPASS US US] presented the first data collected with transversely 
polarised proton and deuteron targets, which showed clear evidence of transverse SSA's on proton. One 
of the main advantages of SIDIS is that the Collins and Sivers effects, as well as the other TMD effects, 
are not mixed, as in hadropro duct ion, but generate different azimuthal asymmetries, which can be 
separately explored. Thus, the Collins and Sivers asymmetries could be extracted analysing the same 
data. 

Another major step in the understanding of the Collins effect occurred from the Belle Collaboration 
studies of the azimuthal correlation between the hadrons in the two jets created in e~^e~ annihilations 
[H]. In the process e~^e~ — > qq the transverse polarisations of the qq pair are correlated, thus the 
Collins effect is expected to cause correlated azimuthal modulations of the hadrons into which the q 
and the q fragment. The high precision of the Belle data provided very accurate measurements of such 
modulations, and a combined analysis has allowed a first extraction of both the Collins function and of 
the transversity distribution [15]. The Sivers and Collins effects are by now theoretically well established 
and the overall picture is essentially in agreement with the still limited set of results produced by the 
SIDIS experiments. 

All this work on transverse polarisation effects eventually opened up the Pandora box of the 
transverse-momentum structure of hadrons. The importance of the intrinsic transverse momentum 



of quarks in hadrons has been acknowledged since many years. The transverse momentum of the 
quarks is responsible for the large azimuthal asymmetries of the hadrons produced in SIDIS processes 
on unpolarised nucleons (the so-called Cahn effect). In a similar way it is largely responsible for the 
azimuthal asymmetries observed in DY processes, namely in the production of a lepton pair in hadron- 
hadron scattering at high energy. When the intrinsic transverse momentum of the quarks in the nucleon 
is taken into account, several new functions are needed to describe the transverse spin structure of the 
nucleon. Transverse spin, in fact, couples naturally to the intrinsic transverse momentum of quarks, 
and the resulting correlations are expressed by various transverse-momentum dependent distribution 
and fragmentation functions, that give rise to a large number of possible single-spin and azimuthal 
asymmetries [311 EHl [291 SHI 113 HHl HH]. Of particular interest are the correlations between the quark 
transverse momentum and the nucleon spin, the quark spin and the transverse momentum of the frag- 
menting hadron, and the quark transverse spin and its transverse momentum in an unpolarised nucleon, 
which give rise to the Sivers function, the Collins function and the so-called Boer-Mulders function re- 
spectively. All these three functions are (naively) T-odd, and all three are responsible for transverse 
spin asymmetries in SIDIS. In particular the Boer-Mulders function [50], measures the transverse-spin 
asymmetries of quarks inside an unpolarised hadron, and contributes to the cos (j) and cos 2(f) azimuthal 
modulations in the cross sections of unpolarised SIDIS and DY processes which have been observed 
since many years and are presently been accurately measured. 

The TMD description of hard processes has been put on a firm basis by the proof of a non-collinear 
factorisation theorem for SIDIS and DY, in the low transverse momentum regime [511 152]. On the other 
hand, it is known that twist-3 coUinear effects, expressed by quark-gluon correlation functions, can also 
produce single-spin and azimuthal asymmetries [53l |54l [55] . This mechanism works at high transverse 
momenta, Pt 3> M. Thus, there is an overlap region where both the collinear twist-3 factorisation 
and the non-collinear factorisation should be both valid. The relation between these two pictures, that 
is, between the T-odd TMD functions on one side and the multiparton correlators on the other side, 
has been clarified in a series of recent papers [561 [SZl [58]. These works have opened the way to the 
derivation of the evolution equations for the TMD functions, a longstanding problem in transverse-spin 
and transverse- momentum physics [SHI [SOI [SD [S2] • 

Before concluding these introductory remarks, a caveat is in order. This review is far from being 
exhaustive. Transverse spin physics is in fact developing so fast that it is nearly impossible to cover all 
the results and the ongoing work. The following pages necessarily reflect the specific competence and the 
preference of the authors. Among the various processes involving the transverse spin and the transverse 
momentum structure of hadrons, we chose to focus on SIDIS. From a theoretical viewpoint these are the 
cleanest and best understood reactions. A related important process, that we will also treat, is hadron 
pair production in e^e~ annihilation, which probes transverse-spin fragmentation functions. As we 
mentioned, a great wealth of data on transverse-spin phenomena come from inclusive hadroproduction. 
We will pay less attention to these processes, because a recent review [63] is largely dedicated to their 
phenomenology. Generalised Parton Distributions will also be only mentioned since they are nicely 
covered in a very recent and comprehensive review [61] . 

A comment on the notation is in order. The proliferation of distribution and fragmentation functions 
involved in transverse-spin and transverse-momentum phenomena makes the issue of notation and 
terminology a very problematic one. Throughout this paper we adopt for the distribution functions 
the Jaffe-Ji nomenclature [23], extended to transverse momentum dependent distributions by Mulders 
and collaborators [IHl [5D] illustrated in detail in Section [3751 . Thus, /i(a;), gi{x), and hi{x) are the 
unpolarised, the helicity and the transversity distribution functions, respectively, with the subscript 
1 denoting leading-twist quantities. The main disadvantage of this nomenclature is the use of gi to 
denote a quark distribution function whereas the same notation is universally adopted for one of the 
two structure functions of polarised deep inelastic scattering. Other common names in the literature 
are q{x) for the unpolarised distribution, Ag(x) for the helicity distribution, ATq{x) for the transversity 



distribution (which is also called sometimes 5q: here we reserve this name to the tensor charge). The 
fragmentation functions are denoted by capital letters: D (unpolarised), G (longitudinally polarised), 
H (transversely polarised). Thus, Di is be the usual leading- twist unpolarised fragmentation function, 
Gi the fragmentation function of longitudinally polarised quarks. Hi the fragmentation function of 
transversely polarised quarks. Note that capital letters are also used for the gluon distribution functions 
[65] and for the generalised parton distributions [66] and the reader should be aware of this possible 
source of confusion. 

The outline of this paper is the following. In Section 2 a brief account of the technical features of the 
main ongoing experiments is presented. Section 3 is devoted to the formal aspects of the transverse-spin 
and transverse-momentum structure of hadrons. In Section 4 we introduce the relevant processes and 
observables, and to the theoretical frameworks that describe them. Section 5 illustrates the experimental 
findings about single-spin and azimuthal asymmetries, and their phenomenological interpretations. In 
Section 6 we discuss the short- and mid-term perspectives of planned and proposed measurements. 
Finally, Section 7 contains some concluding remarks. 

2 The experiments 
2.1 DIS experiments 

Deep inelastic scattering as a tool to unveil the nucleon structure was invented in the late 60 's at 
SLAG, when for the first time a high energy (1 GeV) electron accelerator became available, a wealth 
of e-N scattering data were collected, and eventually it became clear that scattering at large transverse 
momentum could be interpreted as elastic scattering off the nucleon constituents, the "partons" . From 
the dependence of the cross-section on the energy and the momentum transfered to the nucleon it has 
been possible to identify the charged partons with the quarks, and assess the existence of the gluons, 
as carriers of half of the proton momentum. In these experiments, only the scattered electrons were 
detected with suitable magnetic spectrometers, and no coincidence experiments were possible, due to 
the small duty cycle of the intense electron beam. In the subsequent years the Linear Accelerator 
energy was gradually increased, to reach eventually 50 GeV in the most recent experiments, allowing 
to measure at larger and larger Q^, and at smaller and smaller x. Higher energy experiments could 
be performed at GERN and at FNAL by constructing muon beams from vr and K decays: all these 
data and the neutrino-Nucleon data eventually led to the extraction of the parton distribution functions 
(PDF's) and of their dependence. 

In a second generation of experiments polarised lepton beams and polarised targets have been used. 
Sophisticated techniques have been developed to polarise the electron beam at SLAG: in the latest 
experiments (for instance E155), 85% beam polarization was typically achieved by using as source 
the photoelectrons emitted by a gallium arsenide surface. The high energy muon beams on the other 
hand were naturally polarised in the weak decay process: by suitably choosing the ratio between muon 
and pion momentum (0.94 in the EMG experiment at the GERN SPS) muon polarisation of ~80% 
are obtained. The goal of these experiments was the measurement of the structure functions gi and 
g2-, and to verify the Bjorken sum-rule. The beam was longitudinally polarised, while the targets were 
longitudinally polarised to measure gi and transversely polarised to measure • A variety of solid targets 
have been used, butanol, ammonia or ^LiD, which are kept at very low temperatures (< 0.1 K) and in 
a strong magnetic field (typically 2.5 T, but for some small targets even 5 T). In these targets a dopant 
is added, which provides unpaired electrons which in the high magnetic field and at low temperature 
get polarised to almost 100%. By irradiating the sample with microwave of proper frequency one 
can induce hyperfine transitions which flip the proton (or deuteron spin) to the preselected spin state 
(Dynamical Nuclear Polarisation, DNP, method). Only the free protons (deuterons) are polarised, so 



Figure 1: The HERMES polarized target. The atomic beam from the source ABS is focused in the 
sextupole system and diffuses in the storage cell (right panel). Nuclear polarization is obtained by 
inducing hyperfine transitions by RF-pumping. A Breit-Rabi polarimeter (BRP) measures the beam 
polarization. 

that the "figure-of-merit" for a spin-asymmetry measurement, which is ~ l/cr^^^j^, is proportional to 
pPfg^, where / is the ratio between the polarisable nucleons and the total number of nucleons in the 
target, and P^g^ is the free protons (deuterons) polarisation. These experiments allowed to discover the 
"spin crisis", to extract the helicity quark distributions (71, to provide first measurements of the quark 
contribution to the nucleon spin AS, and to verify the validity of the Bjorken sum-rule. 

A third generation of experiments aiming at the study of the nucleon spin structure started in the 
past decade. They still use polarised beams and polarised targets, but complement the detection of the 
scattered lepton with the reconstruction and identification of the hadrons produced in the fragmentation 
of the struck quark, the so-called current jet. A suitable trigger system still allows these experiments 
to record DIS events, but the main objective of the measurement is SIDIS events. These coincidence 
experiments require the disposal of a continuous beam and of a large acceptance spectrometer with full 
particle identification. 

2.2 SIDIS experiments 

At HERA the HERMES experiment was designed to utilise the circulating electron or positron beam. 
At the experiment, the stored beam (27.5 GeV and 40 mA) passes through a cell, a tube 60 cm long, 
coaxial with the beam, in which polarised atoms of hydrogen or deuterium are pumped in from an Atomic 
Beam Source. After diffusing in the storage cell, the atoms are pumped away by a huge pumping system 
before they diffuse into the electron beam pipe. Polarisation is achieved in the Atomic Beam Source by 
Stern-Gerlach filtering followed by radio-frequency transitions to the selected spin state. A schematic 
view of the target system is given in Fig. [1] This target system is particularly attractive when compared 
to the solid polarised targets because there is no dilution of the target polarisation due to the presence 
of the unpolarised nucleons bound in the other nuclei present in the material. Of course, the target 
thickness cannot be increased at will, not to destroy the circulating electron beam, but densities of 10^^ 
nucleons/cm^ were regularly achieved. After the target, a large acceptance magnetic spectrometer based 
on a 1.3 Tm dipole magnet analysed all charged particles up to 170 mrad in the horizontal plane and 
between 40 and 149 mrad in the vertical plane. The reduced acceptance in the vertical plane was due 
to the fact that since both the electron and the proton beam pipes of the HERA collider go through the 
middle plane of the magnet, its gap was divided into two identical sections by a horizontal septum plate 



Figure 2: Side view of the COMPASS polarised target. One can see the three target cells inside the 
mixing chamber, the microwave cavity, the solenoid coil, the correction coils, and the end compensation 
coil. 

that shields the electron and the proton beams from the dipole magnetic field. Charge particle tracking 
is provided by several micro-strip counters, multiwire proportional chambers, and drift chambers located 
before, inside and behind the magnet. Charged particle identification is provided by a RICH Cherenkov 
counter, while electron-hadron discrimination is achieved with a lead-glass calorimeter with a pre-shower 
hodoscope in front, and by a Transition Radiation detector. At the end of the spectrometer, a muon 
hodoscope located after an iron absorber helps the muon identification. The experiment took data with 
polarised targets until 2005. After an upgrade to implement the spectrometer with a recoil detector to 
investigate exclusive channels, it took data on unpolarised targets from 2008 to July 2009, when HERA 
ceased operation, and it has afterwards been dismantled. 

At CERN the COMPASS experiment has been assembled in the Hall 888, where the EMC and 
afterwards the SMC experiments were installed. Its physics program includes not only the investigation 
of the spin structure of the nucleon, but also the search of exotic light-quark hadronic states, like 
glueballs and hybrids. For spectroscopy the experiment uses hadronic beams (mostly pions and protons), 
and data have been collected with a liquid hydrogen target and various nuclear targets in 2008 and 
2009. In the following only the configuration which has been used for the study of the spin structure of 
the nucleon will be described, which uses the high energy muon beam at the CERN SPS. COMPASS 
has been designed to deal with a beam intensity five times larger than that of the previous EMC or 
SMC experiments. Typical intensities are 2 ■ 10^ muons per spill (about 5 s every 18 s) at 160 GeV, 
the momentum at which most of the data have been collected. Since there is no problem of beam 
radiation in the target, the target length has been chosen as long as possible, within the boundaries 
of the complexity and cost of the cryogenic system. The target materials which have been used in 
so far are ^LiD (/ ~ 0.4) as a deuteron target, and NH3 (/ ~ 0.15) as a proton target. The target 
system uses a solenoidal superconducting magnet, providing a highly homogeneous field of 2.5 T over 
a length of 130 cm along the axis, and is schematically shown in Fig. |5J About one kg of material is 
contained in a 4 cm diameter cylinder, coaxial with the beam, over a length of 120 cm, distributed either 



Figure 3: An artistic view of the polarised ^He gas target used in Hall A of Jlab for several experiments. 
For the E06-010 neutron transversity experiment important upgrades have been made, including the 
isntallation of a third pair of Helmoltz coil, enabling to orient the target spins in any direction. 

over two cells (in 2002, 2003, and 2004), or over three cells (since 2006). Nearby cells are oppositely 
polarised, so that scattering data on the two orientation of the target are taken simultaneously to 
minimise systematic effects. A cryogenic system allows to keep the target at temperatures of about 0.5 
K, and to polarise it with the DNP method. Once high enough polarisation values are reached, the R-F 
is switched off, the temperature drops to less than 50 mK, the spins get frozen, and data taking can 
start. A set of two saddle coils allows to get a transverse field of up to 0.6 T which can be either used to 
adiabatically rotate the target polarisation from parallel to antiparallel to the beam, or to set it in the 
transverse mode, orthogonal to the beam direction. At regular intervals, the polarisation orientation 
of the target cells are reversed by changing the frequency of the microwaves, so that possible effects 
due to the different acceptances of the different cells can be cancelled in the analysis. The experiment 
is still running, and in 2010 and 2011 will take more data with a NH3 polarised proton target, 50% 
of the time in the transverse polarisation mode and 50% of the time in the longitudinal mode. Large 
angular and momentum acceptance is guaranteed by a two-stage magnetic spectrometer, 60 m long, 
centred around two dipole magnets with 1 Tm and 4.4 Tm bending power respectively. A variety of 
tracking detectors ensures charged particles tracking from zero to ~200 mrad scattering angle, and 
charged particle identification is provided by a RICH Cherenkov counter. Two hadronic calorimeters, 
two electromagnetic calorimeters and two muon filters complement the particle identification and allow 
the reconstruction of neutral pions. 

At Jefferson Lab ( JLab) many measurements of electron-nucleon scattering and in particular of DIS 
have been performed over the past 10 years using the electron beam of CEBAF (Continuous Electron 
Beam Accelerator Facility), energies from 0.8 to 6 GeV, and polarised targets. In Hall A, the focus 
has been on measurements on a neutron target. Thanks to the very high beam current, a polarised 
^He gas target could be used as a neutron target. The advantage of this target is that to first order 
one can think that its spin is carried entirely by the neutron, since the two protons have their spin 
anti-aligned. The ^He gas fills a pressurised glass vessel (10 atm, typically) and is mixed with Rubidium 
vapour whose electrons can be polarised via optical pumping with circularly polarised laser light. The 
^He nuclei get polarised through spin exchange collisions with the Rubidium atoms. A schematic view 
of the target system is shown in Fig. [31 These targets can stand a much larger beam intensity than a 
solid target, without suffering of radiation damage. With a 15 /xA electron beam on a 40 cm target at 



10 atm a luminosity of 10^® cm~^s~^ is achieved. A series of high precision experiments have provided 
invaluable information on gi and g2 for the neutron particularly at large x values. The measurements 
have utilised a pair of high resolution magnetic spectrometers to measure the scattered electron, and 
by changing the angular settings and the momentum settings the structure functions could be precisely 
measured over a broad {Q^ — W) plane. On the other hand, no SIDIS measurements were possible, due 
to the small angular acceptance of the spectrometers. 

Complementary measurements, using solid polarised targets (both Ammonia, and Deuterated am- 
monia) have been carried on in Hall B, using the CEBAF Large Angle Spectrometer (CLAS). The 
large acceptance of this spectrometer has allowed to detect also hadrons to study SIDIS and exclusive 
events, but the target geometry did not allow to put the polarisation orthogonal to the beam, so that 
no transversity measurements have been possible in so far. 

Recently, in 2009, the first transversity measurements have been performed at JLab, by the E06-010 
collaboration, in Hall A. Another important feature of the ^He target is that the field necessary to 
hold the polarisation is low, so that it is easy with a set of three pairs of Helmoltz coils to rotate the 
polarisation in any direction, in particular to a direction orthogonal to the beam and measure SSA's on 
a transversely polarised target. The experiment has used the High Resolution Spectrometer to detect 
the hadrons in the SIDIS reaction (e — )■ e'vr^), and a new large acceptance (64 msr) spectrometer, 
the BigBite spectrometer, to measure the electrons. Thus the first ever measurement of SIDIS on 
a transversely polarised neutron target should soon be available, to complement the HERMES and 
COMPASS proton and deuteron data. 

The HERMES, COMPASS and JLab experiments have given an important contribution to precisely 
measuring AS, the quark contribution to the nucleon helicity (0.33 ± 0.01 ± 0.01), and to verify with 
some precision the Bjorken sum-rule. Moreover HERMES and COMPASS have provided first estimates 
of AG/G, the gluon contribution to the nucleon spin, from the spin asymmetry of the cross section 
of pairs of hadrons. An important result obtained by the CLAS collaboration is the first evidence for 
a non-zero beam-spin azimuthal asymmetry in the semi-inclusive production of positive pions in the 
DIS kinematical regime. This effect, observed also by HERMES, is not leading twist, and should give 
information on quark-gluon correlation. Most relevant to this report, HERMES and COMPASS have 
been the first experiments to measure SIDIS events on transversely polarised targets, and the data 
they have collected are still the only ones in this field. HERMES took SIDIS data from 2002 to 2005 
and new results from those data are coming and will still come. COMPASS took SIDIS data with 
the transversely polarised deuteron target in 2002, 2003 and 2004, and with the transversely polarised 
proton target in 2007, and more data will be collected in 2010. JLab has carried on measurements with 
the transversely polarised He3 target in 2009, but no data exist yet at the time of this report. 

2.3 Hadroproduction experiments 

A novel attack to transverse spin phenomena is provided by the Relativistic Heavy Ion Collider RHIC, 
at BNL, which has been designed to accelerate and store not only ion beams, but also polarised proton 
beams. As a high-energy polarised proton collider, RHIC by now is the world flagship facility for spin 
effects in hadron physics. The polarised proton beams are accelerated at the AGS to an energy of about 
25 GeV, then transferred into the RHIC rings and accelerated up to the desired energy, typically in the 
range 100- 250 GeV. An overview of the accelerator complex is shown in Fig. |H The polarisation in 
RHIC is maintained by two Siberian snakes, namely two sets of four helical dipole magnets which rotate 
the proton spin by 180°. The snakes are placed at two opposite points along the rings, so that the beam 
deflection between the two snakes is exactly 180° and the spin tune is thus equal to a half-integer and 
energy independent, thus cancelling the effect of both imperfection and intrinsic resonances. In this 
way polarisation values of about 70% have been obtained. Two more sets of the same helical dipole 
magnets before and after the two major collider experiments (PHENIX and STAR), for a total of 8 




10-25% Helical Partial Siberian Snake 



Figure 4: The Brookhaven hadron facility complex, which includes the AGS Booster, the AGS, and 
RHIC. Two snakes per ring and four spin rotators per each of the two large experiments (STAR and 
PHENIX) are also shown. 



spin "rotators" , allow to set the spin direction from vertical to parallel to the beam and then back to 
vertical, to allow the experiments to use either longitudinally polarised protons or transversely polarised 
protons. Polarisation-averaged cross-sections for pion production at y/s = 200 GeV have already been 
measured at RHIC both at mid-rapidity and in the forward region and found to be well described by 
next-to-leading order (NLO) pQCD calculations. 

The major players at RHIC are the three experiments BRAHMS, PHENIX and STAR. Since no spin 
rotators are installed before and after the BRAHMS experiment, this experiment always took data with 
transversely polarised proton beams. The experimental set-up consisted of two movable spectrometer 
arms to measure charged hadrons over a large rapidity and transverse momentum range. The forward 
spectrometer consists of 4 dipole magnets, providing a total bending power of 9.2 Tm, it can be rotated 
between 2.3 and 15 degrees, and has particle identification capability thanks to a RICH Cerenkov 
counter. The second spectrometer uses a single dipole magnet (1.2 Tm bending power), can be rotated 
from 34 to 90 degrees relative to the beam, and covers a solid angle of approximately 5 msr. The 
experiment completed data taking in 2006. The PHENIX experiment is one of the two large ongoing 
experiments investigating the proton spin structure at RHIC. The detector consists of two spectrometer 
arms at mid-rapidity [r] < 0.35) and two larger-rapidity spectrometer arms at 1.2 < 77 < 2.4. The 
mid-rapidity spectrometers identify and track charged particles, and are equipped with electromagnetic 
calorimeters. Two more electromagnetic calorimeters cover the large rapidity region, 3.1 < 77 < 3.7. 
The experiment has been designed to detect rare probes, and it has a sophisticated triggering system 
and a fast data acquisition system. The second large ongoing experiment is the STAR experiment. 
The apparatus has cylindrical symmetry around the mean of the two beams directions, and consists of 
a Solenoid magnet, a large Time Projection Chamber (TPC), coaxial with the magnet, and a barrel 
electromagnetic calorimeter just outside of the TPC. More electromagnetic calorimetry is provided in 



the forward region by an endcap calorimeter and by two pion detectors. Two more TPC's provide 
information on the charged particles in the angular range spanned by the forward calorimetry. 

In the longitudinal mode both PHENIX and STAR have concentrated their experimental programme 
on the study of double-spin asymmetries in the production of mesons and jets at high pt, to probe the 
gluon polarisation AG/G, and first results have already been published, which confirm the smallness of 
AG/G observed in the SIDIS processes. In the transverse spin mode the large single spin asymmetries 
in meson (pion) production in proton-proton scattering already observed in the late 1970's at lower 
energy at Argonne and at CERN, and in the 1990's by E704 at FNAL, have been confirmed to persist 
at the RHIC energies. For inclusive pion SSA's the overall accumulated statistics by now is so large 
that data can be binned in Pt, in rapidity and in xp- The behaviour of the SSA's as a function of the 
various kinematical variables is essential to constrain the parameters of the phenomenological models 
and assess the physical origin of the observed asymmetries. 

2.4 Electron-positron collider experiments 

Experiments at electron-positron colliders play a special role in the extraction of the transversity PDF's. 
It has been seen in the Introduction that the Collins conjecture, namely that in the hadronisation of 
a transversely polarised quark the hadrons of the jet might exhibit a left-right asymmetry relative 
to the plane defined by the quark momentum and the quark spin, might conveniently be exploited 
to measure the transversity distribution. To unfold the measured Collins asymmetry and extract the 
transversity distribution a knowledge of the Collins function is mandatory, but QCD tools are not 
capable to calculate either the quark distribution functions, or the fragmentation functions, and the 
Collins function is no exception. However, such an effect might be detected in high energy e~e~^ 
annihilations into two jets, and it was unambiguously observed by the Belle Collaboration |31], analysing 
data collected at the asymmetric e~e~^ KEKB storage rings, as will be explained in Section Jl]2l The Belle 
detector is a large-solid-angle magnetic spectrometer, based on a superconducting solenoidal magnet 
and many different tracking detectors, calorimeters and Cherenkov counters, which provide excellent 
particle identification. Particle identification is particularly important because the Collins effect which 
is needed for a global analysis with the SIDIS data is the effect occurring in the fragmentation of the 
light quarks. The fragmentation of charm quark or 6-quark is expected to dilute the effect, because 
helicity is only conserved for nearly massless quarks, thus removal from the data sample of all the events 
associated with c- or b- quarks is very important. 

2.5 Drell-Yan experiments 

The measurement of the Drell-Yan cross-sections is not straightforward due to the smallness of the cross- 
section. Up to now only unpolarised hadrons have been used as either beams or targets, and the quark 
PDF's have been extracted from the coefficients of various angular modulation in the cross-section. 

The NAlO experiment [671 [68] impinged a high intensity (~ 2 ■ 10^ particles per burst, 95% vr", 
corresponding to a mean intensity of 2 • 10^/sec and a duty cycle of 30%) negative beam from the 
CERN SPS unto a nuclear target, detecting muon pairs in a spectrometer whose analysing magnet had 
a hexagonal symmetry and produced a toroidal field. The experiment was run at 140 GeV, 194 GeV, 
and 286 GeV beam momenta, and different targets were used: a tungsten target (either 5.6 or 12 cm 
long) at all three beam momenta, and a liquid deuterium target, 120 cm long, placed 2 m upstream 
of the tungsten target, for the 286 GeV runs. To reduce the flux in the spectrometer, the target was 
followed by a beam dump/hadron absorber, placed between the target and the core of the dump, at 
a distance (120 cm) such that there was no contamination in the data from muon-pairs created in the 
dump. After the absorber, the trajectories of the two muons were detected with two sets of multiwire 
proportional chambers, one upstream and one downstream of the spectrometer magnet. In the analysis. 



J/\E' and T events were eliminated from the muon-pair sample by suitable cuts in the pair invariant 
mass, so that the final sample of ~ 0.3 ■ 10^ events contained only DY events in the continuum. Very 
large modulations in cos 20, where is the azimuthal angle between the hadron plane and the lepton 
plane, have been measured, in strong disagreement with the prediction of coUinear QCD. 

Similar results have been obtained at Fermilab by the experiment E615 [69]. The experiment was 
carried out in the Proton- West High intensity Area at Fermilab, and the measured reaction again was 
TT~N — )• /i~/i+X. The negative hadrons (93% 7r~) beam momentum was set at 252 GeV, and had an 
intensity of 2 x 10^ particles/sec, with a duty cycle of 33%. The experimental lay-out is similar to 
that of the CERN experiment, a nuclear target (a 20 cm long tungsten cylinder) was followed by a 
long (7.3 m) absorber, acting both as beam stopper and hadron absorber. The absorber was made of 
light material (beryllium-oxide bricks first, and graphite afterword), inserted between the pole faces of 
a dipole magnet, in order to sweep away from the spectrometer the low energy muons (corresponding 
to low-mass muon-pairs) and at the same time to focus the high mass pairs into the central part of the 
spectrometer. Low-Z material was chosen for the absorber in order to minimise the multiple scattering 
of the muons. In order to extend the measurements to very small muon angles with respect to the 
beam, the absorber was uniform in the plane transverse to the beam, and had no central plug of high-Z 
material, like the uranium/tungsten core of the CERN experiment. After the absorber, a momentum 
analysing spectrometer consisting of a system of wire chambers upstream and downstream of a second 
dipole magnet was used to measure the muon-pair trajectories. The angular distributions have been 
measured in the invariant mass region 4.05 < M < 8.55 GeV, which is free from resonances. 

Recently the E866 collaboration at Fermilab has carried out measurements of the angular distribu- 
tions of DY muon pairs produced by scattering 800 GeV protons on a deuterium and a proton target 
|70] . The Collaboration uses the upgraded Meson-East magnetic pair spectrometer at Fermilab. The 
primary proton beam, with an intensity of ~ 2 x 10^^ protons/spill impinges over one of three identical 
50.8 cm long target fiasks containing either liquid hydrogen, liquid deuterium, or vacuum. A copper 
beam dump located inside the second dipole magnet absorbs the protons that passed through the target. 
Very much as in the other experiments, downstream of the beam dump an absorber wall removes all 
hadrons produced in the target and in the beam dump. The muon trajectories are detected by four 
tracking stations (drift chambers) and a momentum analysing dipole magnet. Extrapolating the tracks 
to a vertex in the target the parameters of each muon track are optimised, the invariant mass of the 
muon pair is evaluated, and the events from the J/\l/ and T region are rejected from the DY final events 
sample. This experiment probes the DF's of the sea antiquarks, so it provides information which is 
complementary to that of the vr beam experiments. 

To summarise, the study of the DY process is very promising and several new experiments are being 
proposed. At variance with all the past experiments, the new experiments will all use polarised beams 
and/or polarised targets. 

3 Transverse-spin and transverse- momentum structure 
of hadrons 

Parton distribution functions (PDF's) and fragmentation functions (FF's) incorporate a large part of 
the information on the internal structure of hadrons that can be probed in hard processes, that is, 
in strong- interaction processes characterised by at least one large momentum scale Q. We will start 
discussing the nature and the formalism of PDF's, and then extend our presentation to the FF's. 

Right after the first DIS experiments at SLAG, Feynman proposed the concept of parton distributions 
as the probability densities of finding a parton with a certain momentum fraction inside a nucleon [71] . 
In its original formulation, Feynman's parton model was based on the observation that the time scale of 
the interaction between the virtual photon and the partons is ~ hence much smaller than the time 



scale of the binding interactions of partons, which is ~ 1/M (M is the nucleon mass) in the target rest 
frame and gets dilated in a reference system where the nucleon moves with a very large momentum (the 
"infinite-momentum frame"). Thus, we can approximately assume that in DIS the electron interacts 
elastically with free partons, and define the PDF's as the single-particle momentum distributions of the 
nucleon's constituents. 

Focusing for definiteness on quarks, fi{x) denotes the number density of quarks carrying a fraction 
X of the longitudinal momentum of the nucleon. It is not difficult to introduce polarisation into this 
simple picture. Consider first a longitudinally polarised nucleon. The helicity distribution function 
gi{x) is defined as the helicity asymmetry of quarks in a longitudinally polarised nucleon, that is, the 
number density f+{x) of quarks with momentum fraction x and polarisation parallel to that of the 
nucleon minus the number density f-{x) of quarks with the same momentum fraction but antiparallel 
polarisation: gi{x) = f+{x) — f-{x). In terms of /± the unpolarised distribution /i is simply the sum of 
the two probability densities: /i(x) = f+{x) + f-{x). The case of transverse polarisation can be treated 
in a similar way: for a transversely polarised nucleon the transversity distribution hi{x) is defined as 
the number density of quarks with momentum fraction x and polarisation parallel to that of the hadron, 
minus the number density of quarks with the same momentum fraction and antiparallel polarisation, 
that is, denoting transverse polarisations by arrows, hi{x) = ff{x) — fi{x). In a basis of transverse 
polarisation states, hi too has a probabilistic interpretation. In the helicity basis, in contrast, it has no 
simple meaning, being related to an off-diagonal quark-hadron amplitude. 

Moving from this intuitive approach to quantum field theory, the PDF's admit a rigorous definition in 
terms of correlation functions of parton fields taken at two space-time points with a light-like separation 
[72l [73] (a modern treatment is given in Ref. [H]). Since DIS probes the parton dynamics on the light- 
cone (see, e.g., Ref. [75]), it is convenient to introduce here some notions concerning light-cone geometry. 
The light-cone components of a four vector are defined as = (a° ± a^) / \/2, and grouped in triplets 
of the form = (a"^,a~,aT), where the transverse bi-vector is = (a^,a^). The norm of a'^ is 
given by = 2a+a^ — a|,. It is customary to define two light-like vectors = (1,0, O-r) and n_ = 
(0, 1,0t), sometimes called "Sudakov vectors", which identify the longitudinal direction and are such 
that n+ • n_ = 1. Any vector can be written as = a'^n'^ + a~n^ + a^, where = (0, 0, ax)- This 
is the four- dimensional generalisation of the familiar decomposition of a three-vector into longitudinal 
and transverse components with respect to a given direction. The reference frame of DIS (or SIDIS) 
is chosen so that the nucleon's momentum is purely longitudinal: ^ P'^n'^, where the approximate 
equality means that we are neglecting the nucleon mass (a legitimate approximation in the deep inelastic 
limit). The infinite momentum frame corresponds to P~^ — )• oo. Dominant contributions to DIS are 
0{P'^), whereas subleading corrections are suppressed by inverse powers of P~^, or equivalently, in terms 
of the momentum transfer, by inverse powers of Q. 

The field-theoretical expression of the quark number density fi{x) is (we postpone the QCD sub- 
tleties) 

/i(^) - I dre'^^"^-(iV|V;(V)(0)^(V)(O|iV) , (1) 

with ^+ = 0, = Ot. The peculiarity of eq. ([T]) is the appearance of the so-called "good" components 
'?/'(+) of the quark fields ip. These admit the general decomposition ip = ?/^(+) + with ip(^±) = 

i -yT j± ^_ 'pj-^g good components ■?/'(+) are the dominant ones in the infinite momentum frame, whereas 
the "bad" components ip^-) are not dynamically independent: using the equations of motion, they can 
be eliminated in favour of "good" components and terms containing quark masses and gluon fields. Due 
to the structure of eq. (JT]), one can insert between the quark fields a complete set of intermediate states 
\X), obtaining a modulus squared: 

Mx) ~ ^(P^ - - |(iV|^(+)(0)|X)p . (2) 

X 



Figure 5: The quark correlation matrix 



This expression confirms in a field-theoretical form the probabilistic interpretation of PDF's: fi{x) is the 
probability to extract from the nucleon a quark with longitudinal momentum xP^, leaving an inter- 
mediate state X of longitudinal momentum P^. A similar reasoning applies to polarised distributions 
and yields 

g,{x) ~ $^5(P+-a:P+-P+){|(iV|Q+^(+)(0)|X)r-|(iV|Q_^(+)(0)|X)r} , (3) 

X 

hix) ~ $^5(P+-a;P+-P+){|(iV|Q,^(+)(0)|X)p-|(iV|Q,^(^)(0)|X)p} , (4) 

X 

where Q± and Q-^i are the helicity and transversity projectors, respectively. 

If the quarks are perfectly collinear with the parent hadron, the three distribution functions we have 
mentioned so far, gi{x), hi{x), exhaust the information on the internal dynamics of hadrons 

at leading twist, i.e., at zeroth order in 1/Q (for an operational definition of twist, see Ref. [75]). If 
instead we admit a non negligible quark transverse momentum, the number of distribution functions 
considerably increases. At leading twist, there are eight of them. In order to understand their origin 
and meaning, it is necessary to adopt a more systematic approach. 



3.1 The quark correlation matrix 

Quark distribution functions are contained in the correlation matrix $ (Fig. E]), defined as 

$,,(A:, P, 5) = 1^ e''^-^ (P, 5|V^,(0) W[0, ^] UO\P, S) • (5) 

where |P, S) is the nucleon state of momentum P'^ and polarisation vector S^, i and j are Dirac indices 
and a summation over colour is implicit. The Wilson line W, which guarantees the gauge invariance of 
the correlator, is a path-ordered exponential of the gluon field (see below) arising from multigluon final 
state interactions between the struck quark and the target spectators. The presence of this gauge link 
introduces in principle a path-dependence in which in some cases turns out to be highly non trivial 
(Section 13.31) . 

Integrating $(A;,P, 5) over the quark momentum, with the condition x = k~^/P~^ that defines x as 
the fraction of the longitudinal momentum of the nucleon carried by the quark, yields 

$(x) = Jd'k^k,P,S)6ik+-xP+) = y^e'^'^^«"(P5|^(0)>V-[0,e]^(0|P,5)|^+=o,s,=o, (6) 

where the Wilson line >V~[0,^] connects (0,0, Ot) to (0,^~,Ot) along the n~ direction and reads (V 
denotes path ordering) 

>V-[0,e] =Pexp (-ig d^-A+(0,^-,Or)j . (7) 



In the light-cone gauge, A'^ = 0, the Wilson link reduces to unity and can be omitted. The situation is 
more complicated in the case of transverse-momentum distributions, which are defined in terms of field 
separations of the type (0,^~,$t)- we shall return to this issue in Section [3T3l 

contains the collinear (i.e., x-dependent, /cT-integrated) quark distribution functions. In- 
troducing the longitudinal and transverse components of the the polarisation vector of the nucleon, 
= (Sl/M) P^^ + S*^, the expression of $(x) at leading twist, that is at leading order in P"*", is 

$(a:) = ^|/i(x)#^ + Si^?i(x)75#+ + /ii(x)75i^^| . (8) 

Here one sees the three distributions already introduced: the number density /i(x), the helicity dis- 
tribution gi{x) and the transversity distribution hi{x)^ first identified by Ralston and Soper p5]. The 
quark distributions can be extracted from ([8]) by tracing $ with some Dirac matrix F. We will use 
the notation <I>[^l(a;) = | Tr [$(a;) V]. The explicit expressions of the leading-twist distributions are (the 
transverse Dirac matrix 77- is either 7^ or 7^): 

h{x) = <l>[^"l(x) = I ^e-^"«-(P,^|V^(0)W-[0,e]7+^(OI^,^)k+=o,S.=o., (9) 

g, (x) = $[^"^^](x) = I ^e'-^"«"(P,S|V^(0)W-[0,e]7^75V^(OI^,5)lc+=o.^.=o. , (10) 

h, {x) = $[^^^-^^](x) = I ^e^^^^«-(P,^|V^(0)>V-[0,e]7+7T75V^(OI^,^)k+=o,4.=o, • (H) 

In QCD the operators appearing in flQl fTT]) are ultraviolet divergent, so they have to be renormalised. 
This introduces a scale dependence into the distribution functions, fi{x) — t- fi{x,fi), etc., which is 
governed by the renormalisation group equations, the well known DGLAP equations [THl EH [TS] . 



3.2 The transversity distribution 

The main properties of the "third" parton density, the transversity distribution hi, eq. ffTTl) . are: 

i) it is chirally-odd and therefore does not appear in the handbag diagram of inclusive DIS, which cannot 
flip the chirality; in order to measure hi, the chirality must be flipped twice, so one always needs two 
hadrons, both in the initial state, or one in the initial state and one in the final state, and at least one 
of them must be transversely polarised ; 

ii) there is no gluon transversity distribution: this would imply a helicity-fiip gluon-nucleon amplitude, 
which does not exist since gluons have helicity ±1 and the nucleon cannot undergo an helicity change 
of two units. 

The DGLAP equations for hi have been worked out at leading order pi], and years later at next- 
to-leading order [791 IMl |8T]. There are two noteworthy features of the evolution of hi. first of all, 
since there is no gluon transversity distribution, hi does not mix with gluons and evolves as a non- 
singlet density [21]; second, at low x, hi is suppressed by the evolution with respect to gi [82j. This 
has important consequences for those observables that involve hi at low x and large Q^, such as the 
Drell-Yan double transverse asymmetry at collider energies [53] . 

The transversity distribution satisfies a bound discovered by Soffer [M] : 

\hiix)\<^[fiix)+giix)]. (12) 

This inequality, which is derived in the context of the parton model from the expressions of the distribu- 
tion functions in terms of quark-nucleon forward amplitudes, is strictly preserved in leading-order QCD 



[521 [55] . At next-to-leading order, parton densities are not univoquely defined, but a regularisation 
scheme can be chosen such that the Soffer inequahty is still valid [81] . 

The integral of over x gives the local matrix element (P, S'|?/'(0)-?/'(0)|P, 5*), which can be 

parametrised in terms of the vector, axial and tensor charge of the nucleon. In particular, the tensor 
charge (that we call 5g, for the flavour q) is given by the matrix element of the operator ipia'^'^'-f^ip, 

(P, ^|^,(0)ia^^75V^,(0)|P, S) = 25q {S^P'' - S'^P^) , (13) 

and is related to the transversity distributions as follows 

/ dx[h\{x) -h\{x)] =Sq. (14) 
Jo 

Note that, due to the charge-conjugation properties of ipia'^'^'j^ip, which is a C-odd operator, the tensor 
charge is the first moment of a flavour non-singlet combination (quarks minus antiquarks). 

An important distinction between transverse spin and transverse polarisation [22j is in order. The 
transverse spin operator, i.e. the generator of rotations, for a quark is T,t = 757o7r, and does not 
commute with the free quark Hamiltonian Hq = azPz- Thus, there are no common eigenstates of Sj- 
and Hq-. said otherwise, in a transversely polarised nucleon quarks cannot be in a definite transverse spin 
state. The distribution related to E^, called grix), is a twist-three quantity that reflects a complicated 
quark-gluon dynamics with no partonic interpretation. On the other hand, the transversity distribution 
hi carries information about the transverse polarisation of quarks inside a transversely polarised nucleon. 
The transverse polarisation operator is IIt- = ^'-/qJ^x, and commutes with Hq, owing to the presence of 
an extra 70. Therefore, in a transversely polarised nucleon, quarks may exist in a definite transverse 
polarisation state, and a simple partonic picture applies to hi. 

The argument above shows that the integral / dx (hl + hl) does not represent the quark + antiquark 
contribution to the transverse spin of the nucleon. A transverse spin sum rule containing the first 
moment of hi + hi has been derived in Ref. within the parton model, but, in the light of what we 
have just said and of other general considerations, is subject to some controversy (see the discussion in 
Ref. [M])- A sum rule for the total angular momentum of transversely polarised quarks in an unpolarised 
hadron [871 [88|, involving the generalised parton distributions, will be introduced in Section [331 



3.3 The transverse-momentum dependent (TMD) distribution functions 

The intrinsic transverse motion of quarks, is a source of azimuthal and spin asymmetries in hadronic 
processes. Taking into account its transverse component, the quark momentum is given by = 
xP'^n'^ + k!^. As we will see later in this Section, at leading twist there are eight TMD distributions: 
three of them, once integrated over k^, yield fi,gi,hi; the remaining five are new and vanish upon 
integration. Integrating ^{k,P,S) over k~^ and k~ only, one obtains the fcy-dependent correlation 
matrix 

$(x,fcr) = jdk+ jdk-<^>{k,P,S)S{k+ -xP+), (15) 

which contains the TMD distribution functions. The field-theoretical expression of ^{x, kx) l89j turns 
out to be quite complicated due to the structure of the gauge link, which now connects two space-time 
points with a transverse separation. One has [90| [91 J 



$(a;, k 



T 



27r J (27r)2 

X (P, ^I^(O) >V"[0, 00] W^[Ot, oor] W^[ooT, ^t] >V-[oo, mO\P, S) |g+=o , (16) 



(0,0) 



C 



(00,00) 
(00,0) 



(00,CXDj 



(00,0 



P 



Figure 6: The gauge-link structure of TMD distributions in SIDIS. 



with two longitudinal Wilson lines directed along n~, from (0, 0, Ot) to (0, 00, 0^) and from (0, 00, ^t) to 
(0, $,t)i and two Wilson lines at = 00 containing the transverse gluon field (Fig. IH]). This 
link structure, with the longitudinal Wilson lines W~ running to = +00, applies to semi-inclusive 
deep inelastic scattering. In Drell-Yan processes, the Wilson line runs to —00 and this may change the 
sign of the distributions, as we will discuss later (Section l3.3.ip . 

It is important to stress in eq. (fT6|) the presence of the transverse links, which survive in the 
light-cone gauge = 0, enforcing gauge invariance under residual gauge transformations. These 
transverse links are responsible for the final-state or initial-state interactions that generate some TMD 
distributions otherwise forbidden by time-reversal invariance (the so-called T-odd distributions). In 
non singular gauges, on the contrary, the gauge potential vanishes at infinity and one is left with the 
longitudinal links. It is known that in this case there are light-cone logarithmic divergences arising in the 
limit z"*" — )■ ^89j due to contributions of virtual gluons with zero plus momentum, i.e., with infinitely 
negative rapidity. One way to avoid these singularities is to use Wilson lines slightly displaced from the 
light-like direction. This introduces a dependence of the TMD distributions on a new scalar quantity, 

= (2P ■ v)^/v'^ {v is a vector slightly off the light-cone), acting as a rapidity cutoff. The light-cone 
divergences now appear as large logarithms of (, which are resummed by the so-called CoUins-Soper 
equation [HH1IS2]- A lucid presentation of this subject is contained in Ref. [71] 

Coming back to the quark correlator, at leading twist $(x, k^) has the following structure 15 uj^ 
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where is the two-dimensional antisymmetric Levi-Civita tensor, with = 1. By tracing $(x, fc^) 



-'^In the "Amsterdam classification" of TMD distributions |48J wliicli we follow in this review the letters f,g,h refer to 
unpolarised, longitudinally polarised, and transversely polarised distributions, respectively (as first proposed by Jaffe and 
Ji [221^3]). The subscript 1 labels the leading twist. Subscripts L and T indicate that the parent hadron is longitudinally 
or transversely polarised. A superscript _L signals the presence of factors in the quark correlation function. 



with Dirac matrices, $'^1 = |Tr (r$), one gets 



= /i(x,4)-^2%^/i-'t(^,4), (18) 
$[7+75] ^ SLgiL{x,kl) + ^^:^^g,T{x,kl), (19) 

- '^'^ % /.^.(x, 4) - %^^(^^ 4) , ^ = 1,2. (20) 

where /ii = hiT+{k'^/2M)h^j.. The three quantities $[7+]^ $[7"'' 75] g^^^j <|)[io''+75] represent the probabihties 
of finding an unpolarised, a longitudinally polarised and a transversely polarised quark, respectively, 
with momentum fraction x and transverse momentum k^- In eqs. (11811201) eight independent TMD 
distributions are present: /i, /fi^, gi, gir, hi, h^i^, /i^j., h^. Upon integration over fcy, only three of these, 
/i(x, k^),gi{x, k^), hi{x, k^), survive, yielding the x-dependent leading- twist distributions fi{x), gi{x), 
hi{x). 

From eq. fl20l) one sees that the spin asymmetry of transversely polarised quarks inside a transversely 
polarised nucleon is given not only by the unintegrated transversity hi{x,k^), but also by the TMD 
distribution hij,{x,k^), which has been given the name of "pretzelosity" , as it is somehow related to 
the non-sphericity of the nucleon shape [93] (for a review of the properties of hij^, see Ref. [Sj). Note 
that, due to the intrinsic transverse motion, quarks can also be transversely polarised in a longitudinally 
polarised nucleon (/i^^.); longitudinally polarised in a transversely polarised nucleon (gir)- 

3.3.1 The T-odd couple: Sivers and Boer-Mulders distributions 

From eq. f|T8|) the probability of finding an unpolarised quark with longitudinal momentum fraction x 
and transverse momentum kr inside a transversely polarised nucleon is 

f,/^,ix, kr) = /i(x, 4) - X ^) ■ 4) , (21) 

where P = P/\P\. Thus the azimuthal asymmetry is 

f,/^,{x, kr) - f,/m{x, -kr) = -2 (-^ >< ^) " 4) ^ (22) 

which is proportional to the so-called Sivers function [511 ES]- ^ '^on vanishing signals that 
unpolarised quarks in a transversely polarised nucleon have a preferential motion direction: in particular, 
> means that in a nucleon moving along +z with transverse polarisation in the +y direction, 
unpolarised quarks tend to move to the right, i.e. towards —x. 

Specularly, the distribution of transversely polarised quarks inside an unpolarised nucleon is 



fqt/Nix^kr) 



2 

and from this we get a spin asymmetry of the form 



f ( (-P X fcy) ■ SgT j^±f ,2\ 

ji{x, Kj:) — n,-^ [X, itnp) 



(23) 



(P X kx) ■ S T 

fqr/N{x, kx) - fqi/N{x, kx) = — — h^ix, kx) , (24) 



which is proportional to h^, the Boer-Mulders distribution ^U\. Positivity bounds for and hj- 
were derived in Ref. [96]. Note that in the literature (see Ref. |63] and bibliography therein) one also 
encounters the notation 

The Sivers and Boer-Mulders functions are associated with the time-reversal (T) odd correlations (P x 
kx) ■ St and (P x kx) ■ SgT, hence the name of "T-odd distributions". To see the implications of 
time-reversal invariance one has to recall the operator definition of these distributions which, in the 
case of the Sivers function, is: 

/f!r(x,4) ~ J dr|d2|^e^-^^«"-^-^-(P,5T|?(0)7+W[0,e]V^(OI^,^T)k+=o (26) 

If the overall Wilson link W is naively set to unity, the matrix element in f l2B]) changes sign under time 
reversal, and the Sivers function must therefore be zero [29]. On the other hand, a direct calculation [36] 
in a spectator model shows that is non vanishing: gluon exchange between the struck quark and the 
target remnant generates a non-zero Sivers asymmetry (the presence of a quark transverse momentum 
smaller than Q ensures that this asymmetry is proportional to M/ kx-, rather than to M/ Q, and therefore 
is a leading- twist observable). The puzzle is solved by carefully considering the Wilson line in eq. (1261) 
[38j . In fact H^[0,^] includes transverse links at infinity that do not reduce to unity in the light-cone 
gauge [91]. Since time reversal changes a future-pointing Wilson line into a past-pointing Wilson line, 
T-invariance, rather than constraining to zero, gives a relation between processes that probe Wilson 
lines pointing in opposite time directions. In particular, since in SIDIS the Sivers asymmetry arises 
from the interaction between the spectator and the outgoing quark, whereas in Drell-Yan production it 
arises from the interaction between the spectator and an incoming quark, one gets 

/{^(x, /cr)siDis = -fixix, A^Ddy • (27) 

A similar relation holds for the Boer- Mulders function hj;. Eq, ( 1271) is an example of the "time- reversal 
modified universality" of distribution functions in SIDIS, DY production and e~^e~ annihilation studied 
in Ref. [97]. The relation (!27|) is a direct consequence of the gauge structure of parton distribution 
functions, and its experimental check would be extremely important. 

Gauge link patterns of hadroproduction processes are more complicated and do not result in a simple 
sign flip of the distributions [981 M\ HDQl [Ml HDS] . For these processes the authors of Refs. [991 HTin] 
suggested that a factorisation scheme should hold with kx distributions containing pro cess- dependent 
Wilson lines. This "generalised TMD factorisation" evidently differs from the standard TMD factori- 
sation, wherein the k^ distributions are fully universal quantities. A recent study [104] . however, has 
shown that even the generalised factorisation scheme is violated in hadroproduction of nearly back- 
to-back jets or hadrons, a process investigated experimentally by the STAR collaboration at RHIC 

m- 

The quark Sivers function has an exact gluonic counterpart, f^, which represents the distribution 
of unpolarised gluons in a transversely polarised hadron. This function is called Gt in Ref. [65], where 
a complete classification of leading-twist fcr-dependent gluon distributions is presented. There is no 
gluonic equivalent of the Boer-Mulders function, but a somehow similar quantity is the distribution of 
linearly polarised gluons in an unpolarised hadron. 

A sum rule for the Sivers function was derived in QCD by Burkardt |106[ 1107] , who showed that the 
sum of all contributions to the average transverse momentum of unpolarised partons in a transversely 
polarised target (that is, the average transverse momentum induced by the Sivers effect), must vanish: 



J2 (feT)|sivers = 

a=q,q,g 



(28) 
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Figure 7: The quark-gluon correlation matrix. 
In terms of the Sivers function, the condition fj28|) becomes |108] 

5^ j\xfti'^\x) = Q, (29) 

where the first A;|.-moment of is given by 

ft^'\x)^ |d^fc^^/^4,(x,4). (30) 

Although some QCD aspects, such as ultraviolet divergences and light-cone singularities, were not 
considered in the original derivation, eq. (!30|) is likely to be valid in general. In Ref. |109] it has been 
shown that the Burkardt sum rule is fulfilled for a quark target in perturbative QCD at one-loop order. 

From a phenomenological point of view, the importance of eq. f l29|) is that one can infer the size 
of the gluon Sivers function from fits to SIDIS observables involving the quark and antiquark Sivers 
functions fTTOlfTTT] . 



3.4 Higher-twist distributions and quark-gluon correlators 

At twist three, suppressed by l/Q, i.e., by 1/P+ in the infinite momentum frame, with respect to 
leading twist, the quark correlator $(x) admits the general decomposition 



(a:;) I twist 3 = |e(x) + gT{xh5$T + Sl hL{x) M±l|_j2^| ^ (31) 

displaying the three distribution functions e{x), grix), hL^x). In particular, grix) contributes to the 
polarised DIS structure function g2{x,Q'^) (see, e.g., Ref. [112] ). Higher-twist distributions do not have 
a probabilistic interpretation. They involve in fact both good and bad components of the quark fields, 
so the procedure leading to expressions such as eqs. (^j^ cannot be applied. 

Higher-twist effects are a manifestation of quark-gluon correlations inside hadrons |113[ 1114] . At 
twist three there are four quark-gluon correlators, which depend on two momentum fractions, x and 
x' (see Fig. ED: Gpix^x'), Gir(a;,x'), Hf{x,x'), Ef{x,x'). In the literature [53l |55l [561 |60l [HS] , these 
correlators are also called Tp, Tp, respectively, with a normalisation varying from one paper 

to another. In QCD the quark-gluon correlation functions acquire a dependence on a scale /i. The 
equations governing the evolution in /i have been recently written down and solved EO] • 

Using the QCD equations of motion and integrating over x' one can show that 1116] 



g^^:,) = m:^ + g^^a:), h^ix) = -2 '-^^^^ + h{x) , e{x) = e{x) , (32) 



where g\rf, h^j^ are the first transverse moments of qit and h-^-^, defined as in eq. ( !30l) . and the tilde 
functions are genuinely twist-three distributions related to the quark-gluon correlators. 

Ignoring the contributions of tilde functions and of quark mass terms one gets the generalised 
Wandzura-Wilczek (WW) approximation, so called in analogy with the original WW relation |117] 
between the polarised DIS structure functions gi{x,Q'^) and g2{x,Q'^) \IIS\ 11121 1119] . The generalised 
WW approximation relates twist-three distributions to twist-two distributions. It has been investigated 
by various authors [120^ 11211 1122j and also applied in phenomenological analyses }123]. 

Coming to the transverse motion of quarks, the structure of the fcT-dependent quark correlator 
$(a;, kr) at twist three has also been studied by various authors [HI EHl I124[ 1125] . It is now known 
that there are 16 twist-three TMD distributions: e, ej^, e^, ct, /t, /l , /t , /"*", gr, gi, gr-, hi, hx, 
h, /ly in the classification of Ref. |126] ). Some of these functions, namely g^,ej^, fr, /t: not identified 
in earlier studies, exist because the Wilson line in the quark correlator provides an extra independent 
vector (n_) for the Lorentz decomposition of ^{x, fcy). If we integrate $(x, kj^) over k^, the only non 
vanishing distributions are the three T-even functions in eq. (13T|) . 

Concerning twist four, the integrated parton distributions were first identified in Refs. [221 ESI ■ 
More recently, the complete expression of the fcr-dependent correlator $(a;, fcy) has been given in Ref. 
[125] , where it is shown that up to twist four there are in total 32 TMD distributions. The unintegrated 
correlation matrix $ is also composed of 32 Lorentz-scalar structures: 12 amplitudes associated to the 
four- vectors k,P,S and 20 amplitudes associated to n_. The number of distribution functions being 
equal to the number of amplitudes of $, all the TMD distributions are independent and there are no 
general relations among them. In earlier studies |120l 1^. some "Lorentz- invariance relations" (LIR's) 
were derived from an expansion of $ that did not take into account the amplitudes associated to the 
gauge link vector n_. Two of these relations are: 

qAx) = gi{x) + ^gf-^^\x) , hL{x) = hi{x) - -^hf^ (x) . (33) 

The presence of the n_-dependent amplitudes invalidate the LIR's, which are not valid in QCD |128[ 
1129] . However, they approximately hold in the generalised WW approximation [122] . 

A general remark about higher- twist distributions is in order. While the distributions e{x), grix), 
hiix) - or, to be precise, the corresponding quark-gluon correlators - enter into the collinear twist- 
three factorisation theorem of QCD ^53], i54j , the /cr-dependent higher-twist distributions are employed 
in factorisation formulae that lack a solid QCD foundation. Thus, they should rather be interpreted as 
a way to model subleading effects. 

3.5 Generalised parton distributions 

The generalised parton distributions (GPD's), which are related to non-forward quark-quark (or gluon- 
gluon) correlators, emerge in the description of hard exclusive processes, such as deeply- virtual Compton 
scattering and exclusive meson production, characterised by a non-zero momentum transfer to the target 
nucleon |130l I13H I132[ [66| I133[ 1134] . Here we will be mostly concerned with the relations existing 
between the GPD's and the transverse spin distributions (for more details, see Ref. [135] ). 

The kinematics of GPD's is represented in fig. [H] (we follow the conventions of |SS])- The momenta of 
the incoming and the outgoing nucleon are p = P — |A and p' = P + |A, respectively. The momentum 
transfer squared ist = A^. The GPD's depend on t and on two light-cone momentum ratios: x = k~^/ P^ 
and ^ = — A+/2P+. The variable ^ is sometimes called "skewness", and the GPD's are also known as 
"skewed parton distributions". At leading twist, there are 8 GPD's |136] : 



H{x,^,t), E{x,^,t), i^(x,e,t), ^(x,e,t), HT{x,^,t), ET{x,^,t), H{x,^,t), E{x,i,t) 



(34) 




Figure 8: Kinematics of GPD's 



The first four are chirally even and are related to the famihar form factors. Integrating H , E, H, E 
over X, one gets in fact the Dirac, Pauh, axial and pseudoscalar form factors, respectively. The quantity 



dxE«(x,0,0) = k'^ 



(35) 



is the contribution of the flavour q to the anomalous magnetic moment of the nucleon, that is, to the 
Pauli form factor F2 ai t = 0. The GPD's H, H, Ht, taken at ^ = t = 0, coincide with the integrated 
quark distributions fi,gi,hi: 



H{x, 0, 0) = /i(x) , H{x, 0, 0) = gi{x) , Ht{x, 0, 0) = h{x) 



(36) 



The original interest in GPD's was prompted by Ji's sum rule relating the total angular momentum of 
quarks (in a nucleon with polarisation vector S) to the second moment of H and E |131] : 



( Jj) =S' jdxx [H{x, 0, 0) + E{x, 0, 0)] . 



(37) 



A similar decomposition for the angular momentum of quarks with transverse polarisation vector Sq in 
an unpolarised nucleon has been derived in Ref. [871 EH] and is given by: 



(j;(5,)) = ^ I dxx [Hrix, 0, 0) + 2HTix, 0, 0) + Erix, 0, 0)] . 



(3J 



Here Ht{x, 0, 0) coincides with transversity, whereas the combination 211^ + E^ appears in the impact- 
parameter description of the Boer-Mulders effect (see Section 13. 6p . 

Note in conclusion that there are no direct and model-independent connections between the GPD's 
and the TMD distributions, as stressed in Refs. |135[ 1137] . GPD's are instead directly related to the 
distribution functions in the impact-parameter space. 



3.6 Distribution functions in the impact-parameter space 

In the impact-parameter space one can get a more intuitive picture of some transverse spin and trans- 
verse momentum effects. To define the impact-parameter distributions (IPD's), one first introduces 
nucleon states localised at a transverse position Rt, by means of an inverse Peierls-Yoccoz projection: 

|P+, Rt, S)=Ar I ^ e-'^-^- |P, S) , (39) 

where A/" is a normalisation factor. The IPD's are light-cone correlations in these transverse-position 
nucleon eigenstates. For instance, the unpolarised IPD is given by 

g(x,6|) = I ^e^"^^^- (P+,Ot;5|V5(z0>V-[zi,Z2]7+^(^2)|P+,Ot;5), (40) 



with Zi^2 = (0'*', =f|-2^~5 &t)- This is the number density of quarks with momentum fraction x and 
transverse position inside an unpolarised hadron. The polarised IPD's are obtained by inserting in 
the matrix element of eq. fHOj) . instead of 7"'", the matrices 7''' 75 and i(j*^75. 

IPD's are Fourier transforms not of the TMD distributions, but of the GPD's. The impact-parameter 
transform of a generic GPD X for ^ = (which implies = — A^) is defined as 

^^""^ = / °' • ^^^^ 

It is straightforward to show [138] that the unpolarised IPD q{x, 6^) coincides with the Fourier transform 
of H{x,0, -At), that is 

g(x,6|) =-H,(x,6|). (42) 

The impact-parameter density of unpolarised quarks in a transversely polarised nucleon (A^^) is |138[ 
[T391IT35] 

q^,{x, br) = n,{x, 6|) + (-^ >< ^) " ^^ S'^{x, b'^) , with S'^{x, h^) ^ ^S,{x, hi) , (43) 

where £{x, 6^) is the Fourier transform of E[x, 0, — A|,). Notice the formal similarity with eq. f l2T|) and 
the correspondence /fi^ —£' . . Due to the E'^ term, which can be regarded as the fo^-space analogue 
of the Sivers distribution, q^t{x,bT) is not axially symmetric and describes a spatial distortion of the 
quark distribution in the transverse plane. Final-state interactions can translate this position-space 
asymmetry into a momentum-space asymmetry. For instance, if the nucleon moves in the +z directions 
and is polarised in the +x direction, a positive £g implies that quarks tend to be displaced in the 
—y direction, and final-state interactions, which is expected to be attractive on average, convert this 
transverse distortion into a momentum asymmetry in the +y direction. This is the intuitive explanation 
of the Sivers effect in the impact-parameter picture |1U61 ll4Uj . A measure of the space distortion is 
given by the flavour dipole moment 

di = Jdx Jd'bTbirqr,rix,bT) = -'^ J dx E,ix,0,0) = -'^^ , (44) 

where k'^ is the contribution of the quark flavour q to the anomalous magnetic moment of the nucleon, 
see eq. (1351) . The argument developed so far is summarised by the following qualitative relation between 
the Sivers function and k"^ |14H 11421 1140] (any quantitative relation between these two quantities 
is necessarily model-dependent |135] 1137] ) 

/it' ~ , (45) 

where the minus sign is a consequence of attractive final-state interactions that transform a preferential 
direction in the 6T-space into the opposite direction in kx- Eq. fH5|) leads to an immediate prediction: 
since the quark contributions to the anomalous magnetic moment of the proton k,p, extracted from the 
experimental value of using SU(2) flavour symmetry, are ~ 1.7, k'^ = —2.0, one expects < 
and > 0. This prediction has been corroborated by the SIDIS experiments (see Section l5.3.ip . 

Consider now the case of transversely polarised quarks inside an unpolarised nucleon. Their impact- 
parameter distribution is |139] 

qH^^br) = \ \H,{x,hl) + ' ^'"^ [S^x, hi) + 2g- ,(x, h^)] | , (46) 
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Figure 9: First x-moments of the densities of unpolarised quarks in a transversely polarised nucleon 
(left) and transversely polarised quarks in an unpolarised nucleon (right) for u (upper plots) and d 
(lower plots) quarks. Quark spins (inner arrows) and nucleon spins (outer arrows) are oriented in the 
transverse plane as indicated. From Ref. |143j . 



The term S' + 2H^ is the analogue of the Boer-Mulders function in the far-space - see eq. (l23l) . Again, 
we see that transverse spin (of quarks, in this case) causes a spatial distortion of the distribution, which 
is at the origin of the Boer-Mulders effect. One can repeat the same reasoning developed for the Sivers 
effect and introduce a transverse anomalous moment k^, defined by 

4 = jdx[E^{x,0,0) + 2H^{x,0,0)]. (47) 

The Boer-Mulders function is expected to scale with this quantity, 

h^^ ~ -4 . (48) 

where the minus sign has the same meaning as before. Unfortunately, no data exist for k^. This quantity, 
however, and the impact-parameter distributions have been calculated in lattice QCD |143i ,144 J and 
are shown in Fig. [9]). The result for is: = 3.0, 4 = 1-9- Thus, at variance with f^, we expect 
the u and d components of to have the same sign, and in particular to be both negative. Moreover, 
assuming simple proportionality between the ratio hj^ / frr k/kt, the u component of should be 
approximately twice as large as the corresponding component of /fi^, while hj^'^ and should have a 
comparable magnitude and opposite sign. These predictions are well supported by a phenomenological 
analysis of SIDIS data [145] as will be shown in Section 15.3.11 

3.7 Model calculations of TMD distributions 

Models and other non-perturbative approaches like lattice calculations play a very important role when 
the experimental information about distribution functions is scarce or lacking at all. So, it is not 
surprising that a considerable effort has been made to compute the TMD distributions in various 
models of the nucleon and by lattice QCD. Here we will be not be able to give an exhaustive account 
of all this work still largely in progress and we will limit ourselves to a general discussion. For a recent 
review of model results see Ref. |146j . 



The first calculation of TMD distributions was performed in a quark-diquark spectator model [147] . 
This class of models, with various quark-diquark vertex functions, has been subsequently used by 
many authors. In particular, using a simple scalar spectator model with gluon exchange is was shown 
explicitly [36] that the Sivers function is non vanishing. Since Wilson links, representing gluon insertions, 
are crucial in order to guarantee the existence of the T-odd distribution functions, these can only be 
computed in models containing gluonic degrees of freedom. Following Ref. [36], more refined calculations 
of the Sivers and Boer-Mulders functions were performed in spectator models with both scalar and axial- 
vector diquarks and various quark-diquark vertices |148[ I149[ I150[ I109[ 11511 1152[ 1153] . Other models 
used to evaluate the T-odd functions include the MIT bag model [1541 11551 11561 1157] . the constituent 
quark model |158[ 1157] and a light-cone model |159] . In Ref. |155] final state interactions were assumed 
to be induced by instanton effects. 

What emerges from models is that the Sivers function, although quite variable in magnitude, is 
negative for u quarks and positive for d quark. A different sign of is however found in the model of 
Ref. |155] . As for the Boer-Mulders function, the general prediction (with the exception of Ref. |150] ) 
is that both the u and the d distributions are negative. These signs for and hj^ are also expected in 
the impact-parameter picture [138^ 1141^ 1140^ 1142^ IS7] , in the large- A^^^ approach [16U] , which predicts the 
isoscalar component of and the isovector component of to be suppressed, and in chiral models 

Spectator models have been also used |162[ 1163] to calculate T-odd twist-3 distributions, in particular 
g-^, which contributes to the longitudinal beam spin asymmetry in SIDIS. 

Models without gluonic degrees of freedom can be used to compute T-even TMD distributions 
only. These distributions have been calculated in a spectator model ^152] . in light-cone quark models 
[164[ 11651 1166] , in a covariant parton model with orbital motion [167] and in the bag model [9H 1168] . 
In particular, Ref. |168] presents a systematic study of leading and subleading twist TMD distributions 
and of the relations among them. 

In any quark model without gluons, the Lorentz-invariance relations, obtained by neglecting the 
amplitudes of the quark-quark correlator related to the gauge link (Section 13. 4p . must obviously be 
valid. There are also a number of other relations that hold in some specific models like [91] 

gi{x, kl) - hi{x, kl) = h^ri^, k^) . (49) 

According to this relation, h^j^ can be interpreted as a measure of the relativistic effects in the nu- 
cleon, which are known to be responsible for the difference between the helicity and the transversity 
distributions |169] . Other model-dependent relations involving the TMD distributions are listed in 
Refs. [TMITBg] . 

Finally, one should keep in mind that models provide a dynamical picture of the nucleon at some 
fixed, very low, scale jj? < \ GeV^ |170[ I17H I172[ 1173] . The quark distributions that one gets are 
therefore valid at this unrealistic scale and must be evolved to the experimental scales. The evolution 
of the TMD distributions has been unknown until very recently and is therefore usually neglected or 
approximated in current phenomenological analyses. 

3.8 Fragmentation functions 

In partially inclusive processes a parton hadronises into a particle h carrying away a fraction of the 
parton's momentum. In the following it is supposed that the fragmentation process is initiated by a 
quark, as is the case in SIDIS and e^e~ annihilation at leading order. The momentum of the fragmenting 
quark is indicated as and z is the fraction of its longitudinal component carried by the final hadron, 
z = / K~ . Since the hadron moves in the opposite direction with respect to the target nucleon, 
the minus components of momenta are the dominant ones. The fragmenting quark has a transverse 



momentum Kt with respect to the final hadron. Conversely, the hadron has a transverse momentum 
Pt = —zkt with respect to the quark. 

The fragmentation analogue of fc^) is: 



^.,[z,zKr) - 2^ 2n J (27r)3 " 



X (0| W[+oo, ^i(0 1^/., X)(P,, 5,; X|^,-(0) >V[0, +oo]|0)|5-=o , (50) 

where each Wilson line includes a longitudinal link along n+ and a transverse link at infinity |126] . In 
the case of fragmentation one has the same gauge structure in SIDIS and in e~^e~ annihilation, which 
means that there is no difference between the fragmentation functions of these processes, i.e. they are 
universal quantities in a full sense |97j. The integrated fragmentation correlator is given by 



E{z) = I d'prE{z,pT) = \ \D,{z)i_ + ^^^1(^)75^, + Hr{z) ^^^'^-J^^ j + h.t. , (51) 

where "h.t." denotes higher-twist terms. Di, Gi and Hi are the integrated leading-twist fragmentation 
functions (FF's): Di is the ordinary unpolarised fragmentation function, Gi is the analogue of the 
helicity distribution. Hi is the analogue of the transversity distribution (it describes the fragmentation 
of a transversely polarised quark into a transversely polarised hadron). 

To compute azimuthal asymmetries the transverse-momentum dependent FF's are needed. For 
simplicity, we limit ourselves to listing the FF's of main phenomenological interest. The traces of the 
fragmentation matrix corresponding to unpolarised and transversely polarised quarks are 



-[--](., p^) = Diiz^pl) + ^^^4^ Di^i^.vl) , (52) 

zMh 

S[-'-^^](z,p^) = Sl^Hi{z,pl) + '^Hi{z,pl) + . . . . (53) 

D^rp is analogous to the Sivers distribution function and describes the production of transversely po- 
larised hadrons from unpolarised quarks. For this reason it is called "polarising fragmentation function" 

3.8.1 The Collins function 

The most noteworthy FF appearing in ( 153|) is H^Zjk^), the so-called Collins function, describing 
the fragmentation of a transversely polarised quark into an unpolarised hadron [29]. The resulting 
transverse-momentum asymmetry of hadrons is 

(k X Pr) ■ S'rj, 

Dh/A^,Pt) - D,/^r{z, -pt) = 2 ^ ^-L-^H^{z^pl) , (54) 

where S'^ is the spin vector of the fragmenting quark. From the structure of the correlation [k^pr)- S'^rp 
one sees that a positive H^ corresponds to a preference of the hadron to be emitted on the left side 
of the jet if the quark spin points upwards. Through this mechanism the transverse momentum of the 
produced hadron with respect to the jet direction acts as a quark polarimeter. 

The Collins function satisfies a sum rule arising from the conservation of the intrinsic transverse 
momentum during quark fragmentation. This sum rule, discovered by Schafer and Teryaev jl75] . reads 

l'dzzHt^'^''{z) = 0, with Ht^'\z)^z' fd'^T^^H^iz^z'^). (55) 



h 




Figure 10: The fragmentation process in the recursive string model |177] . The z direction is along the 
string and the fragmenting quark go is supposed to be polarised in the direction, out of the page. 



A simple qualitative explanation of the Collins effect is provided by the "recursive string model" 
of Ref. [281 1176[ 1177] . which is illustrated in Fig. [TOl Suppose that a quark go? polarised in the +^ 
direction, i.e. out of the page in the figure, fragments into a pion with an antiquark gi created by string 
breaking. If the gigi pair is in a ^Pq state, the orbital angular momentum of the pair is L = 1, and 
the pion, inheriting the transverse momentum of gi, moves in the +x direction. The quark gi, with the 
subleading pion that contains it, moves in the opposite direction. This model predicts opposite Collins 
asymmetries for 7r+ and tt^ assuming u dominance, and a positive (negative) sign for the favoured 
(unfavoured) Collins function. "Favoured" refers to the fragmentation of a quark or an antiquark 
belonging to the valence component of the final hadron, e.g. u — >■ vr"*", c? — )■ vr", c? — )• tt"*", etc.. 

The CoUins function for pions has been computed in various fragmentation models |178[ I179[ I180[ 
1181^ 1182] . What is common to these approaches is that arises from the interference between a tree 
level amplitude and loop corrections that provide the necessary imaginary parts. The differences reside 
in the pion-quark couplings and in the nature of the virtual particles in the loops (pions or gluons). An 
assessment of model calculations of the Collins function is contained in Ref. |183] . 

4 Processes and observables related to transverse spin 

In this section we will present a general description of the processes probing the transverse-spin and 
transverse-momentum structure of hadrons, and of the relevant observables: single-spin asymmetries, 
double-spin asymmetries and unpolarised azimuthal asymmetries. The focus will be on two classes of 
reactions that have clear and well-established theoretical descriptions, namely SIDIS with the related 
process e^e~ annihilation into hadron pairs, and DY production. Hadroproduction will also be de- 
scribed, but in lesser detail. The last subsection contains a sketchy presentation of other processes 
somewhat related to transverse spin. 

4.1 Semi-inclusive deep inelastic scattering (SIDIS) 

SIDIS is the process £(/) + N{P) + h{Ph) + X{Px), where £ {£') is the incoming (outgoing) 

lepton, N the nucleon target, h the detected hadron. The corresponding four-momenta are given in 
parentheses. In the following we will denote by S\\ and S± the longitudinal and transverse component 
of the target spin vector, respectively, and by Xi the longitudinal polarisation of the incident lepton. 

SIDIS is usually described in terms of the invariant variables xb = Q^/2P ■ g, y = P ■ q/ P ■ I, 
Zh = P ■ Ph/P ■ with q = I — I' and = — g^. In the deep inelastic limit, is much larger than the 



Figure 11: Lepton and hadron planes in semi- inclusive deep inelastic scattering. The reference frame is 
defined according to the convention of Ref . |184] . 



mass M of the nucleon and the mass Mh of the final hadron. Hereafter mass corrections are neglected 
unless otherwise stated. 

To parametrise the cross section in terms of structure functions the so-called "7*iV coUinear frame" 
|184] is usually adopted. As shown in Fig. [TT] in this reference frame the virtual photon and the target 
nucleon are coUinear and directed along the z axis, with the photon moving in the positive z direction, 
and the final hadron has a transverse momentum All azimuthal angles are referred to the lepton 

scattering plane: 0^ is the azimuthal angle of the hadron /i, 05 is the azimuthal angle of the nucleon 
spin vector S^. The phase space of the process contains another angle, ip, which is the azimuthal angle 
of the outgoing lepton around the beam axis with respect to an arbitrary fixed direction, which is chosen 
to be given by the target spin. Up to corrections of order M^/Q^ one has (hp ~ d05 [185] . 

In this Section we consider the case of a spinless or unpolarised detected hadron, while leptopro- 
duction of transversely polarised hadrons will be treated in Section I4.1.4[ The SIDIS differential cross 
section in the six variables XB,y, Zh, 4>s, Ph± = \Ph±\ and (ph, is given by 



dV y 



dxBdydzhd(psd(phdP^_^ 8Q^ Zh 



L^.W", (56) 



where L^j, is the usual DIS leptonic tensor and W^" is the hadronic tensor 



X 

X (P, S\J^\0)\X; P,, Sh){X; P,, (0)|P, S) . (57) 



Neglecting for simplicity the M^x^/Q^ corrections the complete SlDlS cross section can be parametrised 



in terms of 18 structure functions as follows |185[ 1126] : 
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+ (2 - v^r^ sin 05 F^^'t'^ + {2-y) sin(20. 
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^cos(2(^h-(/)s) 
LT 



(5^ 



The structure functions F depend on XB,y,Zh and P^^. Their first and second subscript denote the 
polarisation of the beam and of the target, respectively [U = unpolarised, L = longitudinally polarised, 
T = transversely polarised), whereas the third subscript refer to the polarisation of the virtual photon. 

and 



If we integrate (158!) over Ph±, only 5 structure functions survive: Fuu^t-, Fuu,l, F^l, F\ 



cos <ps 
LT ; 



FjjjJ^ . The first two, upon a further integration in z and a sum over all outgoing hadrons, yield 
the unpolarised DIS structure functions Ft{xb,Q^) = '2xFi{xb,Q^) and Fl{xb,Q'^) = F2{xb,Q^) — 
2xbFi{xb, Q^)- The second two lead to combinations of the structure functions gi{xB, Q^) and g2{xB, Q^) 
of longitudinally polarised DIS. The fifth one vanishes. The fact that 



UT 







(59) 



is a consequence of time-reversal invariance |185] and is another way to express the Christ-Lee theorem 
[26], according to which there cannot be transverse spin asymmetries in inclusive DIS. On the contrary, 
in SIDIS no first principle forbids the existence of transverse spin asymmetries. 

In the literature, the spin and azimuthal asymmetries of SIDIS are defined in two different ways: 



- as the structure function ratios: 



^XY l"'-; ^) ^h±) 



XY 



uu 



(60) 



where X and Y label the polarisation of the beam and of the target respectively, w{(f)h,4's) is a 
trigonometric function of its arguments, and Fuu = {1 — y + \y'^)Fuu,T + (1 — y)Fuu,L] 



- as the azimuthal moments of cross sections [184] : 

j d0,, j d05 da(0,„ 05) 
where da is a shorthand notation for the fully differential cross section. 
Notice that the two definitions of asymmetries differ for y-dependent factors. 



(61) 
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Figure 12: Diagram contributing to semi- inclusive DIS in the parton model. 



4.1.1 SIDIS in the parton model 

In the parton model the virtual photon strikes a quark (or an antiquark), which successively fragments 
into a hadron h. The process is represented by the diagram in Fig.[T21 We will take transverse momenta 
of quarks into account and refer to this description as the "extended parton model" . 

For the partonic description of SIDIS we work in a reference frame where the momenta of the target 
nucleon and of the outgoing hadron are coUinear and define the longitudinal direction. In this "/lA^ 
collinear frame", one has P^ = P^n^ and = P~n't, whereas the virtual photon momentum acquires 
a transverse component Qt- The incoming quark momentum is k'^ = xP^ + fc^, with x = k^/P~^; 
the fragmenting quark momentum is k,'^ = Pj^/z + k^, with z = Pf^ /k,^. Notice that the "transverse" 
quantities in the hN frame (labelled by the subscript T) differ from the "perpendicular" quantities in 
the 7*A^ frame (labelled by the subscript ±) by terms suppressed at least as 1/Q. In particular, Qt is 
related to Phi_ by Qt = —Ph±/zh, up to corrections. 

The hadronic tensor corresponding to Fig. [12] is given by 



where the second expression has been obtained by working out the momentum-conservation delta func- 
tion and neglecting 1/Q^ terms. In this case one has z = Zh and x = xb- 

Inserting the expressions of $ and S into fl62|) and contracting W^"^ with L^^, leads to the SIDIS 
structure functions. With the following notation for the transverse momenta convolutions 
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2zh J2 / d'fcT / d^i^T S\kT + qT- i^t) Tr [<^%xb, kr) ^ i^t)Y] , (62) 
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C[wfD] = Y,^l^ / d'fcr / d^i^T5\kT - i^T - PhJ z)w{kT^i^T) rix.kDD'^iz.Kl) , (63) 



the non vanishing structure functions at leading twist are [IHl EHl 1126] 
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(64) 
(65) 

(66) 
(67) 
(68) 

(69) 

(70) 

(71) 



where h = Ph±/\Ph±\- The structure function Fjju^t gives the dominant contribution to the unpolarised 
cross section integrated over (ph- 

The Colhns term p^^'^^^'^^^ couples the Collins function to the transversity distribution hi, 
thus representing one of the privileged ways to access this quantity. Note that in the original paper 
[29] the Collins angle $(7 was defined as the angle between the momentum of the produced hadron and 
the spin of the fragmenting quark, i.e. $c = (p^ — 4>s'g- In terms of the azimuthal angle of the target 
spin (ps, one has = <Ph + <Ps ~ t^- On the other hand, according to the conventions of Ref. [184] . 
the Collins angle is defined as <I>J^ = (ph + (ps- Thus, one gets different signs for the Collins asymmetry, 
depending on which definition of the Collins angle, either $(7 or is adopted. 

Another leading-twist asymmetry source is the Sivers term F^^^p"~'^^\ which couples the Sivers 
function to the unintegrated unpolarised fragmentation function Di. In the transversely polarised 
case, a further angular modulation, of the type sm{3(ph — (ps), involves the distribution function h-^j,. 

In unpolarised SIDIS, a leading-twist azimuthal asymmetry is generated by the structure function 
pcos^2<f)h ^ which couples the Boer- Mulders distribution to the Collins fragmentation function H^. 



Going to twist three, i.e. to order l/Q, it turns out that the leading-twist structure functions (16411701) 
do not acquire any extra contribution, but there appear other non vanishing structure functions [126j . 
Among them, of particular phenomenological importance are those related to the cos^^ and sin^^ 
modulations. Ignoring, in the spirit of the parton model, interaction-dependent terms, i.e. quark-gluon 
correlations, and quark mass contributions (the generalised Wandzura-Wilczek approximation) one 
finds [126] 
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(72) 



(73) 



In the same approximation one gets F^^^'^ = 0. Thus a deviation of the beam-spin siiKp^ asymmetry 
from zero might signal the relevance of interaction effects in the nucleon. One should recall however 



that at high transverse momenta F^^'^^ is non zero in next-to-leading order QCD. The term in F^^'^'^ 
containing the product of the unpolarised functions /i-Di is a purely kinematical contribution arising 
from the intrinsic transverse motion of quarks, with no relation to spin. This contribution was discovered 
longtime ago by Cahn |186tll87] . and the corresponding azimuthal asymmetry is referred to as the "Cahn 
effect". A similar contribution emerges at twist four, that is at order in the cos20/i term: 



/Vf2 

pcos 2(t>h _ ^" p 
^UU,Cahn ~ Q2 



(74) 



All the above parton-model results have been obtained using the most general decompositions of 
the correlation matrices $ and S, and inserting them into the SIDIS hadronic tensor W^^ . There is 
an alternative approach, which relies on the helicity formalism and expresses the cross section as a 
convolution of helicity amplitudes of elementary subprocesses with partonic distribution and fragmen- 
tation functions, taking fully into account non collinear kinematics |188j . Considering for simplicity an 
unpolarised lepton beam and a spinless or unpolarised final hadron, and adopting the 7*A^ collinear 
frame, the basic factorisation formula for the SIDIS cross section in this picture is, up to order l/Q, 



1i \i Ki 



(75) 



where the A's are helicity indexes, fqi{x,k±) is the probability of finding a quark with momentum 
fraction x and transverse momentum k± inside the target nucleon, p\^ x'. helicity density matrix 

of the quark g^, D^^'^f{z, p±) is the fragmentation function of the struck quark qj into the hadron h with 
transverse momentum p± with respect to the fragmenting quark, and dcTA,, a^, ~ ^A^,Ag^.;AfAg.^A^,Ag^;A(.A' 
is the cross section of lepton-quark scattering £qi — )■ i'qf at tree level. Note that, whereas in the collinear 
case the produced hadron is constrained to have Ph± = and the entire process takes place in the 
scattering plane, the intrinsic transverse momentum of quarks introduces a non planar geometry. The 
elementary scattering amplitudes M's take into account this non collinear and out-of-plane kinematics. 
Neglecting 0{kj_/Q^) contributions, no Jacobian factors appear in eq. ( !75|) and one has x = xb, z = z^. 

Despite their apparent dissimilarity, the two parton model approaches described so far, namely the 
approach based on quark correlation matrices and eq. f l62|) and the approach based on the helicity 
formalism and eq. fl75]) . are equivalent as far as parton interactions are ignored. In other terms, all the 
leading-twist asymmetries listed in eqs. f lMlITT]) can be exactly reobtained from eq. f l75]) |189] . whereas 
at twist three the results of the two approaches are identical if one neglects the "tilde" distribution and 
fragmentation functions arising from quark-gluon correlations. 



4.1.2 TMD factorisation in QCD 

So far, we have only considered the extended parton model, i.e. the parton model incorporating intrinsic 
transverse momenta. One may wonder whether the results we have presented have any solid QCD foun- 
dation. The answer to this question is positive, at least in a particular kinematic regime. Semi-inclusive 
processes are characterised by two scales, besides the confinement scale Aqcd- the momentum transfer 
Q and the transverse momentum of the final hadron Ph± or, equivalently, the transverse momentum 
Qrp = \qnp\ of the virtual photon in the hN collinear frame. 

Extending the pioneering work on back-to-back jet production of Ref. [89], a TMD factorisation 
theorem for SIDIS and DY has been proven [521 ED I19U] . The proof is valid in the low transverse- 
momentum region Ph^iQx) ^ Q- In this framework the unpolarised SIDIS structure function is 



Figure 13: Feynman diagrams of the elementary processes contributing to SIDIS at first order in a. 



written as 

Fuu,T{^B,Zh,Q\Q^T) = ^(^l^ l^^^r I^^^t j dHr 5\kT - + h + Qt) 

a 

xH{Q^) f^ixB, kl) Dl{zh, 4) UHl) . (76) 

For simplicity the dependence of the distribution functions on = {2v-PY /v"^ and of the fragmentation 
function on Cl = {2v ■ P)^/?)^ {v and v are vectors off the light-cone) is omitted . The variables ( and 
(h serve to regulate the light-cone singularities, as explained in Section 13.31 if is a perturbative hard 
factor written as a series in powers of a^. The soft factor U arises from the radiation of soft gluons 
(of transverse momentum It) and is a matrix element of Wilson lines in the QCD vacuum. Also not 
displayed in eq. (176|) is the dependence of all quantities on the renormalisation scale fj, and on the 
soft-gluon rapidity cut-off p = a/ (2f ■ v)'^/v'^v'^. Of course, the physical observable F does not depend 
on any of these regulators. 

The generalisation of eq. ( I7B]) to the polarised structure functions, in particular to those generating 
transverse SSA's, has been proposed in Refs. |19H 1192] . The parton model expressions of Section [4. 1.11 
are recovered at tree level, i.e. C(a°), since H^^^ = 1 and U^^\It) = 5^(/t)- 



4.1.3 QCD description at high transverse momenta 

At high transverse momenta, Qt S> Aqcd, SIDIS structure functions can be described in collinear 
QCD. The azimuthal angular dependence of hadrons in leptoproduction was proposed longtime ago as 
a test of perturbative QCD |193] . In collinear factorisation, transverse momenta are generated by gluon 
radiation. At first order in the hard elementary processes shown in Fig. [12] contribute to the four 
unpolarised SIDIS structure functions Fujj and to the two double- longitudinal structure functions Fn. 
Introducing the partonic variables x and z, defined as x = Q^/2k ■ q = xb/x, z = k ■ k'/k ■ q = Zh/z, 
where k and k' are the four-momenta of the incident and fragmenting partons, respectively, and x and 
z are the usual light-cone momentum fractions, i.e. k = xP and k' = Ph/z, one has for the Fuu's at 
leading order in and leading twist [194^ 1195^ 1192] 



'uu[x,Q)= }_^e,XB / ^/ 
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n (f ) Dl (I) ClT'' + ft (f ) Dl (I) CiT'' + n (f ) Dl (I) ClT''] (77) 



and analogous formulae for the Flls. The Wilson coefficients C represent elementary cross sections 
and are listed in Ref. |192] . 

The structure function F^^'^'" encountered in Section [4:. l.lj which produces a beam-spin asymmetry 
and vanishes in the parton model, gets a non zero perturbative QCD contribution at leading twist and 
order [TO 1197] . 



On the contrary, the transversely polarised structure functions Fjjti which vanish at leading twist in 
coUinear factorisation, since there is no chirally-odd fragmentation function, emerge at twist three, as the 
result of quark-gluon correlations. Following the early work of Ref.[T7t [T8| 1198] . a twist-three coUinear 
factorisation theorem valid at large transverse momenta was proven [53l [54t [55] . In this approach the 
cross section for SIDIS with a transversely polarised target has the general form [1911 11991 1200] 



do- ~ Gf{x,x') ® do- (g) Di{z) + hi{x) (g) da' ® Ef{z,z') 



(7J 



where the first term contains a quark-gluon correlation function for the transversely polarised nucleon 
and the ordinary unpolarised fragmentation function for the final hadron, whereas the second term 
combines the transversity distribution with a twist-three fragmentation function. Let us focus on the 
first contribution (twist-three effects in the initial state). The hadronic tensor can then be written as 



/dz 
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(79) 



where the partonic tensor w^i, contributing to the transversely polarised structure functions is (see 
Fig.Mj 



w^^iP, g, P,) = Jd^'k J d^k' Tr [$a(A;, k') H,,{k, k\ q, P,)] 
In this expression is the quark-gluon correlation matrix 
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and H^i, represents the perturbatively calculable partonic hard scattering. By means of the coUinear 
expansion [M] one can get 
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and finally end up with 
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where the quark-gluon correlation matrix defined as 
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(P, ^I^(O) W-[0, r,] gF+''{r,)W-[7], ^(0|P, S) , (84) 



contains the multiparton distributions Gp^Gp, Hp, Ep introduced in Section [331 It is easy to verify 
that, due to the structure of the hadronic tensor receives contributions only from the imaginary 
part of the hard blob, arising from internal propagator poles. 

Considering for definiteness the Sivers contribution to the cross section, its explicit expression is 
[T9n[T99] 



da I 



Siv 



dx f dz 



X 



z 



Q 



— —5^-1- — 



X 



Xb 



1-^ 

Zh 



X 



dG"' (x x) 

X — (ID + Gp{x, x) ac + Gp{x, 0) ap + Gp{x, xb) 



DUz) + ... (85) 




Figure 14: General diagram contributing to SIDIS SSA's in the twist-three factorisation. 



The first two terms represent the so-called "soft-gluon pole" contribution {xg = x' — x = 0), the third 
term is the "soft-fermion pole" contribution [x' = 0), the fourth term is the "hard pole" contribution 
{x' = xb)- The dots represent the contributions of Gp and of the gluonic correlation functions. 

In the intermediate transverse-momentum region, i.e. Aqcd ^ Qr^PhA.) ^ Q^i expects that 
both the TMD and the twist-three pictures should hold. This has been explicitly verified in Refs. [56l |58] . 
The output of these important works is a set of relations that connect the T-odd TMD distributions 
(Sivers and Boer-Mulders functions) on one side, with the quark-gluon correlations on the other side. 



4.1.4 Leptoproduction of transversely polarised spin 1/2 baryons 

The leptoproduction of spin 1/2 baryons in another interesting channel to access transversity. If the 
nucleon target and the detected hadron are both transversely polarised, £ -|- A^^ — )■ + -B''' + X, the 
spin transfer between the initial and the final particle occurs in coUinear kinematics. The cross section 
integrated over the transverse momentum gets in fact a double-spin term proportional to the product 
of the transversity distribution hi and the "transversity" fragmentation function Hi. 

This doubly polarised process has the advantage of being free from the theoretical complications 
related to the transverse motion of quarks: the ordinary collinear QCD factorisation theorem applies 
and a perturbative study is possible, since we know the evolution of both hi (Section 13. 2p and Hi 
|201] . At leading order in QCD, the transverse polarisation of the produced baryon is given by 



= DMy) t:" l i , (86) 



where is the nucleon polarisation and Z)jvAr(y) = (l—y)/ (1 — y + y^/2) is the spin transfer coefficient 
of quarks. 

In the case of A hyperons, information on the spin transfer in the fragmentation process can be 
obtained from the A polarisation extracted from the angular distribution in the weak A — )■ pvr" decay. 
The transverse polarisation of A's produced in hard processes with initial transversely polarised hadrons 
was studied long time ago in Refs. |202l 1203] . and reinvestigated more recently |204l [32] . Phenomeno- 
logical analyses are presented in Refs. |205[ I206[ 1207] . The transversity fragmentation function Hi is 
measurable in e~^e~ production of transversely polarised A pairs. 

Transverse A polarisation can also be observed in SIDIS with an unpolarised target, by measuring 
asymmetries in the transverse momentum distribution of the hyperons. One contribution (T-odd in 
the final state) involves the quark density fi and the polarising fragmentation function -Dj*^. This 



mechanism for A polarisation in SIDIS was studied in Ref. |2U8j . Another contribution (T-odd in 
the initial state) involves the Boer-Mulders function hi coupled to the (unintegrated) fragmentation 
function H^^z^p^). At high P/i_l, this effect has a counterpart in the twist-three approach. The 
initial-state T-odd mechanisms for producing A polarisation have been investigated in the context of 
twist-three factorisation by Zhou, Yuan and Liang [209] . 

4.1.5 Two-hadron leptoproduction 

Another partially inclusive DIS reaction that can provide information on the transverse-spin structure of 
hadrons, and in particular on transversity, is two-particle leptoproduction from a transversely polarised 
target, i{l)+N^{P) {l')+hi{Pi)+h2{P2)+X{Px), with the two spinless final hadrons in the same jet. 
Two-hadron production in SIDIS has been proposed and studied by various authors |210[ [33 tl211U212] as 
a process probing the transverse polarisation distribution in combination with a dihadron fragmentation 
function (DiFF). The idea is to look at an angular correlation between the spin of the fragmenting 
quark and the relative transverse momentum of the hadron pair, without involving the transverse 
momenta of quarks. Integrating over the total transverse momentum of the final hadrons, one gets 
an asymmetry in the azimuthal angle between the two-hadron plane and the scattering plane. This 
asymmetry is determined by a fragmentation function usually called Hf, which does not depend on the 
intrinsic transverse motion of quarks and arises from the interference between different channels of the 
fragmentation process into the two-hadron system. Thus, all the difficulties related to non-collinearity 
are in this process avoided. 

The first authors who suggested resonance interference as a way to produce non-diagonal fragmen- 
tation matrices of quarks were Cea et al. |213] in their attempt to explain the observed transverse 
polarisation of A hyperons produced in pN interactions The unpolarised dihadron fragmentation 
functions appeared for the first time in the context of jet calculus |214[ 1215] . The extension to the 
polarised case was discussed in Refs. |210] 130] 1216] [33] and a complete classification of the DiFF's was 
given at leading twist in Ref. |211] and at twist- 3 in Ref. |217] . 

The kinematics of the process in the 7*A^ frame is shown in Fig. [151 We introduce the total 
momentum of the hadron pair Ph = Pi + P2 (with invariant mass = P^), the relative momentum 
R = (Pi — ^2)72, and the variables z = Zi + Z2 = /k^ + P2 / 1^^ = Ph / (the light-cone fraction 
of the fragmenting-quark momentum carried by the hadron pair) and C = '^R /Ph (which describes 
how the total momentum of the pair is split into the two hadrons). Rt is the transverse component of 
R with respect to Ph-, and (f)R is the azimuthal angle of -Rt in the plane orthogonal to the 7*A^ axis, 
measured with respect to the scattering plane. The azimuthal angle of the target spin vector is (f)s- 
Calling kt the transverse momentum of the fragmenting quark with respect to Ph-, the unpolarised and 
transverse-spin projections of the two-hadron fragmentation matrix A at leading twist are 

Al^^l = Di{z, C, 4. Rt^ kt-Rt) , (87) 
Al-'"^^] = jr^-TT- [4«^T, H^{z, C, 4, P2 , ^^.R^) + ei^RTj Hf{z, (, 4, P^, Kt-Rt)] ■ (88) 

Ml + M2 

These are the probabilities for an unpolarised quark and for a transversely polarised quark, respectively, 
to fragment into a hadron pair. Upon integration over Kt-, the contribution of the Collins-type DiFF 
Hi disappears, and the only remaining transverse-spin term is the one containing the (integrated) 
interference fragmentation function Hi{z,C, M^). 

It is convenient to consider the centre-of-mass frame of the two hadrons (Fig. [T5]l . where Rt = 
\Rt\ = \R\ sin 6* and 6 is the angle between the direction of the hadron emission and Ph (in the 'y*N 
frame). For two alike hadrons of mass m, one has ( = 2 {\R\/Mh) cos e and \R\ = ^^Ml-Am^. 
In terms of these variables, keeping only the unpolarised and the transversely polarised terms, the 



Figure 15: Kinematics of two-hadron leptoproduction. 



leading-twist partonic expression for the cross section of two-hadron leptoproduction reads 



dxBdydzhdcpRdcpsdM'^d cos 9 
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^E^'^ {(l-^ + y) mxB)D^^{zH,Mlcose) 



h 



Here it is z = Zh ^ P ■ Ph/P ■ Q, a relation valid modulo l/Q"^ corrections. At twist three, there 
appear extra terms in the cross section: a cos (j)R unpolarised contribution and a sin (ps transverse-spin 
contribution |217j . 

The first model for the two- pion DiFF H^''^ was presented in Ref. [33], where the phase difference 
between s and p waves was taken from vrvr phase shifts in elastic scattering. The resulting fragmentation 
function (called 6qi in Ref. [33j), changes sign around the p mass. 

In a more recent model |218j the fragmentation functions are expanded in Legendre polynomials of 
cos 6 keeping only the first few terms, corresponding to the lowest values of relative orbital momentum. 
This truncation is expected to be legitimate for not very large M^. Thus one can write 



(remember that multiplies a sin 6* factor in the cross section). The interpretation of these terms, 
signaled by their superscripts, is the following |218j : D° is a diagonal component, receiving contributions 
from s and p waves of the dihadron system separately (the "background" ) ; Dl^ and H^''^^ originate from 
the interference of a s wave and a p wave; D^^ and Hf'^^ arise from the interference of two p waves. 
The main channels contributing to the fragmentation of a quark q into a tt^tt" pair are: 1) incoherent 
fragmentation, q — t- tt+tt^X; 2) fragmentation via a p resonance, q — )■ pX — t- tt+tt^X; 3) fragmentation 
via a uj resonance decaying into three pions, q — )■ uX — )■ tt+tt^tt^X. Pions in channel 1 are expected to 
be mostly produced in s wave; pions in channel 2 come from the two-body decay of a vector meson and 
are in a relative p wave; pions in channel 3 are prevalently in p wave, but a fraction of them may also 
be in s wave. The functions and Hf'"^ arise from the interference of channels 1-2 and 1-3 [21 9J. 

A model based on a more sophisticated analysis of the fragmentation channels |219j predicts a 
completely different behaviour for H^''^^, with a peak at the p mass and a broader maximum at the to 
mass. Its size is about 30 % of the unpolarised fragmentation function as shown in Fig. [161 

A different definition of the relative transverse momentum is proposed by Artru [28j, who uses the 
vector = {z2Pi± — ZiP2±){zi + Z2), which is perpendicular to the 7*X axis, and its azimuthal angle 



D,{zh, Ml cos 6) = D'l{zh,M^) + Dl''{zH,M^) cos^^ + (^^, M^2)i(3 cos' 
H^izH, Ml cos 6) = H^'^'iz,, Ml + H^'^'izh, Ml cos 6 



,^9-1) 



(90) 
(91) 




0.4 0.6 0.8 1 1.2 0.2 0.3 0.4 0.5 0.6 0.7 0.8 



Mh (GeV) z 

Figure 16: Model prediction for the ratio {—\R\Hf'^'^/{MhD°) as a function of (left) and z (right). 
The dotted lines represent the positivity bounds. 
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Figure 17: Left: Configuration of the process e+e^ — >■ hih2X in the jet frame, used for the cos(0i + ^2) 
reconstruction of the Collins asymmetry. Right: The same process in the Gottfried- Jackson frame, used 
for the the cos 20o reconstruction of the Collins asymmetry. 



(j)r. The advantage of (f)r is that it is by construction invariant with respect to boosts along the 7*A^ 
direction, which is not the case of (pR although the two angles are the same in the 7*A^ frame. Using 
the recursive fragmentation string model (see Section 14.1 .51) the "joint pt spectrum" of the first and the 
second rank hadrons has been calculated |177j . By suitably integrating the spectrum, expressions for 
both the Collins FF and the dihadron FF in principle may be obtained. 

4.2 Inclusive production of hadron pairs in e+e^ annihilation 

An independent source of information on the Collins fragmentation function is inclusive two-hadron 
production in electron-positron collisions, e"*" + e~ hi + h2 + X , with the two hadrons (typically 
pions) in different hemispheres. We know that the Collins function produces a cos0 modulation, 
where is the azimuthal angle between the plane containing the quark and the hadron momenta, and 
the plane normal to Sg. Considering a single jet in e~^e~ hadron production, the Collins modulation 
would average to zero in a large event sample. Thus, in e~^e~ annihilation the Collins effect can only 
be observed in the combination of two fragmenting processes of a quark and an antiquark, resulting in 
the product of two Collins functions with an overall modulation of the type cos(0i + ^2), where 0i and 
02 are the azimuthal angles of the final hadrons around the quark-antiquark axis, with respect to the 
e~^e~ — )■ qq scattering plane. 

Two-hadron production in e+e" collisions was studied in Refs. [22D1 [Ml [2221 [2231 [221 [22S]. The 
tree-level differential cross section in the jet frame with respect to the quark-antiquark direction (Fig.[T71 



left) reads 



em 
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+ sin^ 9 cos(0i + 02) H^'^'^z^) i/^"^"(^2)} , (92) 

where dQ = d cos ^ d(j)e {9 is the angle between the lepton axis and the qq axis, in the qq centre-of- 
mass frame, whereas (pe gives the orientation of the scattering plane around the qq axis), and we have 
introduced the one-dimensional moments 



F{z,p't)- (93) 



Eq. (!92|) and the following results refer to the case of photon-mediated e~^e~ annihilation. Z production 
and 'y*Z interference effects have been investigated in Refs. \222\ 1225] . 

From an experimental point of view, the quark-antiquark direction is not directly accessible, and is 
approximated by the dijet thrust axis n, defined by 

T = max ^J%^ , (94) 

where the sum is over all detected particles. The resulting cos(0i + ^2) asymmetry is given by 

sin^^ E.el(Ht^^^z^)Ht^^\z2)) 
^MO, ^1, ^2) = ^ — • (95) 



1 + cos^^ ^^el[Dfiz,)Dr) 



We will refer to this method of extracting the Collins asymmetry by measuring azimuthal distributions 
around the thrust axis as to the "cos(0i + ^2) method" [2211 SI I^Sl [225] . 

There is another way of reconstructing the asymmetry, the so-called "cos20o method", which is 
based on a different geometry and does not require the knowledge of the thrust axis. In this case one 
measures the hadron yields as a function of (pQ, the angle between the plane containing the momentum 
of hadron 2 and the leptons, and the plane defined by the two hadron momenta |222[ I223[ 1225] (this 
frame, shown in Fig. [T7] (right), is similar to the Gottfried- Jackson frame in Drell-Yan processes j227j ). 
The corresponding cross section is 
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where h = Pi±/\Pi±\ and kit, K2T are the transverse momenta of the two fragmenting quarks with 
respect to the hadron directions. The cos 20o asymmetry reads 
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Electron-positron scattering can also allow accessing the fragmentation function Hi of transversely 
polarised baryons [271 12201 1228] . In the case of A's, the specific process is back-to-back AA inclusive 



Figure 18: Geometry of two hadron-pair production in e'^e collisions. 



production, e~^e — )■ A AX, with the hyperon and the anti-hyperon decaying into pir and pvr+, respec- 
tively. It was shown in Ref. |220] that the cross section of this process contains an azimuthal modulation 

proportional of the type sin^cos(0 + (f))H{^{zi)H^{z2) oc H^{zi)H^{z2), where 9 is the angle between 
the collision axis and the qq axis in the jet frame (Fig. [T71 left) and (p, (f) are the azimuthal angles of 
the proton and antiproton in the same frame (with the z direction given by the qq jet and the x axis in 
the plane of the beams and jets. The cos(0 + 0) asymmetry is determined by measuring the difference 
between the number of pp pairs on the same side of the scattering plane and the number of pairs on 
opposite sides. 

Finally, the interference fragmentation function can be extracted from the production of two 
hadron pairs in electron-positron annihilation: e^e~ — > {hih2){h[h'2)X , where the particles in brackets 
belong to two back-to-back jets [216] . The (complicated) geometry of this process is shown in Fig. [181 
The observable quantity is the angular correlation of the production planes, expressed by the so-called 
Artru-CoUins asymmetry. The kinematics of the process is described by doubling the variables intro- 
duced in Section |4.1.5[ If we call 0r and (pR the azimuthal angles of the transverse relative momenta Rt 
and Rt of the two hadron pairs, the Artru-Collins azimuthal asymmetry is the cos{(j)R + (f)fi) correlation. 

In e~^e~ annihilation, the interference fragmentation functions and the Collins function are typically 
probed at much larger scales compared to SIDIS. However, the evolution of H^, differently from that 
of H^, is known |229] . Therefore, a consistent combined analysis of dihadron production in e~^e~ 
annihilation and SIDIS is possible and may provide an alternative way to extract the transversity 
distributions |230j . 

4.3 Drell-Yan production 

Drell-Yan (DY) dilepton production with various polarisations of the two particles in the initial state 
is a very rich source of knowledge on the hadronic structure. The main advantage of this class of 
reactions is that they do not involve fragmentation functions, but only parton distributions. However, 
unless one considers antiproton-proton scattering, or pion-proton scattering, DY processes necessarily 
involve sea x valence products. This means that, while they provide direct information about antiquark 
distributions, which are less determined in SIDIS, their asymmetries are generally small. 

In principle, DY production with two transversely polarised hadrons is the cleanest reaction for 
studying the transversity distribution hi{x) and the pioneering works of Ref. and ^231j were indeed 
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Figure 19: The Collins-Soper frame. The z-axis bisects the angle between P2 and —Pi, the momenta 
of the two initial state hadrons. 



devoted to this process. However, in order to observe sizable double-spin asymmetries and extract hi, 
we probably have to wait for a new generation of experiments with polarised antiprotons |232j . On the 
other hand, unpolarised and singly-polarised DY processes can probe a large variety of TMD's related 
to transverse spin, and are now attracting a wide theoretical and experimental interest. 

4.3.1 Kinematics and observables 

Drell-Yan lepton-pair production is the process Ai{Pi) + ^42(^2) — ^ ^"""(O + ^~(^') + ^5 where Ai and 
A2 are hadrons and X is an undetected system. The center-of-mass energy squared of this reaction is 
s = {Pi + P2)^ — 2 P1P2, having neglected in the approximate equality the hadron masses Mi and M2. 
The lepton pair originates from a virtual photon with four-momentum q = I + I'. In contrast to DIS, 
g is a time-like vector: Q"^ = > 0, and the invariant mass of the lepton pair coincides with Q^. 
The deep inelastic limit corresponds to Q^, s — 00, with r = Q"^ / s fixed and finite. 

The DY cross section is usually expressed in a dilepton center-of-mass frame and can be written as 



dV 



em 



where L^y is the familiar leptonic tensor and W^^ is the DY hadronic tensor. Among the infinite 
dilepton cm. frames, related to each other by a rotation, the most often used is the Collins-Soper (CS) 
frame |233] . characterised by a z axis that bisects the angle between P2 and —Pi as shown in Fig. [T3 
Another common dilepton cm. frame is the Gottfried-Jackson frame [227] . where the z axis coincides 
with the direction of one of the colliding hadrons. 

In the unpolarised case the DY hadronic tensor contains four independent structure functions |234 
Using the classification of Ref. |235] the cross-section becomes 



= ^ 1(1 + cos^ 9) WIjjj + sin^ 9 W^^ + sin 29 cos W^°^^ + sin^ 9 cos 20 W^°^''^\ . (99) 

The double subscript refers to the polarisation states of the two colliding hadrons: U = unpolarised, 
L = longitudinally polarised, T = transversely polarised. In literature [234j . the structure functions 
^uu^^uuy^uu^ ^^uu'^'^ ^Iso called (apart from a common factor), Wt,Wl,Wa,W^a, respec- 
tively. The angular distribution of leptons is often parametrised as 

^ = ^— ^ f 1 + A cos^ 9 + fx sin 29 cos 20 + ^ sin^ 9 cos 20^ . (100) 
A/tot di I 4:71 A + 3 V 2 / 

The three quantities A, /x, and u are related to W^^, ^uu, W^^J"^, and W^^^'^ by 



The so-called Lam-Tung relation [231 ESZ] 

A + 2z/ 



1 



(102) 



which corresponds to W^JJ = 2W^^'^'^, is valid at order as in collinear QCD |238] (see below) and is 
slightly violated at order [239]. 

In the polarised case, a complete analysis of the DY hadronic tensor is contained in Ref. |235j . where 
it is shown that in single-polarised DY, first studied in [23 Ij . the number of independent structure 
functions is 16, whereas the double-polarised hadronic tensor contains 28 structure functions, so that 
altogether the number of independent structure functions in DY is 48. The full cross section can be found 
in Ref. |235] . Limiting ourselves to unpolarised, single-transverse and double-transverse contributions, 
the cross section reads 
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(103) 



05j and 052 are the azimuthal angles of the spin vectors of hadrons A and B, respectively. This angular 
structure is valid in any dilepton cm. frame, but the numerical values of the structure functions are 
frame- dep endent . 

In the parton model and at leading order in QCD the two invariants xi = Q^/2Pi-q and X2 = 
Q^/2P2-q can be interpreted as the fractions of the longitudinal momenta of the hadrons A and B 
carried by the quark and the antiquark that annihilate into the virtual photon. In the cm. frame 
of the two colliding hadrons, which is the most convenient frame to study the partonic structure of 
the hadronic tensor, the photon momentum can be parametrised as = {xiPi , X2P2 Qt) and 
acquires a transverse component Qt. Neglecting terms of order 1/Q^, one has Q'^/xiX2S = 1, that is 
r = Q'^/s = X1X2. The structure functions in eq. fll03p can be expressed in terms of the four variables 
xi,X2, Qt = \qT\,Q. 

Other variables customarily used are the rapidity of the virtual photon, y = \ ln(g+/g^) = | \n.{xi/x2) 
and the Feynman variable xp = 25^/^/5 = xi — X2. In a dilepton cm. frame, y = |(1 + cos^^). The 
relation between (xi,X2) and (t, ?/) is Xi = y/re^ , X2 = y/re'^. The DY cross-section can be variously 
reexpressed in terms of these variables: 

dV 2 dV ^ dV 



d'^g s dxidx2d^qT dydQ'^ d'^qr 
4.3.2 DY asymmetries in the TMD approach 
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In the parton model, calling ki and k2 the momenta of the quark (or antiquark) coming from hadron 
Ai and A2 respectively, the hadronic tensor shown in Fig. [20] is 

(105) 
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Figure 20: The parton-model diagram for the DY hadronic tensor. 



Here $ is the quark correlation matrix for hadron Ai, $ is the antiquark correlation matrix for hadron 
A2, and the factor 1/3 has been added since in $ and $ summations over colours are implicit. It is 
understood that, in order to obtain the complete expression of the hadronic tensor, one must add to 
(11051) a term with replaced by $(A;2) and $(/c2) replaced by $(A;i), which accounts for the case 

where a quark is extracted from A2 and an antiquark is extracted from Ai. In the following formulae 
we shall denote this term symbolically by [1 2]. 

Introducing the longitudinal momentum fractions ^1 = / and ^2 = ^2 ■, and working out 
the delta function of four-momentum conservation, one finds 1^1 = ^2 = 2:2, and the hadronic tensor 
becomes 
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k2T - qT)Tr [<l>(xi, fcir) 7^ $(^2, k2T)Y] + [1^2]. (106) 



Inserting here the explicit partonic expressions of the quark correlators, it is not difficult to get the 
parton model expressions of the DY structure functions [HI I240[ 1235] . Only 24 of these 48 structure 
functions are non vanishing at leading twist. The most relevant ones are, in the Collins-Soper frame 
(but the Gottfried- Jackson expressions differ from these only by subleading terms 0{Qt/Q)) 
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where h = Qt/Qt and C denotes the transverse-momentum convolution of eq. fl63|l with the addition of 



the [1 -H- 2] term. In eq. flllOp we defined the combinations WtijIut)'^^^^'^^^ 
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which correspond to the angular modulations indicated by 



Eqs. f llOTHlll]) contain a series of interesting results. First of all, the Boer-Mulders function hj- 
generates a cos 20 asymmetry in unpolarised DY. At leading twist, the u parameter is given by 
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Since I' ^ and A = 1, the Lam- Tung relation is violated, and this is one of the remarkable consequences 
of the intrinsic transverse motion of quarks. 

Concerning singly-polarised Drell-Yan processes, the Boer-Mulders function combines with the 
transversity distribution in the sm{2(f) — (psi), or sm{2(j) — ^Sj)? asymmetry, whereas the Sivers function 
is probed via the F^^ (or F^j.) structure function associated with the sin^s^ (or sin ) asymmetry. 
Note that upon integration over the lepton angles only the Sivers asymmetry is non vanishing. 



4.3.3 DY double transverse asymmetries 

With two transversely polarised colliding hadrons, the cos (20 — ~4'S2) term provides a direct access 
to transversity. Inserting eqs. (1107^ lllip into eq. fll03p and integrating the cross section over qx, one 
gets 
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+ ^ir^2r sin^ 6 cos(20 - 05, - 05j/ii(xi)/i^(x2)l +[1^2]. (113) 



This parton-model expression can be generalised to QCD by resorting to the collinear factorisation 
theorem, which for the polarised DY process reads [241] 

da = E E /"d6piV,^/a(6,/^')pl'\Ja(6,/i')da,,v,A.A^(g^/i^«.(/x2) (114) 

where ^1 and ^2 are the momentum fractions of the quark (from hadron Ai) and antiquark (from A2), 
p^^^ and p^^^ are the quark and antiquark spin density matrices, dax^y^x^y^ is the cross-section matrix 
of the elementary subprocesses in the quark and antiquark helicity space, fi is the factorisation scale. 
At leading order, i.e. C(a°), the only contributing subprocess is qq -> and ^1 = Xi, ^2 = ^2. In 
the transversely polarised case, one reobtains eq. flllSp . except that all distribution functions acquire a 
Q^-dependence. Thus the LO double transverse asymmetry is 

~ Eaeimx,,Q^)mx2,Q-) + [1 o 2] ' ^^^'^ 

where 

sin^ 6' 

= ^ 2^ cos(20 - 0s, - 05j , (116) 

1 + cos^ u 

is the elementary double-spin asymmetry for qq — )■ We see that a measurement of Af^^ would 

directly provide the product of quark and antiquark transversity distributions, with no mixing with other 
unknown quantities. At next-to-leading order (NLO) the DY transverse cross section gets contributions 
from virtual-gluon (vertex and self-energy) corrections and real-gluon emission, which were calculated 
by several authors with different methods [242^ 1243^ 1244^ [ST] . The NLO double transverse asymmetry 
was investigated in Refs. [243 1213 I 



4.3.4 DY azimuthal and spin asymmetries in QCD 



As in the case of SIDIS, perturbative gluon radiation can generate a non-zero transverse momentum Qt- 
For instance, tlie contribution of the quark-antiquark annihilation process qq — )■ 7*5' to the unpolarised 
angular distribution in the Collins-Soper frame is |238] 



AT,., d^^ - 16. V + + 2 ^ 20 + . . . J , (117) 

where we have omitted the sin 29 cos (f) term which is the only one depending on the quark and antiquark 
distributions. From ( I117p one gets 

and the Lam- Tung relation is fulfilled. The contribution of the qg — )■ 7*g is more complicated, but also 
satisfies the Lam- Tung relation, which holds for the complete leading-order cross section. 

The perturbative QCD approach to the DY angular distribution sketched above holds for Qt ~ Q- 
At small Qti large logarithms of the form \\i{Q'^ /Q\) appear, which must be resummed. This is done 
in the space conjugate to Qt and gives rise to a Sudakov form factor, according to the Collins-Soper 
procedure [89]. The Sudakov resummation for the structure functions of unpolarised DY production, 
including those related to azimuthal asymmetries, has been studied in Refs. |249[I250] and the Lam- Tung 
relation is found to be unaffected by the resummation. 

The perturbatively generated cos 20 asymmetry is suppressed as Qf^/Q^ at small Qt- A further 
contribution to this asymmetry can arise in the twist-three approach from the product of two quark- 
gluon correlation functions Ep, one associated with the quark from hadron Ai, the other with the 
antiquark from hadron A2 |251j . 

In the singly-polarised DY case the situation is again analogous to SIDIS. At large Qt ~ Q ^ Aqcd, 
a DY single-spin asymmetry is generated by the Gp quark-gluon correlator of the polarised hadron. A 
smooth transition from this twist-three mechanism to the Sivers effect occurs in the intermediate region 
Aqcd ^ Qt ^ Q, where both the higher-twist and the TMD factorisations apply ^J. Another SSA's 
arises from the chirally-odd quark-gluon correlation function Ep of the unpolarised hadron coupled 
to the transversity distribution of the transversely polarised hadron, but in the low Qt limit this 
contribution vanishes (after integration over the lepton angles only the Sivers asymmetry survives). 

4.4 Inclusive hadroproduction 

We finally discuss a third class of reactions that probe the transverse-spin and transverse-momentum 
structure of hadrons: inclusive hadroproduction with one transversely polarised hadron in the initial 
state, that + B — )■ h + X, where an unpolarised (or spinless) hadron h is produced with a 

transverse momentum Pt with respect to the collision axis. An interesting variation of this process is 
the production of a transversely polarised hadron, i.e., a A hyperon, from unpolarised hadron-hadron 
scattering, that is A + B ^ + X (Section I4.4.3p . We will limit ourselves to a brief description of 
these reactions, referring the reader for more detail to some reviews [HSl 1252] and to the original papers. 

We first consider hadroproduction with a transversely polarised colliding particle. The measured 
quantity is the single-spin asymmetry 

do"^ — do-'' 

with the cross sections usually expressed as functions of Pp and of the Feynman variable xp = 2Pi^j ^ 
[Pp being the longitudinal momentum of the produced hadron). In terms of the scattering angle 6, 



Feynman's x can be written as xp = 2P'r/ \/sidA\9. Another often used variable is the pseudorapidity 
7] = -lntan(6'/2). 

According to the QCD factorisation theorem |253[ 1241] the differential cross-section for hadropro- 
duction at large Pt can be formally written as 

d^ = J2 Yl PlA'JaM®fb{Xb)^da,^y^x.y^^V'lil,^{z). (120) 
abc XaX'^XcX'^ 

Here fa [fb) is the distribution of parton a (b) inside the hadron A {B), p^^y ^P^^ density matrix 

of parton a, v'l^y is the fragmentation matrix of parton c into hadron h, and da is the (perturbatively 
calculable) cross-section of the elementary process a + 6 — )• c + . . . (a two-body scattering, a + b ^ c + d, 
at lowest order). 

If the produced hadron is unpolarised, or spinless, only the diagonal elements of T>'l^J^, are non zero, 

i.e. T>^^y ~ ^AcA'c D^^'^, where D^^'^ is the unpolarised fragmentation function. Together with helicity 
conservation in the partonic subprocess, this implies = X'^. Therefore, the cross section fll20p carries 
no dependence on the spin of hadron A and all single-spin asymmetries vanish [16]. In order to escape 
such a conclusion one must consider either the intrinsic transverse motion of quarks [511 123 ED], or 
higher-twist effects [TTJ [TSl 1198] [T^ IM]- In the former case, one can probe a number of distribution 
and fragmentation functions, including the transversity distribution (transversely polarised quarks in 
hadron A'^), the Sivers function (unpolarised quarks in hadron A^), the Boer- Mulders function (trans- 
versely polarised quarks in hadron B), the Collins function (transversely polarised quarks fragmenting 
into hadron h). The twist-three single-spin asymmetries involve various quark-gluon correlators, either 
in the initial state (distribution functions), or in the final state (fragmentation functions). The main 
problem is that all TMD or twist-three contributions mix up in a single observable, A^, which makes 
the physical interpretation of the results quite unclear. 

4.4.1 Hadroproduction in the extended parton model 

When the intrinsic transverse motion of quarks is taken into account, the QCD factorisation theorem 
for inclusive hadroproduction is not proven, and actually is known to be explicitly violated in some 
cases |254] 1104] . Nevertheless, one can write a non-collinear factorisation formula in the context of the 
extended parton model, with a tree-level elementary kernel. One must obviously recall that: i) the gen- 
eralisation of this kernel to higher order in as is not a legitimate procedure; ii) the transverse-momentum 
dependent distribution and fragmentation functions appearing in the hadroproduction factorisation for- 
mula are not guaranteed to be universal quantities, i.e., to be the same functions as in other processes. 

The extended-parton model formula generalising eq. fll20p is (we consider the production of a spinless 
hadron) ^ES\ 

d(^ = Yl Yl PXaX'Ja{Xa,kTa) ® fb{Xb,kTb) ^daXaX'^X,X',^®V^iy^{z,PT) . (121) 
abed XaX'^XcX'^ 

where the convolutions ® are now not only on the longitudinal momentum fractions Xa,Xb,z, but also 
on the transverse momenta kpa, kTb,PT- Note that even though h is unpolarised, its pT-dependent frag- 
mentation matrix v'l^J^, is non diagonal. The elementary cross sections have the structure daA^A^^A^A^ ~ 

J2xi,Xa ^-^AcAd.AaAfc^v Ad,A' Afc- ^he amplitudes M refer to the elementary subprocess ab cd (remember 
that we are considering the tree level only). A natural reference frame is the center-of-mass frame 
of the colliding hadrons. The collision axis forms with the direction of the produced hadron a plane, 
that we call the hadronic plane. The tricky point about eq. f ll21|) is that, due to intrinsic transverse 
momenta, the partonic scattering does not take place in the hadronic plane. This non-planar geometry 



gives rise to some non-trivial phases in the distribution and fragmentation matrices. Also, the ampli- 
tudes M appearing in eq. fll2ip must be Lorentz transformed to the canonical amplitudes M° defined 
in the partonic center-of-mass frame, an operation which introduces further phases. This complicated 
structure has been fully worked out in Ref. |188j . where all the details concerning the kinematics and 
the scattering amplitudes can be found. Here we limit ourselves to quoting some general results. The 
contribution to the transverse single-spin asymmetry from the qq — qq subprocess schematically reads 

dAa,,^y, ~ f^rl^ ® fl ® dAa' ®Dl + h1® ® dAa" ® H^' 

+K ® ® dAa'" ®D1 + ® hi^ ® dAa"" ® H^" . (122) 

One recognises the Sivers effect (first term), the Collins effect (second term), the Boer-Mulders effect 
(third term) and a mixed effect (fourth term). The other contributions, qq — > gg, qg — qg, qg — )■ gq, 
gq — )■ gq, gq — )■ qg, gg — )■ qq, gg — )■ gg, are explicitly given in Ref. jl88j . They contain, besides the 
distributions and fragmentation functions of linearly polarised gluons, the gluon Sivers function f^. 
For instance, the gq — >■ gq and gg — )■ qq contributions are 

dAagq-,gg ~ /{f' ® ^ dAff/ ®Dl + ... , dAagg^gg ~ /{^/ ^ /f ^ dA^// ^ + . . . . (123) 

Processes that select these terms and allow accessing the gluon Sivers function are D meson production 
(which is dominated by the gg — )■ cc channel) [255j and pion production at midrapidity (which probes 
in the RHIC kinematics the small Xa region and thus proceeds predominantly via gluonic channels) 
[256] . Other reactions probing the Sivers functions of quarks and/or gluons without contributions from 
the fragmentation sector are prompt-photon production + B ^ ■y + X \257\ \25S\ 1259] . photon-jet 
production + B + jet + X |238l [260] . back-to-back dijet production + B ^ jet^ + jeta + X 
t261J. 

On the contrary, the Collins effect can be singled out by studying asymmetric azimuthal correlation 
of hadrons inside a jet, that is A^ + B ^ iet + X ^ h + X [262] . 



4.4.2 Single-spin asymmetries at twist three 

As pointed out in Ref. [TT j 1198] . non- vanishing single-spin asymmetries can be obtained in perturbative 
QCD at higher-twist level. A twist-three factorisation theorem was proven for direct photon production 
[531IM] and hadron production [55]. This work has been extended to cover the chirally-odd contributions 
[263[ 12641 1265] . Here we limit ourselves to quoting the main general results of these works. The twist- 
three phenomenological studies of data is treated in Section 15.51 

At twist three the hadropro duct ion cross section is formally given by 

da = {^^(a^a, <) ® f'lixk) ® da' ® D'I'^z) + hl{xa) ® E'p{xb, x',) ® da" ® D'I'^z) 

abc 

+ K{xa) ® fl{x,) ® da'" ® E'p{z, z')^ , (124) 

where Gp^Xa^x'^ and Ep^Xa^x'^) are the quark-gluon correlation functions introduced in Section [331 
Ep[z,z') is a quark-gluon correlator in the fragmentation process and da', da" and da"' are cross- 
sections of hard partonic subprocesses. The first term in fll24p corresponds to the chirally-even mecha- 
nism considered by Qiu and Sterman [55]. The second term is the initial-state chirally-odd contribution 
analysed in Ref. [263] . The third term is the final-state contribution studied in Ref. [266] . The details 
and the elementary cross-sections can be found in the original papers. 



4.4.3 A production 

The origin of the large transverse polarisation of hyperons measured since the 70's [H |5] in high-energy 
unpolarised hadron-hadron scattering is a longstanding problem (for reviews see Refs. [2671 1268] ). A 
transverse-momentum mechanism able to produce sizable asymmetries in A + B + X involves 

the polarising fragmentation function D^j, introduced in Section I3.8[ which describes the fragmentation 
of an unpolarised quark in a transversely polarised hadron |174j . In the twist-three approach the 
A polarisation is generated by the chirally-odd spin-independent quark-gluon correlation function Ep 

While in A + B + X the A polarisation must vanish at pseudorapidity rj = for symmetry 

reasons, in A-|-jet production no such constraint exists. The process A + B — )■ jet-|-jet-|-X — A-|-jet+X 
has been proposed as an alternative way of probing through the correlation between the transverse 
momentum and spin of the A with respect to the dijet axis [269] . 

4.5 Other processes 

We conclude our discussion of the transverse-spin effects in hard processes by listing a series of reactions 
that have been proposed as sources of information on transversity and TMD's. 

- Hadron production in transversely polarised lepton-proton scattering: i + p'^ ^ h + X. Note that 
the final lepton is not detected, so this process is similar to p + h + X. Its transverse SSA 
has been calculated in the twist-three factorisation approach [265j and in the extended parton 
model [270] . 

- Dilepton photoproduction: 7 + A^''^ —)■£"'"+ £~ + X. It has been shown [271] that the transverse 
SSA of this reaction involves the transversity distribution multiplied by the chiral-odd distribution 
amplitude of the photon. 

- Exclusive vr electroproduction: e + — )■ e' + 7r° + p'. Using a model for the GPD's and relating 
Ht{x, ^, t) to the tensor charge, the authors of Ref. [272] show that the transverse-spin asymmetry 
of this process can provide information on the tensor charge 6u. 

- Photo- and electroproduction of two vector mesons: 7^*-' + A^^ — ?■ pi + P2 + A^'- If one of the 
mesons is transversely polarised this reaction probes the transversity GPD HT{x,^,t) [273] . 

5 Experimental results and phenomenological analyses 

In the last decade, many transverse-spin effects have been measured in SIDIS on transversely polarised 
targets mainly by the HERMES and COMPASS Collaborations, in hadron-hadron scattering by the 
RHIC spin experiments, and in unpolarised Drell-Yan processes at Fermilab. 

In this section, we review some of the recent experimental findings together with their phenomenolog- 
ical interpretation. The selected data are organised according to the physics information they provide. 
We start with the measurements aiming to access the transversity distribution, including the related 
measurements which are being performed in e^e~ collisions. Then we describe the experimental re- 
sults related to the T-odd TMD's (Sivers and Boer-Mulders function) and we conclude this part with 
a brief discussion of some measurements involving the T-even TMD's and higher-twist PDF's. The 
RHIC hadroproduction results cannot easily be fitted in this scheme, and are presented in a separate 
subsection. 

Since most of the results currently used to access transversity and TMD's come from the SIDIS 
experiments, we feel useful to give in section 15.11 some details on the kinematical ranges of the present 



SIDIS experiments and on the analyses these experiments are doing, which are essentially common to 
all the SSA's extraction. 



Note: Although one should in principle distinguish between the longitudinal momentum fraction 
X = / and the Bjorken variable xb = jlP ■ q-, we have seen that they coincide as far as l/Q^ 
corrections are neglected. For the sake of simplicity, in this section we ignore this distinction and 
write the distributions as functions of x. Analogously, we take z = Pf^/k~ to be the the argument of 
fragmentation functions. In the data plots, obviously x and xb always stay for xb and z and Zh for 
Zh = P- Ph/P ■ q. 

Other notations: = I^tI is the transverse momentum of the initial quark and pt = \pt\ is the 
transverse momentum of the produced hadron with respect to the fragmenting quark. In SIDIS Ph± = 
\Ph± \ is the transverse momentum of the final hadron with respect to the 7*A^ axis. In hadroproduction 
Prp = \Pt\ is the transverse momentum of the final hadron with respect to the collision axis. 

5.1 SIDIS kinematics and SSA extraction 

The SIDIS events are usually identified with standard cuts, with some differences for the different 
channels and for the various experiments due to the different beam energies and thus to the different 
kinematical domain. The DIS events are selected requiring the photon virtuality to be larger than 
1 GeV^. The fractional energy y transfered from the beam lepton to the virtual photon has to be 
larger than 0.1, to remove events affected by poor energy resolution, and smaller than 0.9 (or 0.95), to 
avoid the region most affected by radiative corrections. A minimum value of invariant mass of the final 
hadronic state ~ 2 GeV is also needed to exclude the resonance region. Typically, values of W"^ 
larger than 4, 10 and 25 GeV^ are required in the data analyses of the JLab, HERMES and COMPASS 
experiments respectively. 

The variables x, and W, for the selected events, cover ranges which strongly depend on the 
lepton beam energy. Fig. |2T] shows the regions of the (x, Q^) plane kinematically accessible with lepton 
beams of 160 GeV, 27.5 GeV and 6 GeV momenta, corresponding to the COMPASS, HERMES and 
JLab experiments respectively. In the COMPASS experiment the x range is between 0.004 and 0.3, 
where the upper limit is given by the low luminosity; for the HERMES experiment 0.02 < x < 0.4, 
while the JLab experiments can presently measure with high precision at x > 0.1, in the valence region. 
The average values are also different and there is a strong x — correlation. At x ~ 0.1 the mean 

value is 6.4 GeV^ at COMPASS and about 2.5 GeV^ at HERMES, while at x ~ 0.3 the values 
are about 20 GeV^ and 6.2 GeV^ respectively. In the overlap region the ratio of the mean values 
measured in the two experiments goes from 2 to 3 with increasing x, in spite of the similar mean values 
when integrating over the whole x range. The W"^ values are between 25 and 200 GeV^ for COMPASS, 
between 10 and 50 GeV^ for HERMES, and below 10 GeV^ for the present JLab experiments. The 
differences in the covered kinematical regions make the experiments complementary, and, all together, 
they guarantee a very good coverage of the phase space. 

In addition to the requirements on the inclusive DIS variables, in the data analysis cuts on the 
energy final state hadrons are applied, which also depend on the experiment and on the physics channel 
under consideration. Particle identification implies momenta above the RICH thresholds, which depend 
on the detector used in the experiment. In the single hadron analyses, the relative energy z of each 
hadron has to be between 0.2 and 0.8. The upper limit, usually not required in the COMPASS analyses, 
is chosen to reject exclusively produced hadrons. The lower limit is used to select hadrons from the 
current fragmentation region. To do that, a selection based on the hadron rapidity should be applied, 
or, equivalently, the so-called Berger criterion should be fulfilled. This criterion |274j has been tuned 
on unpolarised SIDIS data and allows to relax the request based on the W value alone by asking for 
each hadron a large enough z. Thus, if > 7.4 GeV, all hadrons belong to the current fragmentation. 
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Figure 21: The x — DIS regions with lepton beams of 160 GeV, 27.5 GeV and 6 GeV momenta (left 
to right). At large x, the region is limited by > 25 GeV^ or y > 0.1, > 10 GeV^ and > 4 
GeV^ for the three beam momenta respectively, while at low x it is limited by the requirement y < 0.9. 
The dashed curves give the — x correlations for equal to 100, 25 and 4 GeV^. 



Going down to ~ 5 GeV only hadrons with z > 0.2 belong to this region, while if ~ 3 GeV one 
should cut at 2 > 0.5 to be safe. 

The azimuthal asymmetries in SIDIS introduced in Section 14.11 are defined as 

A^^P = Fi^^^/Fuu, (125) 

where /($) is a trigonometric function of a linear combination $ of (pn and (ps, and XY refer to the 
beam and target polarisations (U, L, or T). The asymmetries can be extracted from the measured 
distribution of the final state hadrons in the relevant azimuthal angle $. In the following, the methods 
used to evaluate the transverse SSA's in the HERMES and in the COMPASS experiments will briefly 
described. 

In principle, since all the trigonometric functions appearing in the cross-section are orthogonal, 
the amplitudes of the azimuthal modulations (the so-called "raw asymmetry" a) can be obtained as 
twice the mean value of /($), or by fitting the azimuthal distribution with a function of the type 
F($) = const ■ [1 + a ■ /($)]. In practice this procedure requires to correct the azimuthal distribution 
for possible acceptance effects by means of Monte Carlo simulations. This can be avoided in the case 
of the SSA's by collecting data with two opposite spin orientations, indicated with "+" and "-" in the 
following. If the acceptance and the detector efficiencies are the same for the two sets of data, no Monte 
Carlo correction is needed when fitting the function -F($) on the quantities 

N+(^) - rN-(^) 
^ i^j^ — .^26) 

Here A^^ are the numbers of events in a given $ bin, r is the normalisation factor between the two sets 
of data, and $ is always measured assuming the same orientation of the target polarisation for both sets 



of data. In the case of the HERMES experiment, the target spin orientation is flipped every second, so 
that acceptance and overall efficiencies can safely be assumed to be the same for the two sets of data. In 
the case of the COMPASS experiment, in which the target polarisation can be reversed typically only 
after a few days of data taking, the target is divided in cells with opposite polarisation directions. This 
allows to minimise the possible systematic effects due to acceptance variations by using, instead of A{^) 
the so-called "ratio product" quantities [13] which combine the number of events from the different cells 
and with the different target polarisation orientation. These quantities do not dependent on the beam 
flux and on the acceptances, under the assumption that the relative variations are the same for all the 
target cells, and they have a very simple expression in terms of the azimuthal modulation one wants to 
extract. 

The methods described above are simple and direct, still some systematic effect can be relevant. 
In particular, the apparatus acceptance can introduce correlations between the physical asymmetries. 
For this reason different and more elaborated methods have been developed, which include the binning 
of the data in the {(ph, (ps) plane and the fit with a function which includes all the modulations which 
appear in the cross-section. In the most recent analyses, both HERMES [275(1276] and COMPASS |277[ 
12781 1279] have introduced "unbinned" maximum-likelihood methods based on maximum-likelihood fits 
with the data unbinned in (p^ and (ps- The probability distributions include all the expected azimuthal 
modulations, both for the spin independent and the transverse target spin dependent parts of the 
cross-section. The spin independent part turned out not to infiuence the results for the SSA's. In the 
COMPASS case, the same is true for the acceptance of the apparatus, assumed to have the same relative 
variations for all the target cells |278] . and the SSA's obtained with the "unbinned" maximum-likelihood 
method are in very good agreement with those extracted from the simpler methods used previously. 

In order to obtain the final results, here called SSA's, the raw asymmetries obtained with such fits 
have to be divided by the target (and beam) polarisation. The result are the HERMES "moments" or 
"amplitudes", usually indicated with 2 < /($) >. In the COMPASS data, the raw asymmetries are 
divided by the target polarisation, by its dilution factor /, and by the kinematical y-dependent factors, 
usually dependent on the experimental acceptance in y. As an example, the Collins SSA's published by 
COMPASS are obtained by dividing the raw asymmetries a by the target polarisation, by the dilution 
factor /, and by the mean value of the transverse polarisation transfer from the initial to the final 
quark in the elementary lepton-quark scattering -Datjv = (1 ^ — 2/ + The HERMES results 

for the Collins asymmetry (the so-called "lepto-beam asymmetries") are only divided by the target 
polarisation. 

Concerning the evaluation of possible systematic effects, it has to be noted that, even if the methods 
described above allow to extract simultaneously all the SSA's, the systematic uncertainties have to be 
evaluated independently for all the modulations, since they can be different. For this reason, the results 
are not usually published all at the same time. 

A final remark concerns the measurements of the asymmetries as functions of different kinematical 
variables, typically x, z and Phj_- Since in the experiments there is a strong x — correlation, very 
few attempts have been done to measure the dependence of the asymmetries in the various x bins, 
which would be better studied by comparing the results of the different experiments. Usually the SSA's 
are measured binning the data alternatively in x, z or Ph±, and integrating on the other two variables. 
This extraction introduces some correlation between the data, which should be taken into account when 
fitting all the results in a global analysis. To avoid this problem, the HERMES Collaborations is doing 
multi-dimensional analysis, which are not yet possible in COMPASS, due to the limited statistics. 

5.2 Accessing transversity 

Today, the most direct information on transversity is coming from SIDIS measurements with trans- 
versely polarised targets, which are complementary to the DY experiments and have the advantage of 



allowing a flavor separation by identification of the final state hadrons. 

Among the various SIDIS observables related to transversity, the measurements performed so far 
have provided data on three of them: the Collins asymmetry, the two-hadron asymmetry, and the A 
polarisation. They will be presented in the following subsections, after a a brief description of the event 
and hadron selection. 



5.2.1 Collins asymmetry in SIDIS 

The main source of information on the transversity PDF's is at present the Collins asymmetry, which 
couples hi to the Collins fragmentation function H^. The Collins asymmetry has been measured by 
the HERMES HH [275j and by the COMPASS ^ SSI [2801 M Collaborations. 

Before describing these results, it has to be mentioned that the asymmetries measured by the two 
experiments differ for the already mentioned correction by the D^n factor, applied by COMPASS only, 
and for the sign because of the different definition of the Collins angle $c- In HERMES following the 
so-called "Trento convention" |184j it is defined as = 4'h + 'Psi while in COMPASS the original 
definition [29] = 4>h + (f>s ~ ''^ used, as mentioned in Section I4.1.1[ 

The first signal for a non-zero Collins asymmetry came from HERMES in 2005 [41j, when the results 
on the data collected with the transversely polarised target in 2002 were published. The asymmetry 
had values clearly different from zero in the valence region and of opposite sign for positive and negative 
pions, and this was the first evidence that both the transversity and the Collins FF had to be different 
from zero. An interesting feature of the HERMES results is that the size of the asymmetry turned out 
to be roughly the same for positive and negative pions. As suggested in Ref. [JT], this result implied that 
the favoured {u — )■ 7r+, d — )■ vr^) and the unfavoured {u — ?■ vr^, d — )■ 7r+) Collins fragmentation functions 
iJf*"'^^^ and H^'^^^^ should be of the same size. In fact, neglecting the sea contribution (the asymmetry 
is different from zero only in the valence region) the flavour structure of the Collins asymmetry for a 
proton target can be written as 

^Sl - elh^Ht'^- + elhiHt^--' , A^^'oV ~ elKHt^^-' + e'XHt'^^ . (127) 

Due to the weight factor given by the quark charge, the measured asymmetries are not sensitive to 
/if, thus the result \^qo\\ I — l^coii I imphes that Hi'^'^^ ~ —H^'^^^. This finding, unexpected at the 
time, can be understood jH] in the framework of the string model of fragmentation which in its most 
recent version [177] is described in Section 13.8. 1[ If a favoured pion forms at the string end created 
by the first break, a disfavoured pion from the next break will be opposite in charge and will inherit 
transverse momentum from the first break in opposite direction from that acquired by the first pion. 
Also, assuming that the u and d quark contributions add up in the asymmetries, in the same model it 
is expected that their transversity distributions have opposite sign. 

The COMPASS experiment started data taking with the transversely polarised deuteron target, and 
the first results [12] were published almost at the same time as the HERMES results. The measured 
Collins asymmetries were all compatible with zero, both for positive hadrons and for negative hadrons. 
This result is compatible with the HERMES finding. Limiting again the analysis to the valence region, 
the deuteron asymmetries can be written as 

^Su - {h\ + hi) {elH^''- + elHi'-''') , ~ {h^ + hf) {elH^^'- + e^i/^^^^^) . (128) 

The straightforward conclusion from the COMPASS deuteron measurements is that h^ and hf must 
have roughly the same size and opposite sign, very much as in the case of the helicity quark distributions. 

Both HERMES and COMPASS have continued the measurements with the proton and the deuteron 
targets respectively, and have produced results with considerably better statistics, which have confirmed 
the first measurements. 




Figure 22: HERMES results for the Collins asymmetry from the 2002-2005 data collected with the 
transversely polarised proton target |281] . The asymmetries are shown as function of x, z and Ph±- 
The left plots show the asymmetries for pions, the right plots the asymmetries for charged kaons. 

The HERMES results based on the whole data collected from 2002 to 2005 [2HI] are shown in Fig. 
for pions (left) and charged kaons (right). 

The applied cuts are > 1 GeV^, O.K y < 0.95, > 10 GeV^ and all the hadrons with 2 < Ph < 
15 GeV, 0.2 < z < 0.7, and a polar angle with respect to the direction of the virtual photon larger than 
0.02 rad are used in the extraction of the asymmetries. In the figure the bands represent the maximal 
systematic uncertainty which includes hadron misidentification and acceptance and detector smearing 
effects, and is smaller than the statistical one. The scale uncertainty due to the target polarisation has 
been evaluated to be about 8%. The fraction of charged pions and charged kaons produced in vector 
meson decay has also been estimated. As can seen in the figure, the asymmetries for tt~^ and vr" have 
opposite sign, increase from very small values at x ~ 0.03 to about 5% at the highest x values, and 
have a similar magnitude. The ttq asymmetries are compatible with zero. 

The COMPASS results for the Collins asymmetry from all the data collected from 2002 to 2004 
with the deuteron target are shown in Fig. [23] for charged positive and negative hadrons [13]. The error 
bars are statistical only. The systematic errors have been estimated to be negligible with respect to 
the statistical precision, and the overall scale uncertainty is 7.3% including the uncertainties on the 
target polarisation and on its dilution factors. Here the DIS events are selected requiring > 1 GeV^, 
0.1 < y < 0.9, W'^ > 25 GeV^. The hadrons used in the analysis have Ph± > 0.1 GeV and z > 0.2. 
Following early suggestions by Collins and Artru [282], the asymmetries have been extracted also for 
"leading" hadrons, selected as the highest z hadron with z > 0.25. Also in this case, the asymmetries 
turned out to be compatible with zero. 

Again compatible with zero are the asymmetries measured by COMPASS on deuteron for charged 
pions and for kaons. The final results are shown in Fig. [21] as functions of x, z and Ph± for charged 
pions (top), charged kaons (middle) and neutral kaons (bottom). They have been obtained using all 
the 2002-2004 data for K° and the 2003 and 2004 data for the charged hadrons, since in 2002 the RICH 
was not working during the transverse target polarisation data taking. The kinematical cuts are the 
same as for the unidentified charged hadrons, plus the requirement to have charged pion and kaon 
momenta above 3.1 GeV and 10 GeV respectively, and below 50 GeV. The lower limit is due to the 
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Figure 23: COMPASS results for the Collins asymmetry from the 2002, 2003, and 2004 data collected 
with the transversely polarised deuteron target [13]. The asymmetries are shown as functions of x, z 
and Ph± for all positive (full circles) and all negative hadrons (open circles). In the plots the open 
circles are slightly shifted horizontally with respect to the measured value. 



RICH threshold and the upper correspond to 1.5 cr mass separation between the two mass hypotheses. 

The charged pion and kaon asymmetries have been corrected for the purity of the particle identi- 
fication, which, anyhow, is quite good |280] . The overall systematic errors have been estimated to be 
negligible with respect to the statistical errors. 

The quantitative interpretation of SIDIS data on Collins asymmetries and the extraction of the 
transversity distributions require external information on the other unknown quantity of the process, the 
Collins fragmentation function. This has been recently obtained from the inclusive hadron production 
in e~^e~ annihilation, described in the next section. 

5.2.2 Collins effect in e~^e~ annihilation 

The first indication of the Collins effect in e^e~ annihilation came from a study of the DELPHI data on 
charged hadron production at the pole |283] . which gave an estimate of about 10 % for the analysing 
power (H^)/ (Di), with a considerably uncertainty. 

More recently, data on azimuthal asymmetries in inclusive production of back-to-back hadrons from 
e~^e~ annihilation at s ~ 110 GeV^ have been presented by the Belle Collaboration 1226] . In their 
analysis they use both reconstruction methods described in Section 14.21 The measured quantities are 

„ _ iV(0i + 02) „ _iV(0o) . . 



where A^(0i -|- 02) and A^(0o) are the numbers of hadron pairs with cos(0i -|- 02) and cos 20o modulation, 
respectively, and {N12), (Nq) are the total average number of pairs. In terms of the asymmetries 012 
and Oo defined in eqs. (HSl and [HZl R12 and Rq are given by 

-R12 = 1 + Ctl2 COS(01 02) , -Ro = 1 + «0 COS 200 . (130) 

In order to eliminate the contribution of gluon radiation which is insensitive to the charge of the hadrons 
and the acceptance effects, the ratios of the normalised distributions for unlike-sign (U) hadron pairs 
over like-sign (L) hadron pairs are taken, RY2/R12 Rq /Rq. Focusing on the cos(0i-|-02) modulation, 
one finds ^ 

~ 1 + cos(0i + <P2)A'^^{zi, Z2) , (131) 

^12 
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Figure 24: COMPASS results for the Collins asymmetry for charged pions and kaons, and for neutral 
kaons from all the data collected with the transversely polarised deuteron target |280] . The full and 
open circles refer to positive and negative hadrons respectively. The asymmetries are shown as function 
of X, z and Ph±- 



with the asymmetry parameter given by 



1 + C0S2^ I ^ g2 ^2?^.fev[0]^±,fav[0] ^^±,unf[0]^±,unf[0]j g2/)±,fav [0]^±-..nf [nl ' I ) 



_,unf [1] 
.,unf [ 



where the superscripts "fav" and "unf" denote, as usual, the favored and the unfavored FF's, respec- 
tively. Another independent combination of these functions, given by the ratio R12/R12 of unlike-sign 
pairs over all charged (C) pairs, is also determined, following a suggestion of Ref. [284] . This quantity 
can be written as in eq. (I13ip with an asymmetry parameter A^^. 

An analysis similar to the one we have just sketched is performed with the cos 20o method, leading 
to the ratios Rq /Rq and Rq /Rq, and to the asymmetry parameters A^^ and Aq*"- Notice that the 
cos((/)i + and the cos2</)o analyses of the same events are not independent, and thus cannot be 
included together in a fit. 

The Belle results are presented in Fig. [251 A clear rising behaviour of the asymmetries with zi and 
Z2 is visible, suggesting a similar trend for the ratio H^/Di (recall that the FF's are probed by Belle at 
the scale = s ^ 110 GeV^). Due to the quadratic nature of the asymmetries in terms of the Collins 
function, the difference between H^'^^^ and if^*"'™^ is poorly determined. However, combining the Belle 
constraint on the product H^'^^^ ■ H^'^^^ with the SIDIS measurements of the Collins asymmetry, the 
two Collins functions can be separately determined [i5 l 12851 12841 1286j . 
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Figure 25: Aq (left panel) and A12 (right panel) as a function of Z2 for some zi bins [226]. The UL data 
are represented by triangles and their systematic uncertainty is given by the upper error band. The UC 
data are represented by the squares and their systematic uncertainty is given by the lower error band. 

5.2.3 Phenomenology of the Collins effect and determination of transversity 

The first phenomenological analysis of the SIDIS experimental results on the Collins asymmetry was 
performed by Vogelsang and Yuan |287] . who assumed Soffer saturation of transversity, i.e., \hi\ = 
l/2(/i + gi), to extract from the first HERMES measurement ^Tj the favoured and unfavoured Collins 
functions. Given the poor statistics of those data, the uncertainties on H^'^^^ and H^'^^^ were large. 

Efremov et al. |284[ 1286] analysed the same data with the transversity distributions taken from the 
chiral quark-soliton model |288] . The resulting was shown to reproduce satisfactorily also the 
COMPASS deuteron data |l2l HS], and to be compatible with the Collins function determined from the 
Belle data [44J. The main finding about supports the HERMES interpretation of their data, namely 
that the favoured and unfavoured Collins functions are opposite in sign, and that H^'"^^ is surprisingly 
large, being comparable in magnitude to H^'^'^^ at the average scale (few GeV^) of the HERMES 
experiment. Moreover, it explains why the asymmetry is nearly zero. 

A combined analysis of the first SIDIS data from HERMES and COMPASS, and of the e+e" Belle 
data, was performed by Anselmino et al. flS] and led to the first extraction of the u and rf-quarks 
transversity distributions. This analysis has been updated in Ref. |285] using the preliminary HERMES 
data [275], and the COMPASS [2^] and Belle [226] published data. The fit does not include K"^ 
and 7r° data, nor the preliminary COMPASS results with the proton target. The TMD's are written 
as factorised functions of x and kx, and their transverse-momentum dependence is assumed to have 
a Gaussian form. These two simplifying assumptions are supported by recent lattice studies |289l 
1290] . Fragmentation functions are parametrised in a similar way. Thus the unintegrated unpolarised 
quantities fi{x,k'^) and Di{z,p'^) are expressed as 

fi{x,kl) = fi{x) — TTa^, Di{z,p%) = Di{z) — ■ (133) 
The resulting average transverse momentum of the hadron is 



(134) 




Figure 26: The transversity distributions xAj-m = xh^ and xAj-d = xhf at = 2.4 GeV^ from the fit 
of Ref. |285j . The shaded bands represent the uncertainty of the fit. The sohd hne is the Soffer bound. 
Also shown are the hehcity distributions (dashed curves). 

The widths (k^) and {p^) are those obtained in [291], namely: (fc^) = 0.25 GeV^ and {pD = 0.20 GeV^. 
The parametrisation of the u- and d- transversity distribution and of the Collins function adopted in 
Refs.jlSlEHS] is 

hlix, 4) = AT, - xf imx) + gfix)] 'f^'f , Hi\z,pl) = z\l- zf D,{z)e-^T/^l ^ (135) 

where A^„, A^^^, Nj^av^ ^unf^ f^c are free parameters. Antiquark contributions are ignored. All 

the data used as input of the combined analysis, i.e. the COMPASS deuteron, the HERMES proton 
and the Belle e'^e~ data, are very well fitted. 

Note that since the SIDIS and the e~^e~ data are taken at very different (up to 6 and 20 GeV^ for 
HERMES and COMPASS respectively vs. ~ 10^ GeV^), some assumption about the scale dependence 
of is required. The simple hypothesis adopted in all the phenomenological analyses is that has 
the same evolution as Di, so that the ratio H^/Di is the same at all scales. 

In Fig.[26]we show the transversity distributions extracted from the SIDIS measurements in Ref. |285j . 
They have opposite sign, with \hf\ smaller than \hi\. While the magnitude and the intermediate- a; 
behaviour of hi are reasonably well constrained, its high-x tail is not determined by the data. By 
integration of hi, the tensor charges are found to be 6u = 0.54^q22 and Sd = —0.23~^q^I at the reference 
scale Q2 = 0.8 GeV^ The value for u is smaller than the predictions of lattice QCD |143] and of most 
models. However, one should recall that the model scales are very small and usually just guessed, so 
the evolution from these scales to a higher is affected by large uncertainties |292] . 

A general caveat about the phenomenological analyses of the Collins asymmetry is in order. They 
all ignore the soft factor appearing in the TMD factorisation formulae. At tree level, this factor is equal 
to 1, but as rises it increasingly suppresses the asymmetry |293] . This effect, which is not taken into 
account in the present fits, leads to underestimate the Collins function extracted from Belle data and 
consequently to overestimate the transversity distributions obtained by using that function. 

Given the relatively large values of the COMPASS data in the quark valence region, where the 
HERMES data showed the largest values of the Collins asymmetry, a comparison of the COMPASS 
proton data with the HERMES results was regarded as very important to establish the leading twist 
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Figure 27: COMPASS results for the Collins asymmetry for positive (top) and negative (bottom) 
hadrons from the data collected with the transversely polarised proton target [279] . The error bars 
are statistical only. The curves show the calculation by Anselmino et al. based on the global fit of 
Ref. f2E5] . 



nature of the effect and to check the robustness of the overall picture and of the phenomenological 
analysis. 

In 2008 COMPASS produced the first preliminary results [294] from part of the data collected in 
2007 with the transversely polarised proton target. The results obtained using all the available statistics 
have been published recently [279] and are shown in Fig. |27] for positive (top) and negative (bottom) 
hadrons. The applied kinematical cuts for DIS event and hadron selection are the same as for the 
deuteron data. The systematic uncertainties have been evaluated to be about 0.5 the statistical one, 
and include a 5% scale uncertainty due to the target polarisation measurement. At small x, in the 
previously unmeasured region, the asymmetries are compatible with zero. At larges x a clear signal 
develops both for positive and negative hadron. The results are compatible with the HERMES results 
in the overlap region, and in very good agreement with the values expected values expected on the basis 
the global fit by Anselmino et al. [285] . shown by the curves in fig. |271 For the first time, the Collins 
asymmetry has been measured to be different from zero at ~ 10 GeV, and in perspective these data 
should provide information on the evolution of the transversity and the Collins function. 

5.2.4 Two-hadron asymmetries in SIDIS and e~^e~ annihilation 

The transverse spin asymmetry in the distribution of the azimuthal plane of hadron pairs in the current 
jet of DIS have been measured by HERMES [295] with the proton target and by COMPASS both with 
the deuteron [2961 [2971 [2981 |299l [300] and the proton [278] target. 

Also for this asymmetry there are some differences between the analysis performed by the two 
experiments. The first difference concerns the azimuthal angle of the two-hadron production plane: with 
reference to the definitions introduced in Section 14.1.51 the angle (pR is used in the HERMES analysis, 
while (f)R is used in the most recent COMPASS analysis (see f.i. Ref. [301] ). following the suggestions 
of Ref [282[ 1216] . The two angles, however, coincide in the 'j*N system. Also, HERMES measures 
the amplitude of the modulation in the angle (pR + 05, while COMPASS measures the modulation in 
</'r + 0s + TT, in line with the definition of the Collins angle. Thus the asymmetries measured by the two 
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Figure 28: The two hadron asymmetries versus M/^, x, and z measured by the HERMES experiment with 
the transversely polarised proton target. The bottom panels show the average values of the variables 
that were integrated over. The bands represent the systematic uncertainty. 

experiments are expected to have opposite sign. 

Other differences are in the treatment of the kinematical factor D^n, as in the Collins case, and 
in the kinematical cuts. Finally, the HERMES extraction of the asymmetries is based on a Legendre 
expansion of the dihadron fragmentation functions, as suggested in Ref . |217] . 

In HERMES [295] the events are selected requiring > 1 GeV^, 0.1 < y < 0.85, > 10 GeV^ 
and the missing mass larger than 2 GeV to avoid contributions from exclusive two pion production. 
Also a minimum pion momentum of 1 GeV is required for pion identification. For each event all the 
possible combinations tt+tt^ have been used, labelling as 1 the positive particle. 

If one adopts the partial wave expansion of Ref. |218j , a convenient observable is the asymmetry 



J d cos 9 J d(f)ji J d(f)s sm{(f)R + (f)s) d'^cr / sin 9 
f d cos 9 f d(j)R f dcps d''cr 

which selects the interference fragmentation function Hf'^^. Thus, in HERMES, the spin asymmetries 
defined in eq. f ll26p and divided by the target polarisation, have been measured in each (0/?, 9') bin, 
and fitted with the function asin0Rsin^'/[l + 6(3cos^6'' — 1)]. Here 9' = \\9 — n/2\ — n/2\, a is the 
free parameter and b has been varied to take into account the unknown dependence on cos^ 9 of the 
unpolarised dihadron fragmentation function, where ^^""^fl**"^ = a is the free parameter and b is varied 
to take into account the unknown dependence on cos^ 9 of the unpolarised dihadron fragmentation 
function. The published asymmetries j\^J^r^'^'^^ g^j-e the fitted values of a. The final results from the 
data collected from 2002 to 2005 with the proton target |295j are shown versus M/j, x and z = Zi + Z2. 
The bottom plots give the average values of the other two variables that were integrated over. The 
asymmetries in bins of x and z have been evaluated requiring 0.5 < Mh < 1.0 GeV. The systematical 
errors, which also include the uncertainty due to the value of b and acceptance effects, are given by the 
band centred at -0.01. The scale uncertainty due to the target polarisation uncertainty was about 8%. 
As apparent from the figure, the asymmetries are different from zero, indicating that the spin-dependent 
part of the dihadron fragmentation function is different from zero. Also, the data show clear trends in 
each of the three kinematical variables, and the signal have the same sign and smaller values than that 
of the Collins asymmetry. 
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Figure 29: The COMPASS proton preliminary results for the two hadron asymmetries versus x, z and 

In the COMPASS analysis, the event selection is similar to that described for the single hadron 
asymmetries, namely only events with > 1 GeV^, 0.1 < y < 0.9, W"^ > 25 GeV^ are accepted. 
Only hadrons with zi^2 > 0.1 are used in the analysis and for each pair of hadrons z < 0.9 is required. 
The cuts xp > 0.1 and tt > 7 MeV are also applied in the more recent analysis. In the extraction 
of the asymmetries no dependence on 6 is taken into account. This is justified even in the framework 
of Ref. | 217] . since in the COMPASS kinematics the sin^ distribution is strongly peaked at one (< 
sin^^ >= 0.94) and the cos^ distribution is symmetric around zero. The asymmetries are extracted 
from the azimuthal distributions using the same methods described for the Collins asymmetry. In 
particular the results from the transversely polarised target data were obtained with the "ratio product 
method" . 

Preliminary results using the data collected with the deuteron target have been produced looking at 
different selections for the hadron pair. The asymmetries have been evaluated using all the combinations 
of the positive and negative selected hadrons |296] : taking only the two hadrons with higher transverse 
momentum with different charge combinations |297] : taking only the two hadrons with the highest 
Zi, again for the different charge combinations [297]; taking all the possible combinations of identified 
charged pions and kaons |299j : taking only the charged pions and kaon with higher Zi |300] . All the 
corresponding asymmetries turned out to be compatible with zero and no clear signal could be seen. 
This result could have been expected on the basis of the null result on the measured Collins asymmetry. 

Also for this SSA there was a strong interest for the COMPASS measurement with the proton 
target, since not much variation was expected going from the HERMES to the COMPASS energy. The 
preliminary results from all the data collected in 2007 |278j are shown in Fig. |29] versus x, z, and Mh for 
all the combinations of positive and negative hadrons, selected as in the deuteron case. The systematic 
uncertainties have been evaluated to be not larger than one half the statistical errors. As can be seen, 
the asymmetries are clearly different from zero. They have the same sign of the Collins asymmetries 
on proton, the behaviour in x is very similar and the absolute values at large x are even larger. Also, 
there is no clear indication for a structure, in particular as a function of the invariant mass. 

Very recently, preliminary results on the Artru-Collins asymmetries described in Section 14.21 have 
been produced by BELLE |302j . They show asymmetries different from zero, which depend on Zh and 
Mh- Both these results and the COMPASS proton results are fresh, and no attempt to perform global 
analysis has been done yet. 

On the contrary. The HERMES data on ^^^(<^«+'?^s)sm6» ^^^^^ been recently analysed |230j in terms 
of the model for Di{z,Ml) and H^^p developed in Ref. |219] . It turns out that, in order to get a fair 
description of the HERMES asymmetry, using for the transversity distributions the parametrisation 
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Figure 30: COMPASS results for A and A polarisation versus x from the data collected with the 
transversely polarised proton target |304] . 

of Ref. |285j . the interference fragmentation function H^,.p of Ref. [219] must be reduced by a factor 
3. The predicted dependence of the analysing power Hf^^/Di, with the typical bumps at the u 
and p masses, is not incompatible with the data, within their large errors. However, the interference 
function that fits HERMES data largely undershoot the asymmetry measured by COMPASS and, once 
evolved to high by means of the DGLAP equations for dihadron fragmentation functions [229j . 
yields values of the Artru-CoUins asymmetry which are much smaller than those found by Belle |299j . 
Moreover, the invariant mass behaviour predicted in Ref. |230j does not match the Belle findings. A 
reconsideration of interference fragmentation functions in the light of the recent HERMES, COMPASS 
and Belle measurements seems to be necessary. 



5.2.5 A polarisation 

As shown in Section 14.1.41 detecting a transversely polarised spin 1/2 hadron in the final state of a 
semi-inclusive DIS process with a transversely polarised target probes the leading-twist combination 
hi{x) Hi{z). The typical example of such processes is A (or A) production [32]. The A polarisation is 
in fact easily measured by studying the angular distribution of the A — )• pvr decay. 

From the phenomenological viewpoint, the main problem is that, in order to compute the A po- 
larisation, one needs to know the fragmentation functions H^^'^{z), which are completely unknown. 
Predictions for have been presented by various authors |205l 12061 1207] and span a wide range of 
values. In the calculation of Ref. |205] . where the transversity distributions are assumed to saturate 
the Soffer bound, lies in the interval ±10% at x ~ 0.1. In particular, in the SU(6) non relativistic 
model, the entire spin of the A is carried by the strange quark and one therefore expects a polarisation 
close to zero, due to the smallness of hf. 

The only existing results for this experimentally difficult channel are from the COMPASS exper- 
iment. On deuteron, the preliminary analysis gave polarisation values for A and A compatible with 
zero within the non negligible statistical errors |303j . The same indication comes from the polarisation 
measured with the proton target. Fig. [301 shows A and A polarisation versus x from part of the 2007 
proton data |304] . 

As in the case of the Collins function, independent information on Hi can be obtained from e~^e~ an- 
nihilation [271 1220[ 1228] . The specific process is back-to-back AA inclusive production, e~^e~ — )■ A AX, 
with the hyperon and the anti-hyperon decaying into pn^ and p7r+, respectively. It was shown in 
Ref. [220] that the unpolarised cross section of this process contains an azimuthal modulation propor- 
tional to H^^'^H^^^ = {H^^'^y, which is selected by the asymmetry between the number of pp pairs on 
the same side of the scattering plane and the number of pairs on opposite sides. The measurement of 



such asymmetry was attempted by the ALEPH Collaboration at LEP [3U5], but the scarce sensitivity 
of data did not allow getting any significant result. 

5.3 Accessing TMD distributions: T-odd leading twist functions 

SIDIS represents, at this moment, the best source of knowledge on the two T-odd distribution functions: 
the Sivers function, involved in transversely polarised SIDIS, and the Boer- Mulders function, which 
generates asymmetries in unpolarised SIDIS. Some information on the Boer-Mulders distribution comes 
also from unpolarised Drell-Yan processes. DY can probe the Sivers function as well, if one of the two 
colliding hadrons is transversely polarised, but this reaction so far has not been experimentally explored. 

5.3.1 Sivers effect in SIDIS 

As in the case of the transversity PDF's, the only measurements which give today a clean access to 
the Sivers function are the SSA's in SIDIS on transversely polarised targets. The relevant quantity is 
the so-called Sivers asymmetry Agiv = F^^'^^ '^'^^ /Fjju, where F^^'^'^ and Fuu are the structure 
functions introduced in eq.(l58l). This SSA couples /ft^ to the unpolarised fragmentation function Di. 
Till now it has been measured only by the COMPASS and the HERMES experiments. The same data as 
for the measurement of the Collins asymmetry have been used, and the same kind of analysis, described 
in Section 15. has been performed. In this case the relevant modulation is that in the azimuthal angle 
{(f)h — 'Ps)i and, at variance with the Collins case, the Sivers asymmetry is defined in the same way and 
with the same sign in the COMPASS and in the HERMES experiments. It has to be noted that in the 
HERMES papers the Sivers asymmetry is indicated as 2 < sin(0/j — (ps) >■ 

The first results on the Sivers asymmetry have been produced by the HERMES experiment using the 
data collected in 2002 with the transversely polarised proton target [H] , and showed large positive values 
for the vr"'" and while for tt~ and K~ the asymmetries were compatible with zero. The preliminary 
results from the 2002-2005 data j275j and the final results published recently |276] confirmed with better 
statistical precision the previous measurement. Fig. [3T] shows the final Sivers asymmetries versus x, z 
and Ph± for pions and charged kaons. The error bars are statistical only. The systematic uncertainty 
is shown by the bands and include possible contributions due to the longitudinal component of the 
target spin as well as acceptance and smearing effects, radiative effects, and effects due to the hadron 
identification. The additional scale error due to the target polarisation uncertainty is quoted to be 
7.3%. As can be seen, the tt"*" asymmetry is of the order of 5% almost over all the measured x range, 
compatible with zero for vr", and slightly positive for ir^. Slightly positive signals can be seen also for 
K~ , while the values for are quite large, up to 10%. The pion results can naively be explained in 
the framework of the quark model. With the assumptions used in eq. f ll27p . they can be understood 
as due to a the ci-quark Sivers function of roughly twice the size of that of the u quark, and of opposite 
sign. The K results are more difficult to be explained, and further studies on the difference between 
TT~^ and results are quoted in |276j . Also, the asymmetry in the difference of the distributions of 
vr"*" and vr^, which should be more related to the u and d quark Sivers functions in the valence region, 
has been measured. Finally, a study of the dependence has been performed, without finding clear 
effects indicating sizable 1/Q^ effects. It has to be reminded, however, that the mean values of are 
all in the range between 1 GeV^ and about 7 GeV^. 

As in the Collins case, COMPASS has measured for the first time the Sivers asymmetry on the 
deuteron. First results for charged hadrons from the 2002 data were published in 2005 ^\ , and later on 
final results from all the collected deuteron data have been produced for charged hadrons |13] and for 
identified pions and kaons |280] . The measured Sivers asymmetries for identified hadrons are shown in 
Fig. 1321 The errors are statistical only. The quoted systematic errors are negligible with respect to the 
statistical ones. All the measured values are compatible with zero within the small statistical errors. 




Figure 31: HERMES results for the Sivers asymmetries on proton from the 2002-2005 data [276] as 
functions of x, z and Ph^-, for pions (left) and charged kaons (right). 



This result for pions can be interpreted naively in the framework of the parton model [13] as due to 
opposite Sivers functions for the u and d quarks. 

After the first phenomenological study [108] of preliminary HERMES data on the P/j^-weighted 
Sivers asymmetry [306] . various theoretical groups [29H I287[ 13071 ISUS] extracted the Sivers distributions 
(or their moments) from the HERMES measurement of J^^^^^ '^^^ [H]. The fits were then extended 
[29n 1308] to higher-precision HERMES preliminary data and to the COMPASS deuteron data 
[12] . A comparison of the results of these analyses [310] shows a certain qualitative agreement and 
some common features: a negative j^T and a positive /f!^, as predicted by the the impact-parameter 
approach [141] . with comparable magnitudes, as expected in the large- iVc limit [160] or in chiral models 

USD. 

More recent fits [1101 HI] take into account the new HERMES [276] and COMPASS deuteron data 
[280] . The surprisingly large values of the asymmetry and of the — 7r+ difference call for a 
careful reconsideration of the sea. The Sivers functions are factorised in x and /ct, with a Gaussian 
dependence on fcy Taking as an example the parametrisation of Ref. [110] . the functional form of /{l^ 
is 

/{^/(x, 4) = N,x^(\ - x)^mx) e-'r/^^l , (137) 

where Ng, aq, (3 and are free parameters, the last two being taken to be the same for all flavors. 
It turns out that the exponent /3 governing the high- a: tail of the distributions is not well constrained 
by the data which extend up to x ~ 0.3. The overall quality of the fit is rather good and the first 
moments of the extracted Sivers functions for the u and d quarks displayed in Fig. [331 As expected, a 
non negligible strange sea is required to reproduce the K'^ data. The u and d distributions are more 
uncertain, even in their sign. 

The two fits of Ref. [110] and Ref. [Ill] indicate that the Burkardt sum rule is approximately satu- 
rated by the quark and antiquark distributions, thus little room is left for the gluon Sivers component 
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Figure 32: Final COMPASS results for the Sivers asymmetry on deuteron against x, z and P^^ for 
charged pions and kaons [280] . 

and the orbital motion is restricted to valence quarks. 

In this situation, in which all the existing experimental data on the Sivers asymmetry could be 
explained coherently, the COMPASS measurements with the transversely polarised proton target came 
as a surprise. The preliminary results for charged hadrons from part of the 2007 data were released in 
2008 |294j and showed small asymmetries both for positive and negative hadrons, compatible with zero 
within the statistical errors. The analysis of the complete set of data was concluded only recently |279j , 
and the final results are shown in Fig. [Ml For positive hadrons the data indicate small positive values, 
up to about 3% in the valence region. These values are somewhat smaller than but still compatible 
with the ones measured by HERMES at smaller . The systematic errors is estimated to be ^.%astati 
plus a ±0.01 systematic uncertainty in the absolute scale due to a systematic difference in the mean 
values of the asymmetry which was found between the first and the second parts of the 2007 run. 

Given the importance of the Sivers function the COMPASS Collaboration has decided to remeasure 
SSA's on NH3 in 2010 with an improved spectrometer and better statistics. 

5.3.2 Boer Mulders effect in SIDIS 

The existence of cos 0/^ and cos 20/i asymmetries in unpolarised SIDIS is experimentally well established. 
They have been investigated many years ago by the EMC and ZEUS experiments |31H I312[ 1313] in 
the large region where they are dominated by perturbative QCD effects. It is indeed known that 
these asymmetries are perturbatively generated by gluon radiation |193i I194[ 13141 1195[ I315[ 1316] and 
at high and high P^^ they are dominated by these effects. The recent HERMES [3I7], COMPASS 
|318 ] and CLAS |319] results cover a kinematical region (small and P^^ < 1 GeV) where gluon 
emission is negligible |320] and the asymmetries can be described in terms of TMD's and higher- twist 




Figure 34: COMPASS results for the Sivers asymmetry on proton against x, z and Phi_ for positive 
(dosed points) and negative hadrons (open points) |279] . 



contributions. 

The unpolarised azimuthal asymmetries have been measured by the COMPASS and HERMES 
Collaborations using part of the data collected with polarised targets and combining them in order to 
cancel the net target polarisation. 

The COMPASS preliminary results |318] have been obtained from part of the data collected in 2004 
with the ^LiD target polarised both transversely and longitudinally with respect to the muon beam 
direction. The data have been combining in such a way to cancel the net target polarisation. To reduce 
the acceptance effects, only the events with a vertex in the downstream target cell have been used, and 
the statistics for each polarisation orientation chosen in such a way to have a zero net polarisation. The 
events are selected requiring the usual cuts > 1 GeV^, mass of the final hadronic state W > 5 GeV, 
0.1 < y < 0.9. For the final state hadrons it is required that 0.2 < z < 0.85 and 0.1 < Ph± < 1.5 GeV. 
The data have been binned alternatively in x, z and Ph± and in each bin the measured (ph distribution 
has been corrected for the acceptance of the spectrometer, evaluated with a Monte Carlo simulation. 
Since the COMPASS beam is longitudinally polarised, a sin^^ modulation is also possible, so the 
corrected (ph distributions have been fitted with a function containing the cos(f)h, the cos20/i and the 
sin (ph modulations. This last amplitude turned out to be always compatible with zero. The preliminary 
results for the cos 2(f)h asymmetries are shown in Fig. |35] as functions of x, z and Ph^ for positive and 
negative hadrons. The measured asymmetries have been corrected for the corresponding y dependent 
kinematical factor appearing in the cross-section. The errors are statistical only. The systematic errors 
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Figure 35: The preliminary results for the cos 20 azimuthal asymmetries versus x, z and Phi. for 
positive and negative hadrons measured by the COMPASS experiment on ^LiD [318J. 



are of the order of 2% in both cases and are largely dominated by the acceptance correction. Both 
the COS0/1 and the cos20/j asymmetries are quite large, with a strong dependence on the kinematical 
variables. Also, for the first time, the asymmetries have been produced separately for positive and 
negative hadrons and the measured difference points to different contributions of the u and d quarks to 
the underlying mechanisms. 

The HERMES experiment has produced results for the cos0ft and the cos20/i asymmetries in hydro- 
gen and deuterium using the data collected in 2000, 2005, and 2006. The cuts applied in the event and 
hadron selection are: x > 0.023, 1 < < 20 GeV^, 10 < VV^ < 45 GeV^, 0.3 < y < 0.85, xp > 0.2, 
0.2 < z < 0.75, and 0.05 < Ph± < 0.75 GeV. To correct for acceptance of the spectrometer, detector 
smearing, and QED radiative effects, the data were analysed in a 5-dimensional grid in the variables 
X, y, z, P/j^ and (ph- A 10- dimensional smearing matrix was populated by Monte Carlo simulation and 
incorporated into the fitting procedure, which has been extensively tested with Monte Carlo data. The 
released asymmetries are one dimensional projections of the asymmetries in which the other four vari- 
ables have been integrated over. The measured cos (ph asymmetries for protons and deuterons are shown 
in Fig. [36] for positive and negative hadrons versus x, y, z, and Ph±- The error bars are statistical. 
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Figure 36: The preliminary results for the spin-independent cos (j)h azimuthal asymmetries versus x, y, 
z and Ph± for positive and negative hadrons measured by the HERMES experiment with the 28 GeV 
electron beam |317] on p (top) and on d (bottom). 
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Figure 37: The preliminary results for the cos 20 spin-independent azimuthal asymmetries for deuteron 
from COMPASS (left panel) and HERMES (right) as functions of x, z and compared with the fits 
to the data of ^145j . 

and the bands show the systematic error. As can be seen in the figure, the proton and the deuteron 
asymmetries are very similar, giving a hint for alike sign for the u and d Boer Mulders functions. The 
difference between positive and negative hadrons is remarkable, as in the case of the COMPASS results. 
At variance with the COMPASS data shown here, the HERMES results incorporate the y-dependent 
kinematical factor appearing in the cross section. 

The CLAS results at JLab for the cos20/i asymmetry for tt"*" [319j agree with the HERMES mea- 
surements at large z. A striking feature of their data is the large values measured for z < 0.2, which 
increase with Ph± and for which there is no comparison with other experiments. 

Phenomenologically, the Cahn contribution to the cos (ph asymmetry was studied by Anselmino and 
coworkers [291] . Using the EMC data |31H 1312] in the region Ph± < 1 GeV, they extracted the average 
values of quark momenta: (k^) = 0.25 GeV^ in the distribution functions, {p"^) = 0.20 GeV^ in the 
fragmentation functions. 

The cos2(j)h asymmetry has been analysed in detail in Refs. |320l 11451 1321] . In Ref. |320j . which 
anticipated the measurements, it was predicted that A^^j^'^'' is of the order of at most 5%, and that the 
7r~ asymmetry should be larger than the tt"*" asymmetry, as a consequence of the Boer-Mulders effect. 
These predictions have been substantially confirmed by the experimental results. A fit to the HERMES 
and COMPASS preliminary data has been recently presented in Ref. |145j . It assumes that Ajjjj ^ 
can be described by the leading-twist Boer-Mulders component and by the twist-4 Cahn term (which 
is however only part of the full twist-4 contribution, still unknown). The available data do not allow 
a complete determination of the x and dependence of h^. Thus, the Boer-Mulders functions are 
simply taken to be proportional to the Sivers functions of Ref. |110j . h^'^ = Xqf^, and the parameters 
Xq are obtained from the fit (the Boer-Mulders sea is not constrained by the data and is taken to be 
equal in magnitude to the Sivers sea). The result is 

/i^" = 2.0/fir", hi'' = -l.lf^T\ (138) 

and the comparison with the data is shown in Fig. [37] Since is negative and is positive, the 
u and d Boer-Mulders distributions are both negative. This is what one expects in large-A^^c QCD 
|160] and in some other models |152[ 1157] . The results are also consistent with the predictions of the 
impact-parameter picture |87| combined with lattice calculations |143] . which indicate a u component 
of hj- larger in magnitude than the corresponding Sivers component, and the d distributions with the 
same magnitude and opposite sign. A very recent model calculation of |322] is in good agreement 
with the findings of Ref. |145] . 



5.3.3 Boer-Mulders effect in DY production 



As shown is Section 14.3. II the cross-section for the DY process on unpolarised nucleon is given by eq. |99] 
where the quantities A and v are related by the Lam- Tung relation \ + 2v = 1. 

The NAlO data on vr^A^ DY off a tungsten target show that the angular distribution for the DY 
events do not show any cm. energy dependence nor any nuclear dependence. The data show that the 
value for A is close to 1, as expected by the naive parton model (for massless quarks and no intrinsic 
quark momentum the angular distribution should just be (1 -|- cos^6'). Also, the value for is close to 
expectation: it is essentially consistent with zero, indicating that the annihilating partons contribute 
equally to the transverse momentum of the muon-pair. Both the values of A and of /i are essentially 
independent of any kinematical variable. But the most striking result from this experiment is the large 
value they obtain for v and the strong dependence of v on the dimuon transverse momentum. 

Similar results have been obtained at FermiLab by the experiment E615 which investigated the 
same ti'N ji'^n'X DY process. As shown in Fig. [38] was found to be as large as 30%, and steeply 
rising with Qt-, an effect which is not explained by pQCD and was interpreted as a manifestation of the 
Boer-Mulders mechanism |24U] . Also, by projecting the data points on the axis, the pion valence 
structure function Ft,{xt,) could be precisely determined in the x^^ range from 0.21 to 1, and found in 
good agreement with the extraction of NA3 and NAIO. 
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Figure 38: The amplitude v of the cos 20 modulation in tt A^ Drell-Yan process |323j . The Collins-Soper 
frame is adopted. The curve is the prediction of Ref . [240] . 



Recently, the E866/NuSea Collaboration at FNAL has presented data on the angular distributions 
of dimuon production in pd |324] and pp [70] collisions. Much smaller u values (less than 0.05) than in 
T^N DY are found, as shown in Fig. |39l and the Lam- Tung relation is satisfied. These data could in 
principle give some information about the antiquark Boer- Mulders distributions |325[ 1326] . However, 
in the Pt region above 1-1.5 GeV they are expected to be described by pQCD |238] . Thus the only 
points that have likely to do with the Boer-Mulders effect are those below Pt ~ 1.5 GeV. 

We conclude this section by mentioning that a theoretical study of the DY azimuthal asymmetries 
at small and moderate Pt in the context of twist-three factorisation has been performed in Ref. |251] . 
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Figure 39: Left: the v coefficient in the Colhns-Soper frame for three DY measurements (dots: 
E866/NuSea (324j : squares: E615 [69]; triangles: NAlO [323]). The curves are fits similar to that 
of Ref. [240] ■ Right: the v coefficient for pd [324] and pd DY [70]. The dot-dashed curve is the pertur- 
bative QCD contribution. The sohd and dotted curves are the calculations of Ref. [325] for pp and pd, 
respectively, based on the Boer- Mulders effect. 

5.4 Accessing the TMD distributions: leading-twist T-even functions and 
higher- twist functions 

At leading twist, besides the Collins, Sivers and Boer-Mulders terms, there are other three angular 
modulations in the SIDIS cross section, which probe T-even TMD distributions: 

- Unpolarised beam and longitudinally polarised target: sin 20/1 modulation, involving h^j^. 

- Unpolarised beam and transversely polarised target: sin(30/i — (ps) modulation, involving h^rp. 

- Longitudinally polarised beam and transversely polarised target: cos(0/j — (ps) modulation, in- 
volving giT- 

From the data collected with the deuteron target, the COMPASS experiment has produced prelim- 
inary results on these asymmetries for charged hadrons [32 7[ 1328] . All in all, it is hard to find a signal 
in any of these observables. It will be interesting to look at the corresponding results with the proton 
target. 

The HERMES collaboration has measured the sin20/i moment, both with a proton [329] and a 
deuteron target [330] , finding it to be compatible with zero. Signals of a non vanishing sin 2(j)h asymmetry 
have been recently reported by the CLAS collaboration [331] . 

For the moment, all these data can only be confronted with model predictions. Focusing on h^rp, 
which has attracted some attention for its interesting physical content, a calculation of A!^™^^'^^ '^^^ 
based on positivity bounds [121] gives limits of about ±0.03 for the asymmetry on a deuteron target 
and a slightly larger value for a proton target. The COMPASS deuteron data lie within these bounds. 

The light-cone constituent quark model of Ref. [165] predicts an asymmetry ^^^'1'^''^'^^ smaller than 
0.01, also consistent with the COMPASS findings. All the asymmetries related to T-even TMD's are 
calculated in Ref. [165] and found to be generally close to zero, hence compatible with the COMPASS 
findings. 



At subleading twist, that is at order 1/Q, the situation is much more involved. The SIDIS structure 
functions of eq. f l58|) have have in fact the general form: 

F ~ l.t. TMD ® h.t. FF + h.t. TMD ® l.t. FF 

where "l.t." = leading twist, and "h.t" = higher twist. They contain both leading-twist and twist-three 
TMD distributions and fragmentation functions. Thus, the phenomenological interpretation of these 
observables is rather intricate. In SIDIS with unpolarised (U) and/or longitudinally (L) and transversely 
(T) polarised beams and targets there are 8 twist-three modulations: 

UU : coscph] LU : sincph; UL : siiKph, sin 2(1) LL : cos^/^; UT : sin 05-; 

UT : sm{2(f)h - (ps); LT : cos 0s, cos(20/j - 0s) . 

COMPASS has measured all these quantities on a deuteron target [32 7[ I318[ 1328] and found them all 
to be consistent with zero. HERMES presented results on the sin20/i and siiKph modulations with 
a longitudinally polarised proton |329j and deuteron target [330]. While the sin20/i asymmetry was 
found to vanish, the sin0/i asymmetry showed a large signal (up to 4% for proton and 2% in deuteron), 
incompatible with zero for positive and neutral pions, originally interpreted in terms of large transversity 
PDF's. 

A remark about the definition of the target polarisation is now in order. Experimentally, the target 
polarisation is defined with respect to the a longitudinal (transverse) polarisation with respect to the 
beam axis has a transverse (longitudinal) component with respect to the virtual photon axis, which 
is kinematically suppressed by a factor 1/Q |112j . Thus, any measured w{(f)h,4>s) modulation with a 
"longitudinally" ( "transversely" ) polarised target is mixed with a transverse (longitudinal) modulation 
of the same type, suppressed by 1/Q. This may be relevant in some cases. For instance, the sin0ft 
asymmetry measured with a target longitudinally polarised with respect to the beam axis gets a twist- 
three contribution from A^^Jj''', but receives also contributions from the leading- twist Collins and Sivers 
asymmetries, A^^'^'''^'^^'^ and A^,^'^'^~'^^\ which are multiplied by a kinematical factor ~ In prin- 
ciple, all these contributions might be equally relevant. The HERMES analysis of the sin0/i asymmetry 
|332] has shown that the Collins and Sivers contributions to (sin0/i) are small, and therefore the large 
measured asymmetry is a genuine subleading- twist effect. 

Another indication of the relevance of twist-three effects comes from the beam-spin asymmetry 
^Ti/'^y which has been measured by CLAS at 4.3 GeV for positive pions |333] and by HERMES at 27.6 
GeV for charged and neutral pions |334j . The asymmetry for positive pions is large and positive. The 
CLAS and HERMES results nicely agree if one rescales the HERMES data by the mean value of 
and the y dependent kinematical factor, as shown in Fig. HOI 

5.5 Inclusive hadroproduction 

Sizable SSA's in polarised inclusive hadroproduction have been reported since the 70s [H [21 [3l 1335] . 
for center-of mass energies in the range 5-10 GeV. In the same years, Fermilab experiments discovered 
that A hyperons produced in unpolarised pp collisions have a large transverse polarisation with respect 
to the production plane |H [5]. These findings provoked a certain theoretical interest, as it was widely 
held that large transverse polarisation effects could not be reproduced in the framework of perturbative 
QCD [16]. On the other hand, the small values of ^/s and Pt explored by those experiments left the 
door open to interpretations based on soft physics. 

In 1991 the E581/E704 experiment at Fermilab extended the investigation of transversely polarised 
hadroproduction to higher energies and found remarkably large transverse SSA's in the forward region 
[HI E]- These results were confirmed by the RHIC measurements, which moved the energy frontier 
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Figure 40: Comparison of the kinematically rescaled sin 0/j asymmetries between the HERMES (circles) 
and CLAS (triangles) measurements. The full square represents a previous HERMES measurement 
averaged over the range 0.2 < z < 0.7. The outer error bars represent the quadratic sum of the 
systematic uncertainty and the statistical uncertainty (inner error bars) |334] . 



one order of magnitude upwards. On the A polarisation front, the striking effect discovered in the 
early experiments was observed also at higher Pt values |336l 13371 1338] . Thus, it has become clear 
that there must be some hard mechanism behind the transverse polarisation phenomena observed in 
hadropro duct ion . 

5.5.1 SSA's in inclusive hadroproduction 

As seen in Section 15. the left-right asymmetries associated to an azimuthal modulation of a cross- 
section, are best measured by comparing data taken with up- and down-polarised beam or target. 
Integrating over 0, a global left-right asymmetry Aj^i can then be defined as 

1 Nl-Nr 

Am — , 

with V the beam or target polarisation. 

The modern era of experimental (and theoretical) work on SSA's in hadroproduction was inaugurated 
by the E704 investigation of pN and pN collisions with transversely polarised secondary proton and 
antiproton beams at the cm. energy ^/s ~ 19.4 GeV [El [9]. Two kinematical regions were covered: 1) 
the beam fragmentation (or forward) region, 0.2 < xp < 0.6, with in the range 0.2 — 2.0 GeV; 2) the 
central rapidity region, [3:^1 < 0.15, with Pt up to 4 GeV. In pion production, large SSA's were found 
at high a; [H [7l IHl [10]: the asymmetries are nearly zero up to ~ 0.3 and then start rising with xp, 
reaching 15 % for 7r° and 30-40 % for tt^. In p'^N collisions is positive for 7r+ and negative for tt~ , 
with about the same size. Signs are reversed in ^'''A^ scattering. The asymmetry for is roughly half of 
that of charged pions and always positive. As for the Pt dependence, A^r is zero below Pt ~ 0.5 GeV, 
and above this value increases in magnitude with Pt- In the central rapidity region, where larger values 
of Pt are reached, the asymmetries turn out to be consistent with zero |339] . E704 has also measured 
sizable SSA's in inclusive A and t] production [340] 1341] but the rather low transverse momentum of 
the A's does not allow a safe perturbative QCD analysis. 



(139) 



The E704 findings have been substantially confirmed by other fixed-target experiments at lower 
energies, at IHEP (Protvino) [3l2l |3l3] and at the BNL-AGS [3MIM5] . 

On the phenomenological side, the E704 pion asymmetries have been interpreted in terms of the Sivers 
effect [3461 [391 EST] and of the Collins effect [HEl [40l [347] . A recent reassessment of the situation [348] 
has shown that, contrary to a previous prediction of a strong suppression of the Collins effect [349] . both 
Collins and Sivers mechanisms may give sizable contributions to the SSA's. The E704 results have also 
been studied in the context of twist-3 factorisation, considering quark-gluon correlations in the initial 
state [531 (SSI 12631 1264[ 1350] and in the final state [265] . All these approaches are able to reproduce at 
least qualitatively the data, thus showing that many different physical mechanisms may be at work in 
polarised hadroproduction. Thus, it is impossible for the moment to draw definite conclusions as to the 
dynamical source of single-spin transverse asymmetries. 

The E704 measurement might have left the doubt that transverse SSA's would disappear at collider 
energies. Studying the p^p — )■ 7r° X reaction at a/s = 200 GeV in the first polarised collisions at RHIC, 
the STAR Collaboration showed that this is not the case: the large effects found by E704 persist at 
an order of magnitude higher energy [TT|. As shown in Fig. |?T] STAR measured a large positive 
above xp ^ 0.3 in the transverse-momentum range 1.0 < {Pt) < 2.4 GeV. More recently, the negative 
xp region has been explored, finding an asymmetry consistent with zero, and the Pt dependence of 
the SSA has been determined [351] . In Fig. [42] one sees that the rise of the SSA's at large xp is fairly 
well reproduced by the Sivers mechanism in the generalised parton model with the Sivers function 
extracted from the HERMES SIDIS data [2571 1352] and by the twist-3 factorisation scheme [350] . On 
the contrary the Pt behaviour of the data, showing a clear tendency to increase at fixed xp, contradicts 
the theoretical expectations, which predict a decrease of with Pp. 

Concerning the description of hadroproduction SSA's in terms of TMD's taken from SIDIS analysis, 
one should recall that RHIC asymmetries scan the parton distributions over a wide range of the Bjorken 
variable, including the large-x region, whereas the SIDIS data are limited to a; < 0.3 and do not constrain 
the tails of the transversity distribution and of the Sivers function. Thus, the generalised parton model 
predictions of the SSA's are quite uncertain at high xp [353] . 




Figure 41: The asymmetry Ai\i{p^p — )■ n^X) measured by STAR [TT] at ^/s = 200 GeV. The curves are 
the predictions of Ref. [10] (solid), Ref. [32] (dotted), Ref. [SB] (dot-dashed), Ref. [2SS| (dashed). 

Measurements of forward charged pions by BRAHMS at ^/s = 200 GeV [354j and ^/s = 62.4 GeV 
[355] shown in Fig. US] confirm the asymmetry pattern observed by E704, with the mirror effect of vr"'" 
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Figure 42: The asymmetry A^^p^p — > vr^X) measured by STAR [351] at a/s = 200 GeV, as a function 
of xf for two different values of the average pseudorapidity (rj) (left) and the same asymmetry as a 
function of Pt at fixed xp (right). The curves are the predictions of Ref. [25 7[ 1352] (dashed line) and 
of Ref. [350] (sohd line). 

and TT^ and large absolute values of A^. The SSA's for tt^ and tt^ at ^/s = 62.4 GeV are plotted in 
bins of Pt in Fig. HH A clear rise with Pt is visible up to a transverse momentum of about 1 GeV, 
where the asymmetries reach magnitudes of about 0.3. BRAHMS has also measured kaon production 
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Figure 43: Charged pion asymmetries measured by BRAHMS at ^/s = 200 GeV ]354] and y/s = 62.4 
GeV [SSS], and by E704 at = 19.4 GeV [8]. 

at ^/s = 62.4 finding that the asymmetry has the same sign and approximately the same magnitude 
as the K~ asymmetry [355j . The results are fairly well reproduced both by the twist-3 approach and 
by the Sivers effect in the generalised parton model 

In the midrapidity region (|?7| < 0.35). the PHENIX measurements of neutral pion and charged 
hadron production at y/s = 200 GeV |ll] at y/s = 200 GeV, showed transverse SSA's consistent with 
zero. This result is in agreement with the fixed-target E704 finding and extends it to higher Pt, up to 
5 GeV. 
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Figure 44: Transverse SSA of charged pions at ^/s = 62.4 GeV in bins of Pt, measured by BRAHMS 
[355] . 

5.5.2 Spin-averaged hadroproduction cross sections 

The QCD-based descriptions of the E704 and RHIC resuhs on hadroproduction SSA's have been criti- 
cised |356] on the ground that, whereas the spin-averaged cross sections at a/s = 200 GeV [T2[ 13571 fT5] 
are well described by next-to- leading order perturbative QCD |358] , the cross sections at lower energies 
(for instance at y/s = 19.4 GeV (339] ) are not. The conclusion of Ref. [356] is that the single-spin 
asymmetries discovered by E704 and those found at RHIC are different physical phenomena. However, 
it has been recently shown that the resummation of large logarithms arising from soft gluon radiation 
in the limit xt = IPnpj^ — )■ 1 significantly improves the agreement of perturbative QCD calcula- 
tions (limited however to the cross sections integrated over rapidity) with data at low and intermediate 
energies |359[ 1360] . Predictions for GSI-FAIR and J- PARC kinematics are given in Ref. |361] . The in- 
clusion of intrinsic transverse momentum effects has been also shown to reduce the discrepancy between 
QCD-parton model predictions and cross section data in the fixed-target regime [257]. 

Thus, the hope to achieve a consistent QCD picture of both polarised and spin-averaged hadropro- 
duction phenomena seems to be well founded. 

6 Future measurements 

In this section account is given of near future and more distant future SIDIS and DY experiments. The 
main goal of the SIDIS experiments is to measure the SSA's over a detailed grid of the kinematic variables 
X, P/i_L, and 2, facilitating the extraction of the DF's and of the FF's. In addition, it should be possible 
to study sub- leading twist effects by probing their l/Q dependence, and to explore the transition from 
non-perturbative small transverse momentum, typically lower than 1 GeV, to the transverse momentum 
large regime. In the case of the DY measurements, the main goal is to perform for the first time 
experiments with polarised nucleon targets and/or polarised beams, to test the predicted test of sign 
of the T-odd DF's. 

6.1 SIDIS experiments 

COMPASS at CERN: The analysis of the 2007 transversely polarised proton data is still ongoing, 
and several results on SSA's have still to be obtained. Many more data on the same target (NH3) will 
be collected in the long 2010 run, so that in the near future a large amount of data is expected. 
For what concerns a more distant future, the COMPASS Collaboration is presently preparing a pro- 
posal for measurements aiming to study chiral perturbation theory, generalised parton distributions via 
Deeply Virtual Compton Scattering (DVCS), and TMD parton distributions via Drell-Yan processes. 



The DVCS measurements will be performed using a liquid hydrogen target and a 190 GeV muon beam. 
In parallel with these measurements, SIDIS data will be collected to extract with high precision the 
unpolarised cos(j)h and cos20/i azimuthal asymmetries as well as the beam helicity dependent sin0/i 
asymmetry. Such information cannot be extracted from the data collected with the transversely po- 
larised target because of the complications of using a nuclear target. 

JLab experiments: In the near future also the CLAS Collaboration in Hall B will take SIDIS data 
on a transversely polarised target. For transverse running, the use of a novel HD-ice target is planned, 
which in a frozen-spin state requires only small holding fields. The use of the HD-ice target by the 
E08-015 Collaboration has many advantages: being a solid target, it can be short, a few cm, and thanks 
to the smallness of the holding field it can be located in the centre of the detector, thus increasing the 
acceptance of the spectrometer. In addition, the HD target has almost no dilution, which maximises 
the figure-of-merit, and being of low atomic number, comparatively few bremsstrahlung photons will be 
produced in the target. The experiment should run in the second half of 2011, and an upgraded version 
of the detector has already been proposed for JLabl2. 

In a medium term range, Jlabl2 GeV upgrade could meet the requirements to study TMD's in the va- 
lence region, thus covering a complementary kinematic region with respect to COMPASS. The Clasl2 
experiment in Hall-B is designed to achieve a very broad kinematic coverage while increasing by a factor 
10 the luminosity with respect to the current 6 GeV setup. In particular, the forward spectrometer 
comprises a 2 T toroid with improved geometry to minimise the not-active azimuthal coverage and 
a RICH detector is under study to extend the hadron identification over the full energy range of the 
experiment. The spectrometer is complemented by a central detector embedded in a 5 T solenoid. 
Also, an upgraded version of experiment E06-010 (PR09-018) has already been proposed and condi- 
tionally approved to run in Hall A. The experiment aims to measure the SSA's of the SIDIS process 
e'^n — > e'hX, where h is either a tt or a The experiment will use the large-solid-angle Super BigBite 
Spectrometer as hadron arm, the BigBite Spectrometer as electron arm, and a novel polarised ^He 
target that includes alkali-hybrid optical pumping and convection flow to achieve very high luminosity. 
Thanks to the large acceptances of the electron and hadron arms, an electron - polarised nucleon lumi- 
nosity at the level of 4 x 10^^ cm~^s~^, and a target polarisation of 65%, the experiment should collect 
in a two-month run about 100 times more statistics than that obtained by the past experiments. 
e — N and e — A future colliders: In depth studies of hadron structure can be best performed at 
a high energy polarised electron-polarised proton collider. Large Collaborations at BNL and JLab are 
elaborating proposals which are well advanced and are being encouraged by the USA agencies. As 
written in the NSAC 2007 Long Range Plan, "the allocation of resources arc recommended to develop 
accelerator and detector technology necessary to lay the foundation for a polarised Electron Ion Collider 
(EIC). The EIC would explore the new QCD frontier of strong color fields in nuclei and precisely image 
the gluons in the proton" . To carry out a rich and diversified physics program the recommended energies 
for the electrons are between 3 and 10 GeV, for the protons between 25 and 250 GeV, and for the heavy 
ions between 25 and 100 GeV. The luminosity in the case of the e — p collider should be 10'^'^ — 10'^'^ 
cm~^s~^, i.e. about 100 times the luminosity of the HERA collider. Recently preliminary ideas for a 
polarised electron-nucleon coUider (ENC) at GSI, Darmstadt, have been discussed mostly amongst the 
German community. 

The advantage of the Collider configuration over fixed target experiments are manifold: 

- it provides a large range of Q"^, x, W and Phj_- 

- the figure of merit for asymmetry measurements is very much better. For ammonia (NH3) / ~ 0.15, 
thus the figure of merit when scattering on a pure proton beam is better by a factor ~ 50. Needless to 
say, the comparison is done assuming the same number of collected events, so a high luminosity for the 
EIC is a prerequisite. 

- It provides access to the interaction region, so that modern vertex tracking systems can identify short 
living particles, like produced in the interaction. 



Having access to the interaction region, exclusive reactions are at reach. This opens up the whole field 
of GPDs, which to-day are the only way to quantify how the orbital motion of quarks in the nucleon 
contributes to the nucleon spin. Also, it allows measurements in the target fragmentation region, 
which is presently poorly known due to the difficulty of measuring slowly moving hadrons in fixed- 
target experiments, opening a window to the study of spin- independent and spin-dependent "fracture" 
functions. 

The US groups of RHIC and JLab are proceeding jointly to the formulation of two different proposals 
for two different colliders, eRHIC and ELIC, based in the two different laboratories. The RHIC project 
clearly foresees the use of the highly polarised proton and nuclear existing beams. Two accelerator 
design options are being worked upon, both aiming at high brightness 10 GeV electron beams. A 
Ring-Ring option, which requires a new electron storage ring for polarised electron or positron beams, 
is technologically more mature, and could provide a peak luminosity of 0.5 x 10^^ cm~^s~^. The second 
option is a Linac-Ring option, which offers higher luminosity (by a factor of 5) and possibly higher 
energy, but requires intensive R&D for the high-current polarised electron source. 
The starting point of the JLab project is the availability of the 12 GeV electron beam from the upgraded 
CEBAF. The proton complex has to be built from scratch, so it is being designed taking full advantage 
of the expertise matured at RHIC and other laboratories on acceleration and storage of polarised proton 
beams. The design goal for the collider luminosity is very ambitious, 3 x 10^'^ cm~^s~^ for beam energies 
of 10 and 250 GeV for the electron and protons respectively. 



Quite recently, in the summer of 2008, discussions starteqj, about a possible low-cost realisation of an 
ENC at GSI. The central idea is to use the 15 GeV high energy storage ring HESR, which is planned 
to store an antiproton beam for the PANDA experiment (and possibly PAX) as the ring where to store 
the polarised proton beam. By constructing a 3-3.5 GeV electron ring, a "low energy" ENC could be 
realised. The cm energy would be 14 GeV, i.e. in between the HERMES and the COMPASS energies. 
To inject polarised protons in HESR a new 70 MeV p-linac will be needed. The protons would then be 
injected into the existing SIS18 ring, accelerated up to 1.4 GeV, and transferred then into HESR. New 
hardware for the spin manipulations will be needed in SIS 18 and in HESR, but it is the same which will 
be necessary for the PAX experiment. The electron complex has to be constructed from scratch. An 
e-linac and an electron synchrotron will accelerate the electron beam, which will be stored into a new 
storage ring of about the same length as HESR, and housed in the same tunnel. Preliminary machine 
studies indicate that a luminosity of at least 10'^^ could be achieved, as well as large polarisations ( 80%) 
for the two beams. To further reduce the cost of the project, it is proposed to use the PANDA detector, 
and to operate the collider in time sharing with the PANDA Collaboration. 



As described in section 14.31 the DY process in transversely polarised hadron scattering is theoretically 
a very clean and safe way to access transversity. The original suggestion of measuring DY in p^p^ 
scattering, which could be done at RHIC, turned out to be difficult because of the small value expected 
for the asymmetry, of the order of 1 — 2% [5 ^1245(1246] . The measurement will also require external input 
to disentangle the quark and the antiquark distributions. It will be done when the RHIC luminosity 
will be increased. 

These problems can be circumvented by studying DY production with polarised antiprotons at 
moderate energies, which is the ideal process to observe a sizable double transverse asymmetry [83l I362[ 
13631 1247] , dominated by the valence distributions. Such a measurement has been proposed by the PAX 
Collaboration |232j at the FAIR complex to be built at GSI. Since the production rate for relatively 
large dilepton masses M (> 4 GeV) might be too small to allow an easy measurement , it has been 

^Private communication from D. von Harrach. 




6.2 Drell-Yan 




Figure 45: Left: the prediction of Ref. |364] for the Drell-Yan Sivers asymmetry to be measured by 
COMPASS in np colhsions. Right: the double spin asymmetry integrated over in transversely 
polarised proton-antiproton DY for different values of ^/s |247j . 

proposed [362\ to exploit the J/tp peak to measure the asymmetry. 

Finally, we recall that the Sivers effect can also be observed in Drell-Yan processes with a transversely 
polarised proton, where it gives rise to a sin(0^ — (ps) asymmetry. No measurement has been made so 
far, but many experimental collaborations worldwide plan to investigate this class of reactions in the 
near future. 

6.2.1 The proposed experiments 

In the following we describe very briefly the DY proposed experiments. 

COMPASS: Among the proposed measurements for a second phase of the COMPASS experiment 
an important issue is the possibility to investigate for the first time a vr" induced DY process on a 
transversely polarised proton target. The high mass of the COMPASS target (about 1 kg of NH3) 
and the excellent performance of the COMPASS spectrometer make this measurement feasible, and the 
number of events collected in two years of running would allow to check the expected change of sign 
of the Sivers function. Assuming for the magnitude of the Sivers function the value extracted from the 
HERMES measurements in SIDIS (see Fig. HSl left), the significance of the measurement is expected 
to be 3 to 4 cr. 

PAX and PANDA at GSI: As already mentioned, the PAX Collaboration has proposed to measure 
DY processes in pp scattering at FAIR. An asymmetric collider is proposed, consisting of HESR, where 
polarised protons will be stored, and of a new storage ring for the polarised antiprotons, which could be 
the existing COSY Storage Ring, suitably modified. The predictions for Af^^ in p^p^ collisions at GSI- 
FAIR are shown in Fig. HS] (right). One sees that asymmetries of the order of 10-15 % can be expected 
at PAX. Polarisation of the stored antiproton beam will be done using the "spin filtering" technique 
|365j . The antiprotons beam traverses a polarised proton storage cell, and a beam polarisation builds 
up by repeatedly passing through the cell as long as the cross section for parallel spin is different from 
that of antiparallel spins. The method has been proven to work for proton beams, but in the case 
of antiprotons the spin dependent cross-sections are not known and corresponding measurements have 
been proposed at the Antiproton Decelerator at CERN. 

The PANDA collaboration also envisage the measurement of DY process in pp scattering, where both 
particles are unpolarised. The option to put a transversely polarised proton target in the PANDA 
detector is very interesting, but it is technically very difficult and has presently been discarded. 



Fermilab Experiment E906: This experiment is scheduled to run in 2010 for 2 years of data col- 
lection. It will extend DY measurements of E866 (which were done with 800 GeV protons) using an 
upgraded spectrometer and the 120 GeV proton beam from the main injector. The use of the lower 
beam energy gives a factor 50 improvement of luminosity with respect to E866. To cut down costs, it 
will use many components from E866, and data will be taken with Hydrogen, Deuterium and Nuclear 
Targets. The main goal of the experiment is the study of the structure of the nucleon, in particular the 
d/u ratio at high x. 

In the future, there are plans to measure SSA's on a transversely polarised target, and check the change 
of sign of the Sivers function with respect to SIDIS. 

J— PARC: Two proposals for DY experiments have been submitted: 

- P04: measurement of high-mass dimuon production at the 50 GeV proton synchrotron; 

- P24: polarised proton acceleration. 

The advantage of an experiment at JAPRC is the high proton beam intensity, and consequently the 
high luminosity. The disadvantage is that at the same invariant mass the cross section is smaller at 
lower energy. The transverse polarisation program is best carried on by P24, but is clearly conditioned 
by the realisation of the polarised proton beam, which is not yet approved. 

STAR and PHENIX at RHIC: According to the present accelerator schedule, which foresees a long 
longitudinal run for W-physics, a DY program with transverse spin at PHENIX and STAR will not 
start before 2015. However, since in the intersection regions IP-2 and IP-10, where the PHOBOS and 
the BRAHMS experiments were installed, there are no spin rotators, the polarisation of the beams is 
always transverse and ideas are being put forward to prepare both a collider experiment and a fixed 
target experiment for DY measurements. 

RHIC internal target: Quite recently ideas have been put forward for a DY experiment at RHIC 
scattering one beam off an internal target. With the 250 GeV beam, the kinematic range explored 
would be xi = 0.25 — 0.4 (x2 = 0.1 — 0.2), which would be ideal to investigate the change of sign of the 
Sivers function. To achieve the necessary luminosity, the use of a pellet target is being investigated, in 
two different scenarios according to the target thickness (parasitic running or dedicated experiment). 

NIC A: To investigate the hadron structure a Nuclotron-based Ion Collider fAcility (NICA) is being 
planned at the JINR in Dubna, based on the existing proton synchrotron Nuclotron. The accelerator 
complex will require new ion and polarised proton sources, a new linear accelerator, a new booster 
synchrotron and the two new superconducting storage rings of the collider. Both polarised proton and 
polarised deuteron beams should be available in the two rings. The main physics objectives will be the 
study of elastic processes and of Drell-Yan processes. 



6.2.2 Summary 

The relevant parameters for all these projects are summarised in Table [TJ 

Theoretical predictions for COMPASS (vr/), PAX (pp^), RHIC (p^ at = 200 GeV) and J- 
PARC (p^p at ~ 10 GeV) have been presented by various authors [TOHl IM |367l EMI ESHl ETQl l364l 
13711 1372] . The conclusion one can draw from these analyses is that the future experiments will largely 
be complementary to each other. 

As for p^p DY, the RHIC data in the negative xp region will probe the contribution of the sea Sivers 
function, while experiments at lower energies, like J-PARC (operating at ^/s ~ 10 GeV) will provide 
information on the large- a; behaviour of f^. A comprehensive discussion of all future DY measurements 
of Sivers asymmetries can be found in Ref . |364] . 
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Table 1: Compilation of the relevant parameters for the future planned DY experiments. For RHIC, 
IT stays for Internal Target 

7 Conclusions and perspectives 

The original finding of the EMC collaboration, that the quark spin does not account for the total spin of 
the proton, has been a strong motivation for in-depth studies of the QCD structure of the nucleon and for 
a new generation of experimental investigation of hard scattering processes on polarised nucleons. The 
growing interest in the contribution of the quark and gluon orbital angular momentum to the nucleon 
spin naturally led to an increased attention to transverse spin and transverse momentum phenomena. 

In this context, the most important experimental finding has been the discovery that there is a 
correlation between the spin of a transversely polarised quarks and the Pt of the hadrons created in 
the quark hadronisation process. Convincing evidence for this correlation has been provided by both 
SIDIS processes on transversely polarised nucleons and high energy e'^e~ annihilations into hadrons. 
Thanks to this correlation, it is now possible to measure the transversity distribution function. Global 
analysis of the existing SIDIS data and of the e^e~ data have already provided first rough information 
of this two new function. 

A second important discovery of the recent years is that there is also a non-zero correlation between 
the spin of a transversely polarised nucleon and the intrinsic transverse momentum of the quarks. 

In polarised proton-proton scattering the most impressive result in transverse spin physics is the 
confirmation that large SSA's for inclusive mesons production persist at centre of mass energies which by 
now arc more than one order of magnitude greater than those of the previous fixed target experiments. 

On the theoretical side, the main achievement has been the discovery that the Wilson line structure 
of parton distributions, which is necessary to enforce gauge invariance, has also striking observable 
consequences, allowing for single-spin asymmetries that would otherwise be forbidden by time-rcvcrsal 
symmetry. At leading order the most general descriptions of SIDIS and Drell-Yan processes have been 
revisited and a number of structure functions have been introduced to take into account all possible 
correlations among the transverse momentum and spin of the quarks and the spin of the nucleon. Many 
QCD studies have been performed to understand the properties and the gauge structure of these unin- 
tegrated distribution functions, which have been named transverse momentum dependent distribution 
functions. Non-collinear factorisation schemes have been developed and extended to polarised pro- 
cesses. In an alternative approach, twist-three effects have been evaluated and compared to the TMD 
description in the intermediate Pt region. A third line of attack which is vigorously being pursued is 
QCD computation on the lattice. Recent refined lattice QCD results shed fight on the fine structure 
of TMD's. This approach is particularly interesting for transversity, because the tensor charge, the 
first moment of the transversity distribution, which is an all-valence object, is believed to have been 
evaluated with good accuracy on the lattice, thus a good measurement of this quantity could provide a 
good test of the correctness of the calculation. 



In spite of these achievements, the amount of work which is still needed is not small. 

On the experimental side, the available SIDIS data are only a glimpse to a new territory, and many 
more data are needed to obtain the Pt and Q'^ dependence of the asymmetries in the different x-bins, a 
prerequisite to a model independent extraction of the TMD functions. The situation is worse however 
in the DY sector, since no polarised DY data exist at all. The existing unpolariscd DY data allow only 
to access the Boer-Mulders function, but the full exploitation of the DY potential is not even at the 
horizon. 

On the theoretical side, the present fits to the data do not make use of the evolution schemes 
which are already now available for the TMD's and for the quark-gluon correlators, and which clearly 
must be integrated in the calculations. Also , a better understanding of higher-twist contributions to 
SIDIS observables (which are significant since (Q^) is rather small) requires more insight and more 
effort. 

In the near future more data will be collected in the SIDIS sector by the COMPASS and by the 
JLab experiments at quite different energies. In the proton-proton sector, the luminosity of the RHIC 
coUider should be enough to allow for measurements of DY processes. 

In a more distant future, COMPASS should providing new measurements of azimuthal asymmetries 
in SIDIS on a liquid hydrogen target and a the first measurement of polarised DY process in 7r~p^ 
scattering. Also, higher energy SIDIS data will come from JLab upgraded at 12 GeV. 

In an even more distant future, many projects have been proposed. The PAX experiment at FAIR 
aims to investigate DY pairs in polarised antiproton-polarised proton scattering, a very clean way to 
address the transvcrsity functions, and a terribly difficult experiment. Fixed target DY experiments 
scattering polarised protons on polarised protons are being planned at JPARC, in Japan, and at NIC A, 
in Dubna. It is fair to say, however, that the future of the field will depend in a crucial way on the 
pending decisions to construct a polarised electron-polarised proton collider. Ambitious projects are 
being pursued at BNL and at JLab. In Europe, the electron-proton collider has a long story, which needs 
not to be summarised here, we only mention the most recent proposal for a polarised collider, which is 
tailored to the new accelerators complex FAIR presently being realised at GSI. All these projects are 
based on existing laboratories and existing infrastructures so that optimism is mandatory. 
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